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Preface 


This  report  describes  the  development  of  a  real-time 
digital  simulation  of  the  motions  of  a  modular  escape 
capsule  released  from  a  submerged,  distressed  submarine. 
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under  Contract  Number  N00014-70-C-0247 . 
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I.  INTRODUCTION 


The  modern,  nuclear  submarine  is  one  of  the  most 
effective  weapon  systems  in  use  today  but  it  suffers  from 
one  major  flaw  —  in  case  of  disaster  the  crew  has  very 
little  chance  of  survival.  Future  submarines  will  be  able 
to  operate  deeper  and  deeper  and  hence  the  problem  cannot 
help  but  become  more  critical. 

In  the  aircraft  industry,  the  methods  of  crew  escape 
have  evolved  from  the  parachute,  to  the  ejection  seat,  to 
the  modular  ejection  capsule.  The  parallels  to  this  pat¬ 
tern  for  crew  escape  from  submarines  are  illustrated  below: 


Development  Stage 

Aircraft 

Submarine 

A.  Early  Development: 

no  escape  provision 

no  escape 
provision 

B.  Early  Operational 
Vehicles : 

Parachute 

Moms on  Lung 
(Free  Ascent} 

C.  Recent  Operational 
Vehicles: 

Ejection  Seat 

Steinke  Hood 

D.  Recent  Developments: 

Escape  Capsule 

? 

The  diving  bell  and  the  DSRV  have  no  parallels  in  air¬ 
craft  escape  and  are  not  applicable  unless  the  men  survive 
the  emergency.  This  requires  operation  in  water  depths 
less  than  the  crush  depth  of  the  vehicle.  In  addition, 
free  ascent  is  limited  to  relatively  shallow  depths  as  is 
illustrated  in  Figure  1-1. 


Figure  1-1.  Depth  Limits  for 
Various  Rescue  Techniques 
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Thus,  it  is  quite  logical  that  the  next  step  in  the 
development  of  escape  systems  for  distressed  submarines 
involves  a  direct  means  of  escape  utilizing  a  modular 
escape  capsule. 

The  interface  requirements  with  the  submarine  involve 
many  questions:  Can  the  submersible  be  controlled  from 
the  capsule?  What  is  the  arrangement  of  men  and  materials 
inside  the  capsule?  How  is  the  module  separated  from  the 
submersible?  What  are  the  static  and  dynamic  stabilities 
of  the  complete  submersible  and  of  the  separated  capsule? 

The  answers  to  these  questions  will  greatly  influence 
the  size  and  shape  of  the  escape  module;  but,  there  are 
also  considerations  in  the  ascent  to  the  surface  and  in  the 
survival  requirements  on  the  surface  which  will  influence 
the  size  and  shape  of  the  capsule.  An  operational  capsule 
will  therefore  represent  a  compromise  between  the  various 
requirements  for  submarine  interface,  ascent  requirements, 
and  surface  survival  requirements. 

In  the  present  study,  the  hydrodynamic  behavior  of  an 
escape  capsule  during  the  separation  and  ascent  phase  has 
been  analysed  and  simulated  in  a  real-time  digital  system 
called  MODSEC  -  Modular  Submarine  with  Escape  Capsule.  The 
MODSEC  system  can  simulate  the  behavior  of  a  submarine  with 
a  forward  escape  capsule  as  the  submarine  sustains  damage, 
the  capsule  separates,  ascends  to  the  surface,  and  floats 
on  the  surface  awaiting  rescue.  Figure  1-2  illustrates,  in 
a  condensed  fashion,  the  capabilities  of  the  simulation. 

The  MODSEC  system  has  been  designed  as  a  real-time 
system  utilizing  a  time-shared  computer.  A  Typagraph 
Model  3  terminal  is  used  as  input/output  device.  This 
terminal  allows  the  user  to  obtain  high  quality  graphical 
output  very  quickly.  The  MODSEC  system  incorporates  a  com¬ 
prehensive  command  language  which  enables  the  user  to 
change  input  parameters,  select  output,  and  analyze  various 
phases  of  the  capsule  motions  at  his  option.  Disc  file 
storage  is  incorporated  into  the  system  to  allow  the  user 
to  store  lengthy  weight  distributions  permanently.  This 
feature  prevents  costly  data  loss  due  to  computer  crashes, 
communications  failure,  or  user  mistakes.  The  MODSEC  system 
has  been  designed  so  that  personnel  untrained  in  the  com¬ 
puter  sciences  can  easily  use  the  system.  No  knowledge  of 
computer  languages  is  required  and  no  tedious,  complicated, 
data  entry  is  necessary. 


2 


Figure  1-2.  MODSEC  Simulation. 
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The  MODSEC  system  is  presently  a  sophisticated  inter¬ 
active  system  for  investigating  the  hydrodynamic  behavior 
of  truncated  ellipsoids  of  revolution.  Because  of  its 
modular  design,  MODSEC  has  great  potential  as  a  tool  for 
the  design  of  underwater  or  flying  bodies  of  virtually  any 
geometry  moving  in  up  to  6  degrees  of  freedom.  It  is  an 
easy  job  to  change  the  equations  of  motion  routine  or  the 
geometry  routine  to  accomodate  other  geometries  and  environ¬ 
ments  . 
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II.  PRELIMINARY  CONSIDERATIONS 


The  analysis  of  a  vehicle  having  six  major  degrees 
of  freedom  as  well  as  numerous  internal  degrees  of  free¬ 
dom  caused  by  elastic,  gyroscopic,  and  control  effects 
is  usually  approached  in  the  following  manner: 

(a)  The  vehicle  is  first  regarded  as  a  rigid  body 
and  the  equations  of  motion  are  derived  with  respect  to 
a  set  of  axes  fixed  in  it.  These  are  Euler's  equations. 

(b)  Modifications  required  by  gyroscopic,  elastic 
and  control  effects  are  evaluated. 

(c)  Due  to  moments  and  products  of  inertia  appear¬ 
ing  in  Euler's  equations,  a  choice  must  be  made  as  to 
whether  the  coordinate  system  be  earth  fixed,  in  which 
case  these  terms  become  functions  of  body  orientation; 
or  whether  body  axes  are  to  be  adopted,  which  avoids 
variable  moments  and  products  of  inertia,  but  introduces 
the  difficulty  that  a  complicated  coordinate  transforma¬ 
tion  must  be  introduced  to  relate  body  and  space  fixed 
axes.  The  latter  alternative  is  usually  adopted. 

(d)  The  resulting  trajectory  equations  are 
simultaneous  integro-dif ferential  equations,  and  can¬ 
not  be  conveniently  solved  without  extensive  linearization 
and  approximations. 

For  the  class  of  problems  under  consideration  here, 
however,  the  trajectory  analysis  is  considerably  simp¬ 
lified  since  the  motion  takes  place  in  a  plane,  and  the 
longitudinal  axes  of  the  vehicle  remains  in  the  plane. 

Thus  there  is  no  coordinate  variation  of  moment  of 
inertia  since  the  only  moment  of  inertia  relevant  to  the 
problem  is  that  which  lies  about  the  axis  normal  to  the 
plane  of  motion.  No  elastic,  gyroscopic  or  control  ef¬ 
fects  are  to  be  considered,  and  then  the  motion  may  be 
referred  to  a  single  set  of  earth  fixed  axes  by  a  simple 
transformation . 

This  is  not  to  say,  however,  that  a.  resulting 
equations  are  easily  solved:  for  the  e»-  nee  of  this 
problem  is  the  inclusion  of  large  angul  displacements, 
which  prohibits  linearizations  to  a  larg  extent.  Fur¬ 
ther,  the  effects  of  induced  maes  and  mo.nents  of  inertia 
must  be  included  and  finally,  the  modeling  of  the  various 
hydrodynamic  forces  and  moments  for  large  angular  dis¬ 
placements  requires  considerable  engineering  judgement. 
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III.  EQUATIONS  OF  MOTION 


Consider  a  rigid  body  (either  the  escape  capsule  or 
the  remaining  submarine)  moving  in  an  earth  fixed  x,  y 
coordinate  system.  It  can  be  shown  that  bodies  with 
axial  symmetry  moving  under  the  influence  of  hydrodynamic# 
hydrostatic#  gravitational#  and  inertial  forces  will  move 
in  a  single  plane.  The  problem  is  now  reduced  from  a 
c  general  6  degree  of  freedom  problem  to  one  with  3  degrees 
of  freedom.  The  body  motion  can  be  completely  described 
by  the  x#  y  coordinates  of  the  center  of  mass  and  the 
angle  Q  #  expressed  as  a  function  of  time.  See  Figure 
III-l. 


Figure  III-l.  Coordinate  System 


The  forces  on  the  body  roust  be  expressed  as  functions 
of  the  geometry  of  the  body#  its  velocity#  and  its  orientation. 
These  forces  are  of  three  types t  hydrodynamic#  hydrostatic# 
and  inertial. 

Referring  to  Figure  III-2  we  assume  that  the  lift 
force,  L#  and  the  drag  force,  D,  act  through  the  hydrodynamic 
center  h#  the  lift  force  perpendicular  to  the  velocity  vector 
V#  the  drag  force  parallel  to  V. 
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Resolving  the  lift  and  drag  into  their  x  and  y  components 
we  have 

=  -L  *»n  (©-<)  —15  ccri 


S’ 


L  co*  (.©  -°0  ~  t>  »iv\  Ce  ~«A 


where  F  and  F.  are  the  x  and  y  components  of  the  hydro- 
dynamic  force.”  The  moment  of  the  hydrodynamic  force 
about  the  center  of  the  mass,  g,  can  be  written  (counter 
cloc)cwi8e  moments  positive) 


jv-hJu  s\n  ©  —  d^Cose]*  S  [  (  dji ~  4  )  M4 

where  Md  is  the  hydrodynamic  damping  moment. 

It  can  be  shown  that  H4  -  K3  ©  where 

In  the  above  expression,  C*  is  the  crossflow  drag  co¬ 
efficient  for  a  circular  cylinder  and  <^{k)  is  the  equation 
of  the  ellipsoid  in  the  coordinate  system  of  Section  IV. 

The  hydrostatic  forces  include  the  buoyant  force 
acting  upward  on  the  vehicle  through  the  centroid  of  the 
underwater  volume,  b,  and  the  weight  of  the  vehicle  acting 
downward  through  the  center  of  mass  (See  Figure  III-3) . 


o 
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Figure  III-3.  Hydrostatic  and  Buoyant  Forces. 


Resolving  the  hydrostatic  forces  into  their  x  and  y 
components  we  get 

^  =  B-hr\C| 

Taking  moments  about  g  we  get 

Ms  ~  B  [  (  Cos©  *•  d^x  si*\e] 

The  inertial  forces  and  moment  are  proportional  to 
the  accelerations  and  can  be  written 

f  L%  =  *Wx\x 

f  iy  -  -  C  ™  •*  ™yy)  4 

Mt  «  *  (  I  +Iy)  S  .. 

where  »  Iv  are  the  virtual  m^sseo  and  moment  of  inertia 

The  3  equations  of  motion  can  now  be  written 


-Pi*  +  ‘  0 

f  s +  S  ♦  ^  '  0 

Mfc  V  -  M;_-0 


or 

.(ft)  +Tn,/x)  X  *  -L  *in  (.©-*}  -1)  cos  (©-<.} 

Cm  +  »v\rtj^  =  "B  —  Wi^  -»-  L  COS  (©-<*.)  -b$m  (&-d£) 

C^  +  Iv)^  s  BL(4gi'40  Cos  ©•  +  (PsmoC  +  Lsm»l)(4gr40 

-f  (Lsin«C-  t>cot.'*.N)  4  ~*  H<i 

Non-dimensionalizing  equations  1,  2f  and  3  by  non** 
diroensionalizing  all  lengths  by  the  body  length*  1*  and 
writing 


V* 

-  i/m 

t* 

3-  ^ 

“V^Zt  "t 

ft*- 

,  .^e» 

,  m* 

_  Znn 

5 

Pl$ 

T*  = 

ZI 

pp  s 

P= 

B/P4  _ 
AS 

pnsmd^c- 

CQeM’\c\e+& 

and  introducing  the  lift  and  drag  coefficients 

n  — — -  .  £  - — tL.  - 

■  >■  iyv's  J  *  tf»’s 

we  get  the  non-dimensional  equati one  of  motion  for  x*  * 
y*,  and  0. 


=- 

gtfl 


1"^.  =  ]*[(C»  ■'CLM8^)(d*-dt1) 


+(cus*n«t-  C*  Cof&-c)  4 *.^  2^  Ug~<0<^e  *4gz  s»*©J  -  Kj  0l0| 
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If  we  now  assume  that  the  virtual  masses  can  be  written 


Then 


m'*-  m*  ♦  2P 

3  * 

and 

where  kx  ,1^,  and  K$are  the  virtual  mass  and  inertia 
coefficients  in  the  x,  y,  and  0  directions.  These  co¬ 
efficients  are  discussed  in  more  detail  in  Section  V-2. 


The  first  step  in  the  solution  of  these  three  coupled, 
second  order,  ordinary,  differential  equations  is  to  obtain 
an  equivalent  set  of  six  first  order,  differential  equations. 
This  is  accomplished  as  follows t 

Let  __V(*  y* 

Y,  *  0 
V4  3 

Y„  = 

Then  define: 

F,  »  d  Y,  /  4  x  *  V4. 

,FI-  4Y*./  4k,  »  if 

Fy-  4vr/ax  •  a^/^x 

f.  «  a*#/4t,a’ 

w 
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The  preceding  six  equations  are  six  first  order  equations 

with  six  dependent  variables  (yi,  y2,  -  y§)  and  one 

independent  variable  (x) .  Initial  conditions  for  the 
six  dependent  variables  are  set  by  specifying  values  of 
¥1(0) ,  Y2  (0)  ,  —  yg(0) . 

We  use  a  fourth  order  Runge-Kutta  technique  to  solve 
simultaneously  the  six  first  order  differential  equations 
of  the  form 

#  »  ^»(  S>  i  MmN«i 

m 

The  algorithm  vised  to  calculate  y^(x+h)  given  ,y^(x)  is 

-  V  ■%(*  )  *>*,•••  ^ 

^i(  >*" ^Jn+ 'O 

v<+a  ^ 

where  h  is  the  non-dimensional  time  increment. 

O 
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IV.  SIMULATED  VEHICLE  GEOMETRY  AND  WEIGHT  DISTRIBUTION 


A  thorough  investigation  of  representative  sub¬ 
marine  hull  shapes  revealed  that  the  intact  vehicle  can 
be  mathematically  represented  by  an  ellipsoidal  body  of 
revolution.  This  selection  provides  a  reasonably  ac¬ 
curate  geometric  model  which  is  relatively  simple  to 
manipulate  mathematically.  Although  some  inaccuracies 
are  introduced  in  the  after  portion  of  the  submarine# 
the  study  is  primarily  concerned  with  the  forward  section 
—  the  capsule  —  which  is  very  closely  approximated  by  a 
portion  of  an  ellipsoidal  body  of  revolution.  Figure  IV-1 
illustrates  the  plan  view#  with  the  appropriate  axes. 


V 


Figure  IV-1.  Ellipsoidal  Body  of  Revolution 
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Vehicle  Geometry 


It  is  assumed  that  the  major  axis  of  the  ellipsoid 
is  the  axis  of  revolution  and  that  this  corresponds  to 
the  longitudinal  reference  axis  of  the  vehicle.  The 
equation  of  an  ellipse  in  the  7,  7  axis  system  (Fig*.  IV-1)  is 


where  a  is  the  semi-minor  axis,  and  b  is  the  semi -major 
axis.  Noting  that 

X  »  X.  -v  b 

-  V  ^ 

the  equation  of  the  ellipse  in  the  x,  y  coordinate  .system 
is 

L^lL  +  .Hi  =i  u> 

b-1  ax 

The  geometric  quantities  that  are  of  interest  in  the 
analysis  of  the  dynamics  of  the  vehicle  are: 

volume , 
surface  area, 
volume  moment  of  inertia, 
radius  of  gyration, 

longitudinal  center  of  gravity  of  the  vehicle  shell, 
longitudinal  center  of  buoyancy. 

The  geometry  subroutine  calculates  each  of  these  parameters 
for  the  intact  vehicle,  the  capsule,  and  the  remaining 
vehicle . 


Volume 

The  volume  of  any  segment  of  the  ellipsoid  of  revolu¬ 
tion,  say  between  stations  x»  and  xj,  is  obtained  from 

f*1  ( 

Solving  Eq.  (1)  for  y2 i 
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(3) 


Substituting  this  into  Eq.  (2) : 

v-  (  irf»- atjA*. 

%  «•  b1  -» 

Integration  yields : 

where  the  length-to-diameter  ratio  was  taken  as  L/D  =  b/a. 
Non-dimensionalizing  the  volume  with  respect  to  the  length: 


r»  $  • 


where  _  X4*  -  V(x/t> 

X,4  -  *7  b 

Surface  Area 


The  surface  area  of  any  segment  of  an  ellipsoid  of 
revolution  between  stations  xi  and  x*  is  given  by: 

A* 

X, 

Differentiating  Equation  (3) : 


Ju  _  /j.\  (b~*^ 

dy.  bbj  (_zb*-v^Kt 

Prom  Eqs.  (3)  and  ( 4 ) : 

A  -  23Lpt  [  |D 


(4) 


(5) 
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Integrating  Eq.  (5)1 


*  ■* 


Noting  that  as  a  -*■  b  (or  a/b  1)  ,  the  second  term 
involving  the  Arcsin  becomes  indeterminate ,  we  use  the 
fact  that  as  x-to,  e JX _ 

Thus  from  Equation  (5a)  ,  as  _a/b_=*_i_} 


7  * 

4-  bl4x-b)  f 

^  ^  J  X, 

Returning  to  the  full  equation  and  non-dimensionalizing. 


il.-ltH*  L  H  UM  , 

Volume  Moment  of  Inertia  *i 

The  volume  moment  of  inertia  (I)  of  any  segment  of  the 
ellipsoid  of  revolution  about  the  point  jc«0  is  given  by: 

_x_=  Loxitis _ 

where  j3  is  the  density  of  the  ellipsoid.  The  differential 
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volume  of  a  slice  of  the  ellipsoid  is 

and,  hence,  the  moment  of  inertia:  '  •* 

1=  S  J[Nyx*  * 

Using  Eq.  (3)  allows  writing  this  as 


X  s  lfo> 

f 


Integration  yields: 


which,  upon  non-dimensionalizing  with  respect  to  the 
vehicle  length  (2b)  to  the  fifth  power,  yields: 


Using  the  parallel  -  axis  theorem  allows  calculation  of 
the  moment  of  inertia  about  any  axis  if  the  distance 
from  x=0  to  the  axis  is  known;  i.e.. 


3-2*  ^C<j  +  (  m &.*>)(,  diiUwte  Ioc^wc^k  Ajc-es^2, 

As  an  example,  the  moment  of  inertia  about  the  center 
of  gravity  is 


or 


_r  _ 
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(6) 


Radius  of  Gyration 

The  radius  of  gyration (fg)  is  given  by 

V  i - 

«  m 


(7) 
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Thus,  from  the  previous  results: 


The  value  may  be  obtained  from  previously  calculated  re¬ 
sults.  The  radius  of  gyration  may  also  be  transferred  to 
any  other  axis  using  the  parallel-axis  system.  From 
Eqs.  (6)  and  ( '7  }  : 


Longitudinal  Center  of  Bouyancy 

The  longitudinal  center  of  buoyancy  of  any  segment 
of  the  ellipsoid  of  revolution  is  given  by  the  first 
moment  of  volume  about  the  point  x»0 ,  divided  by  the 
volume: 


■tee 


,  L  'E'x'jMx. 


Using  Eq.  (3),  the  integral  in  the  numerator  of  Eq.  (8) 
may  be  written: 

.>(l 


iTy*  ^  x[  i-  <*<*• 

1 


Integration  yields : 


k  ?  L  3!o  4- Id*  ] 

Thus,  the  longitudinal  center  of  buoyancy  of  any  segment 
of  the  ellipse: 


_ _ 

^[■3  b1  +£c  \-W } 


(8) 


Non-dimensionalizlng  with  respect  £.0  the  length  (2b)  : 


2.b 


][-*  us-*:)  o<i-o’] 


Centroid  of  the  Projected  Area 


The  centroid  of  the  projected  area  of  a  segment  of  an 
ellipsoid  of  revolution  cam  be  written 


With  the  help  of  Equation  0 )  this  cam  be  written 


Carrying  out  the  integrations  in  both  the  numerator  and 
denominator  and  non-dimensionalizing  as  before,  yields 
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Longitudinal  Center  of  Gravity  of  the  Shell 


The  longitudinal  center  of  gravity  of  the  capsule 
and  remaining  vehicle  shell  is  required  for  the  sub¬ 
sequent  determination  of  the  center  of  gravity  of  the 
intact  vehicle,  the  capsule,  and/or  the  remaining  vehicle. 
To  do  this  the  total  weight  of  the  vehicle  shell  must 
first  be  calculated.  We  assume  that  the  capsule  and  the 
remaining  vehicle  have  constant,  but  not  necessarily 
equal,  shell  thickness  and  may  even  be  constructed  of 
different  materials.  Since  the  shell  thickness  is  small 
compared  to  the  average  radius  of  curvature,  the  weight 
of  the  capsule  shell  may  be  closely  approximated  by 
the  surface  area  of  the  shell  timesthe  shell  thickness 
times  the  specific  weight  of  the  shell  material,  and, 
similarly,  for  the  remaining  vehicle.  The  total  shell 
weight  for  the  intact  vehicle  is  then  the  sum  of  these 
two  weights.  Note  that  these  are  the  ellipsoidal  shell 
weights  and  do  not  include  the  weight  of  any  end  plate 
required  to  close  off  the  capsule  and  remaining  vehicle 
from  the  surrounding  environment  upon  separation. 


The  longitudinal  center  of  gravity  of  any  segment 
of  the  ellipsoidal  shell  can  now  be  written  as 


C.  d»  =■ 
'UlaeU 


From  the  above  discussion,  the  total  weight  of  the  shell 
is  toW0S.  Thus,  the  center  of  gravity  becomes i 
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& 


Non-diraensionalized  with  respect  to  the  length  (2b)  this 
yields: 

£a,wi .  i  { (E)*+  ?[>-  (?fj  x*- 


#ol 


Weight  Distribution 


X," 


The  weight  of  the  intact  vehicle,  capsule,  or  remain¬ 
ing  vehicle  is  composed  of  the  weights  of  the  individual 
components,  the  shell  weight,  and  the  weight  of  the  closure 
plate.  The  weights  of  the  individual  components  are  assumed 
to  be  known.  The  weight  of  the  ellipsoidal  shell  is  assumed 
to  be  the  surface  area  of  the  shell  times  the  appropriate 
thickness  times  the  specific  weight  of  the  appropriate 
material.  The  closure  plates  are  assumed  to  be  flat  plates 
of  the  same  thickness  and  material  as  the  capsule  or  re¬ 
maining  vehicle  shell.  The  weight  of  the  closure  plate 
is  then  the  area  of  the  plate  times the  thickness  times  the 
specific  weight  of  the  appropriate  material. 

Mathematically, 

W-M.i  •  £  *4  --'‘4m  +  '•Jfl.+c 

where  the  are  the  individual  component  weights. 

The  longitudinal  center  of  gravity  is  given  by  the 
summation  of  weights  times  the  moment  arm,  divided  by 
the  total  weight: 

_£e£,  *■  v  Xca»M  »  WpUfe  X>c» 

Similarly,  the  vertical  or  radial  center  of  gravity: 

m,  .  "4  'a 

It  should  be  noted  that  there  is  no  contribution  to  the 
vertical  c.g.  from  either  the  shell  or  the  closure  plate. 

The  subroutine  WGTS  calculates  the  weights  of  the 
intact  vehicle,  the  capsule,  and  the  remaining  hull,  and 
their  respective  logitudinal  and  vertical  centers  of 
gravity. 
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V.  HYDRODYNAMIC  COEFFICIENTS 


Accurate  prediction  of  the  hydrodynamic  coefficients 
is  a  requirement  for  an  accurate  dynamic  simulation  of 
the  motions  of  any  underwater  vehicle.  It  is  unfortunate 
that  there  is  such  a  paucity  of  experimental  data  and/or 
theoretical  results  which  can  be  used  in  studies  of  the 
type  undertaken  here.  The  vehicles  which  are  being  model¬ 
led  are  truncated  ellipsoids  of  revolution  (see  Section 
IV)  with  L/D  ratios  between  1  and  10.  We  assume  that  the 
intact  submarine  is  an  ellipsoid  of  revolution  and  thus 
both  the  separated  capsule  and  the  remaining  vehicle  are 
truncated  ellipsoids.  Both  the  capsule  and  the  remaining 
vehicle  may  operate  at  very  large  angles  of  attack,  from 
0  to  360  degrees. 

Because  of  the  lack  of  hydrodynamic  data  on  truncated 
ellipsoids  of  revolution,  considerable  engineering  judge¬ 
ment  has  been  necessary  in  order  to  formulate  the  equations 
of  motion.  Where  good  data  was  available,  it  was  used. 

Where  the  data  is  questionable,  simplified  estimates  were 
used. 


V— 1  The  Lift  and  Drag  Coefficients 

The  resistance  of  the  body  may  be  thought  of  as 
composed  of  three  parts: 

a.  Skin  Friction 

b.  Pressure  Drag 

c.  Base  Drag 

Assuming  that  the  flow  is  turbulent  the  two-dimen¬ 
sional  skin  friction  coefficient  is  given  by 


where 

logarithm. 


and 


— T 


3.4U  \r\iZ cV 

.  lb  denotes  the  natural 

-UvlS 

is  the  Reynold's  number 


based  on  the  length  cf  the  body.  Skin  friction  coefficients 
for  three-dimensional  bodies  are  larger  because  the  3-D 
relief  effect  thins  the  boundary  layer  and  hence  increases 
the  shear  stress  at  the  wall.  For  turbulent  flow  the  in¬ 
crease  in  due  to  three-dimensional  effects  can  be  written 
(See  Reference  1) 


o.ooii,  %  [e*7"VsL 
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The  pressure  drag  for  ellipsoidal  hodie  operating  at 
Reynolds  numbers  above  7.5  X  10s  can  be  ^pressed  as  a 
polynomial  in  d/£  (Reference  1) 

0.U97  02.MJ  4-0.l39t.SI  l  4) 2 
If  we  define  q,  &  Cf  -7  where  SB  is  the  base  area/ 


then  the  base  drag  for  3-D  bodies  can  be  written  (Reference 
1) 

<>  *  -2-  <*•».  . 

M>«  r- 

VCc7 

For  bodies  with  a  base  diameter  smaller  than  the 
diameter  of  the  maximum  cross-section  we  have 


o  .cn°> 


(itS 


'•TcZZ 


The  total  drag  at  zero  angle  of  attack  is  then 


-A 

where  ' 

A  =  &  *  Q,  *•  e„s  . 

The  lift  coefficient  may  initially  be  approximated 
by  the  cross  flow  principle  for  circular  cylinder*'  in 
laminar  flow  (see  Figure  V-l) s 

.  CL*  C*  SmV  £©*«<.  - _ — - 

C*  ^ 

where  Cfi  is  the  side  force  coefficient  based  on  the  pro¬ 
jected  area 


•gVg.  -fort-c 


w,  z 
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Figure  V-l 


We  can  now  write  a  relationship  for  the  lift  and 
drag  of  the  capsule  which  depends  on  the  length-diameter 
ratio  and  the  angle  of  attack: 

—  -  C5  (  SIW2oi  0 ... 

~  £j>|  +■  Cfc  cosV. 

The  above  lift  and  drag  coefficients  are  admittedly 
rough  approximations  but  better  estimates  would  require 
that  a  comprehensive  experimental  program  be  undertaken. 


V-2  Added  Mass  Coefficients 


When  either  the  escape  capsule  or  the  remaining  vehicle 
is  subjected  to  external  forces  and/or  moments  (buoyancy, 
hydrodynamic,  etc.)  not  only  must  the  mass  of  the  vehicle 
itself  be  accelerated  but  also  the  mass  of  seme  of  the 
surrounding  water.  Newton's  second  law  appears  as 

_  ?  -  _ 

for  linear  motions; 

H  -  (i-iv'ie 

for  angular  motions;  where  my  and  Iy  are  termed  the  virtual 
moment  of  inertia,  respectively.  In  the  notation  of  Section 
III,,  Newton's  second  law  becomes  for  motion  in  the  x,  y  plane 


(Vl  ♦  x  »  F* 

9  =  s 

where  M  and  M  are  the  virtual  masses  in  the  x  and  y 
directions.  The  virtual  mass  is  expected  to  be  proprotional 
to  the  displaced  mass  of  fluid 

^vx  ’  <)  ~  y ^ 

where  K  and  K  are  the  virtual  mass  coefficients.  Similarly, 
the  virtual  moment  of  inertia  is  assumed  to  be  proprotional 
to  the  moment  of  inertia  of  the  displaced  mass  about  an 
axis  through  the  center  of  buoyancy;  hence 

-Tv=  f^'ke 

where  rg  is  the  radius  of  gyration  of  the  displaced  volume 
and  is  the  virtual  moment  of  inertia  coefficient. 

Ke  will  not  depend  on  the  body  orientation,  but  K 

and  X  definitely  will.  We  will  assume  in  this  study 

that  X  and  K  can  be  written 
x  y 

k*  =  cos  <|  +-  Kc  1  . . 

-  «  14^1  sinotj  +  1  coso<.i 

where  K  and  K  are  the  axial  and  cross-flow  added 
mass  coefficients.  c 

To  the  best  of  our  knowledge,  values  of  K  and  K 
for  truncated  ellipsoids  of  revolution  are  noS  available 
either  in  the  experimental  or  theoretical  literature. 

Even  if  they  were,  the  results  would  be  of  questionable 
value  in  this  work  because  the  motions  of  the  capsule 
are,  for  the  most  part,  nonosci 11a tory  and  nearly  all  of 
the  published  data  is  for  small  amplitude,  oscillatory 
motions  where  potential  flow  can  be  assumed.  The  virtual 
mass  coefficient  of  a  blunt  body  when  the  flow  is  separated 
is  much  harder  to  analyze  and  measure  because  it  not  only 
depends  on  the  shape  of  the  body  but  also  on  the  velocity 
and  acceleration  in  a  complicated  way.  Also,  because  of 
viscous  effects,  the  virtual  mass  not  only  depends  on 
the  instantaneous  acceleration  and  velocity  but  also  on 
the  motion-time  history.  Despite  the  difficulties  in 
estimating  the  true  virtual  masses  it  appears  that  the 
potential  flow  values  do  give  reasonable  results,  at  least 
for  the  initial  acceleration  from  rest.  However,  even 
potential  flow  virtual  mass  coefficients  are  not  available 
in  the  literature  for  truncated  ellipsoids  of  revolution. 

Lamb  (Reference  2}  derives  the  K  ,  K  ,  and  K9  for  intact 
ellipsoids  of  revolution.  Thesexresults  follow: 
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-Kc>  - 


U-  ex)  L  z<L2-  fe  -ex)Gfo-4oYj 


— <£  «  ecceiA  ^nc\4»j  a  iv  ~a-*  =^t_  ^B.y. 

-x*' 

f-  ■  -fc  - 

and  a  and  b  are  the  semi-minor  and  semi-major  axes 
of  the  ellispoid,  respectively.  As  an  initial  approximation, 
the  virtual  mass  coefficients  of  a  truncated  ellipsoid 
are  assumed  to  be  the  same  as  the  coefficients  of  a  whole 
ellipsoid  of  revolution  with  the  same  L/D  ratio. 

This  assumption  of  course,  neglects  the  effect 
of  the  truncation,  but  hopefully  includes  the  effect  of  L/D 
on  the  virtual  mass  coefficients.  Figure  V-2  shows  how  K  , 


Kc,  and  vary  with  L/D. 


Hydrostatic  Parameters 


V-2  shows  how  K  , 

A 


The  motion  of  the  capsule  will  be  influenced  to  a 
large  extent  by  the  buoyant  forces  and  moments.  The  buoyant 
force,  B,  acts  upward  through  the  centroid  of  the  dis¬ 
placed  volume,  b.  The  weight,  W,  acts  downward  through 
the  center  of  gravity,  g.  The  magnitude  of  the  ratio 

C'  IT  '  7v- 

determines  whether  the  capsule  will  rise  or  sink.  If 
C,  <  l  the  capsule  will  rise.  If  C,>i  the 

capsule  will  sink.  If  C,  »i  the  capsule  is  neutrally 
buoyant  and  will  not  move  vertically  except  under  the  in¬ 
fluence  of  hydrodynamic  forces.  The  hydrostatic  moment 


-M**  cos©  -v 
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.  a  function  of  e  VS’fSS*  “ 

ertically.  This  angle  denoted  hy 


Figure  V-3 


is  the  angle  which  the  capsule  would  assume  if  the  hydro- 
dynamic  forces  were  zero  and  it  was  free  to  trim.  The 
capsule  should  assume  an  angle  close  to  after  the 
angular  oscillations  have  damped  out  and  the  capsule  is 
in  steady  motion  long  after  it  has  separated  from  the  sub 
marine. 
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VI.  SIMULATION  OF  THE  MOTIONS  OF  THE  CAPSULE  ON  THE  SEA  SURFACE 


The  attitude  of  the  capsule  when  it  reaches  the  sur¬ 
face  will,  together  with  the  sea  state,  determine  the 
magnitude  of  the  heave  and  pitch  motions  experienced  by 
the  capsule.  The  first  task  is  to  determine  the  attitude 
the  capsule  will  assume  when  floating  on  the  surface  waiting 
for  rescue.  It  will  rise  until  the  displacement  is  just 
equal  to  the  weight  and  will  trim  until  the  centroid  of  the 
displaced  volume  is  vertically  colinear  with  the  center 
of  gravity  and  the  metacentric  height  is  positive  (See 
Figure  VI- 1) .  Since  there  may  be  more  than  one  stable 
attitude,  we  will  assume  that  eventually  the  capsule  will 
attain  the  most  stable  orientation,  that  is,  the  orientation 
which  yields  the  largest  metacentric  height. 

i 


Figure  VI-1.  Capsule  Floating  on  the  Surface 


Having  reached  equilibrium  on  the  surface,  it  will  be 
excited  in  heave  and  pitch  by  the  waves.  In  order  to 
simplify  the  difficult  problem  of  predicting  the  heave 
and  pitch  for  arbitrary  capsule  attitudes,  several  assumptions 
will  be.  made: 

a.  The  capsule  is  assumed  to  be  slender  enough  so 
that  it  will  be  stable  in  one  of  three  orientations: 
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1.  vertical,  bow  up 

2.  vertical,  bow  down 

3.  nearly  horizontal 

In  other  words  we  are  assuming  that  the  capsule  will 
behave  similarly  to  a  spar  buoy. 

b.  If  the  capsule  is  most  stable  vertically  we  assume 
that  the  pitch  angle  is  always  zero  and  only  heaving 
motion  is  significant. 

c.  If  the  capsule  is  most  stable  horizontally  we 
assume  that  the  pitch  angle  is  equal  to  the  wave  slope 
and  the  heave  amplitude  is  equal  to  the  wave  amplitude. 

With  the  above  assumptions  the  problem  becomes  one 
of  determining  the  most  stable  attitude  of  the  three  pos¬ 
sible,  and  then  calculating  heave  and  pitch  amplitude 
spectrums  from  the  specified  sea  state. 


The  Most  Stable  Attitude 

We  imagine  the  capsule  to  be  enclosed  in  a  rectangular 
box  as  in  Figure  VI-2. 


Figure  VI-2.  Digital  model  of  capsule. 
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The  box  is  now  divided  into  10 6  small  rectangular 
cells,  each  having  dimensions  SL/ioo*  rb/\oo  x  rfe/joo 
Wherever  the  capsule  occupies  any  part  of  a  cell  it  is 
considered  filled  with  material.  We  now  have  a  repre¬ 
sentation  of  the  capsule  volume  which  is  not  exact  but 
is  sufficiently  accurate  for  our  purposes  and  is  very 
simple  to  manipulate  in  the  computer. 

The  capsule  can  be  represented  by  a  three-dimensional 
matrix  V(l,l,l)...  V (100 , 100 , 100)  the  elements  of 
which  are  1  if  the  capsule  occupies  the  cell  and  0  other¬ 
wise.  The  matrix  is  arranged  so  that  the  first  index 
represents  locations  on  the  x  axis,  the  second  index  repre¬ 
sents  the  y  axis  and  the  third  the  z  axis.  For  example, 
v(l,  50,  50)  is  the  cell  at  the  center  of  the  separation 
plane;  V(100,  50,  50)  is  the  cell  at  the  nose  of  the 
capsule. 

The  equilibrium  draft  of  the  capsule  and  other  hydro¬ 
static  parameters  necessary  to  compute  the  metacentric 
height  at  each  of  the  three  allowable  attitudes  are  easy 
to  calculate  once  the  digital  capsule  model  has  been 
generated. 


The  draft  is  found  by  counting  the  occupied  cells 
starting  at  a  base  plane  which  is  parallel  to  the  water- 
plane.  Occupied  cells  are  counted  layer  by  layer  until 
the  number  of  cells  multiplied  by  a  scale  factor  just 
exceeds  the  weight  of  the  capsule.  The  waterplane  is 
then  located  at  the  layer  of  cells  at  which  the  displace¬ 
ment  most  nearly  equals  the  weight.  At  the  same  time  the 
equilibrium  drafts  are  found,  the  first  moments  of  the 
three  underwater  volumes  about  the  baseplane  parallel  to 
the  waterplane  can  be  calculated  by  multiplying  each 
element  by  its  distance  from  the  appropriate  baseplane. 


The  drafts  for  the  three  attitudes— horizontal 
(d,  ) ,  vertical  bow  down  (d  .) ,  vertical  bow  up  (d  _) 
can  be  written 
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where  nvi,n,  ,  and  n  ,  are  found  from  the  following  ex¬ 
pressions: 
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The  vertical  moments  M»,  M  . ,  and  My-  can  be  written 
in  a  similar  fashion  once  n”,  n^£,  and  n^2  are  known: 
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The  heights  of  the  centroid  from  the  appropriate  base- 
plane  can  now  be  computed: 
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The  metacentric  radii  of  each  waterplane  can  be 


written 
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And,  finally,  the  roetacentric  height  for  the  capsule 
floating  at  each  of  the  three  attitudes  can  be  written: 

=  +'oK-^v,  - - 

2mVi= 


where  g.  ,  a  , ,  9v,  are  the  distances  of  the  center  of 
gravity  r rom  the  ^appropriate  base  line  planes. 

The  most  stable  attitude  and  hence  the  most  probable 
attitude  is  the  attitude  which  yields  the  largest  . 


The  Pierson-Moskowitz  Spectrum  is  widely  accepted  by 
naval  architects  since  its  introduction  about  1965  as  a 
realistic  representation  of  fully  developed, 
deep-water  irregular  waves.  The  spectrum  is  based  on  the 
analysis  of  several  hundred  sea  wave  records  recorded  by 
British  scientists.  Mathematically,  the  spectrum  can  be 
written 

-6/u)* 


where 
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S(W) 

W 

A 

B 
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=  Spectrum  ordinate  (ft  -sec) 

=  Circular  wave  frequency  (rid/sec) 
=  .0081  g2 


=  33.56/(5 

=  Acceleration  of  gravity  _ 

=  Significant  wave  height  =*4fE, 
=  Cs(W)  dw  <*;  energy  in  the 


sea  state 


To  obtain  the  spectrum  as  a  function  of  the  sea  state, 
we  use  the  following  tabular  relationship  between  L  and 
the  sea  state:  ^ 


Sea  State 
0 
1 
2 

3 

4 

5 

6 

7 

8 
9 


0.10-0.15 

0.50-1.20 

1.50- 3.00 

3.50- 5.00 
6.00-7.50 
8.00-12.00 

14.00-20.00 

25.00-40.00 

45.00-60.00 

70.00-100.00 


The  maximum  wave  slope  of  deep  water  waves  can  be  written 


where  ^  is  the  wave  amplitude.  But  for  deep  water  waves, 
X-  Zir^/oo1  .  Therefore  /<% 

and  the  ordinates  of  the  slope  spectrum  can  be  obtained 
from  the  amplitude  spectrum: 

S  M=  s  Cu)") 

We  arc  now  ready  to  calculate  the  pitch  and  heave 
response  of  the  capsule,  since  the  attitude  is  known.  If 
the  capsule  is  stable  in  the  upright  attitude  then  the  pitch 
angle  can  be  assumed  to  be  zero  and  the  heave  amplitude  can 
be  written  t 

k'1-/  \  (c-cuox)V  . 

=  — - — — -r— — _e  q  s _ 

(c  -  w  to  -  <uo^Sb\ol 


where 


C  =  pg  A  waterplane 
a  =  added  mass  a  .5m 
b  =  damping  coefficient 
m  =  mass  of  capsule 
d  =  draft  of  capsule 
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If  the  capsule  is  roost  stable  in  the  horizontal  attitude, 
the  pitch  angle  is  assuroed  to  be  equal  to  the  wave  slope 
and  the  heave  amplitude  is  assumed  to  be  in  phase  and 
equal  to  the  wave  amplitude.  With  these  assumptions  we 
can  write 


(co)  *  5  1  Clo) 

Cio)  -  5  Cuo) 


and 


4.0 


6 


A|  *  4.0 

T  . 


M  Auo  . 
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VII.  TEE  DYNAMIC  SIMULATION 


One  of  the  primary  objectives  of  this  study  is  to  de¬ 
velop  a  dynamic  simulation  of  a  modular  submarine  which 
would  allow  the  researcher  to  study  in  a  realistic  fashion 
the  effects  of  changes  in  the  physical  configuration  of  a 
modular  submarine  on  the  dynamics  of  the  separation  and 
ascent  of  the  forward  capsule. 

The  simulation  was  required  to  have  several  important 
features : 

(a)  The  simulation  must  be  capable  of  handling 
arbitraiy vehicle  weight  distributions. 

(b)  It  must  be  able  to  handle  various  vehicle 
geometries. 

(c)  Various  types  of  casualties  must  be  able  to  be 
simulated. 

(d)  It  must  be  user-oriented. 

(e)  It  must  operate  in  as  close  to  a  real-time 
environment  as  possible. 

To  satisfy  the  above  requirements  a  digital  simula¬ 
tion  model  was  developed  for  a  time-shared  PDP-10  com¬ 
puter.  Input  and  output  is  provided  via  a  Typagraph 
interactive  terminal.  The  simulation  model  takes  ad¬ 
vantage  of  the  Typagraph' s  graphical  capability  and 
provides  graphical  output  of  either  the  capsule  or  the 
remaining  vehicle  motion  after  separation  occurs. 

The  simulation  model  appears  to  the  user  as  simply 
a  list  of  commands  which  he  types  on  the  terminal.  When 
a  command  is  typed  the  computer  interprets  the  command 
and  performs  certain  calculations  or  input/output.  After 
a  command  has  been  executed  the  computer  types  two 
asterisks  (**)  to  indicate  that  it  is  ready  for  another 
command.  The  simulation  is  started  by  typing 

RUN  DSKsMuOSEC^ 

Note:  |  is  a  carriage  return  and  is  typed  after  each  line. 

The  following  section  describes  the  simulation  com¬ 
mands.  They  are  listed  in  the  approximate  order  that  they 
might  be  used  in  a  sample  simulation  (see  Section  VII 
for  some  sample  runs) .  Each  command  may  be  abbreviated  by 
the  underlined  letters. 


( 


Command 


Description 


( 


**HELP}  Prints  the  contents  of  the  data  file 

HELP.DAT  which  contains  a  list  of  com¬ 
mands  with  brief  descriptions  as  well 
as  hints  for  running  the  simulation. 

Also  included  in  the  HELP . DAT  file  is 
a  list  of  modifications. 

** INPUT  XXXX}  The  input  commands  are  two- word  com¬ 

mands  used  to  input  data  for  a  simu¬ 
lation.  The  two  words  must  be  separated 
by  a  blank  and  each  word  can  usually 
be  abbreviated  by  its  first  letter. 

All  input  commands  may  be  repeated  in¬ 
dividually  during  a  simulation.  This 
feature  gives  the  user  a  convenient 
method  of  updating  or  changing  single 
input  parameters  without  having  to 
enter  all  of  the  variables  each  time. 

**INPUT  GEOMETRY }  The  computer  asks  for  the  L/D  ratio, 

overall  length,  and  location  of  the 
separation  bulkhead.  The  L/D  ratio  is 
the  L/D  of  the  intact  submarine,  the 
length  is  in  feet,  and  the  location 
of  the  separation  bulkhead  is  given 
as  a  decimal  fraction  of  the  overall 
length.  That  is,  a  number  between  0 
and  1.0. 


**INPUT  WEIGHTS}  Allows  the  user  to  either  type  in  a  new 

weight  distribution  consisting  of  a 
series  of  component  weights,  together 
with  the  horizontal  and  vertical 
locations  of  each  weight  or  retrieve 
a  previously  entered  weight  distri¬ 
bution  from  a  disc  file. 


**^NPUT  VOLUMES}  Allows  the  user  to  type  in  a  list  of 

floodable  spaces  which  will  be  used 
later  to  simulate  a  casualty.  Casualties 
are  simulated  by  flooding  user-specified 
volumes.  As  with  the  weight  distribution 
the  option  of  retrieving  a  previously 
stored  volume  distribution  from  a  disc 
file  is  provided. 
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fcsfc,;  si*#1  wf**'  “  '*• 


“UPDATE* 


*  ‘INPUT  DEPTHS 


“INPUT  SPEED* 


“ INPUT  TRIM* 


* ‘WEIGHTS* 


“INPUT  SEAS* 


“DAMAGE* 


“SEPARATE* 


This  command  gives  the  user  the  capa¬ 
bility  to  update,  change,  and  add  to 
the  current  weight  or  volume  distribu¬ 
tion.  This  command  prompts  the  computer 
to  ash  a  series  of  questions  regarding 
the  exact  nature  of  the  updating  to  be 
done;  these  questions  are  self  explanatory. 

See  the  example  problems. 

Allows  the  user  ..to  type  in  the  depth  at 
which  the  casualty  and  separation  is  to 
occur.  This  number  is  not  used  in  any 
calculations . 

Allows  the  user  to  type  in  the  initial 
horizontal  velocity  and  vertical  velocity 
in  ft/sec.  This  is  the  speed  of  the 
submarine  at  the  instant  of  separation. 

Allows  the  user  to  type  in  the  trim  angle 
of  the  submarine  at  the  instant  of 
separation.  The  command  TRIM  should  be  used 
to  guide  the  user  in  his  choice  of  initial 
trim  angle. 

Calculates  the  various  vehicles'  weights 
and  the  longitudinal  and  vertical  centers 
of  gravity.  This  command  is  given  fol¬ 
lowing  INPUT  WEIGHTS  or  in  the  case  of 
a  damaged  vehicle  DAMAGE,  but  always  prior 
to  ASCEND  or  HULL. 

Allows  the  user  to  type  in  the  surface  sea 
state  as  an  integer  from  0  to  9. 

This  value  is  used  to  calculate  a  sea 
spectrum  which  in  turn  is  used  to  calcu¬ 
late  the  pitch  and  heave  response  of  the 
capsule  on  the  surface. 

This  command  allows  the  user  to  specify 
which  of  the  spaces  he  wishes  to  flood. 

This  command  simulates  various  casualty 
types  by  flooding  the  specified  spaces. 

Separates  the  capsule  from  the  submarine 
and  initializes  parameters  for  the  ascent 
phase.  This  command  does  not  produce 
any  output  or  require  any  Input . 
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**TRIM) 

**ASCEND> 

**HULL) 

**GRAPH) 

**MOTIONS-> 

**PRINT> 


Calculates  the  net  buoyant  force  on 
the  intact  vehicle  and  the  trim  angle 
the  intact  vehicle  would  take  if  it 
were  allowed  to  reach  equilibrium. 

This  command  uses  the  intact  vehicle 
in  its  current  condition  (either 
damaged  or  undamaged)  for  the  calcu¬ 
lations  and  prints  out  the  results  on 
the  terminal. 

Calculates  the  trajectory,  angular 
orientation,  horizontal  velocity  of 
the  capsule  as  a  function  of  time. 
t=0  is  assumed  to  be  the  time  that 
separation  occurs.  Before  the  cal¬ 
culation  is  done  the  user  must  input 
the  maximum  time  in  seconds  he  wishes 
to  calculate  data. 

Essentially  the  same  as  "ASCEND"  ex¬ 
cept  that  the  trajectory  of  the  re¬ 
maining  hull  is  calculated. 

Generates  an  annotated  plot  on  a  Typa- 
graph  Model  3  terminal  of  the  horizontal 
and  vertical  distance  from  the  separation 
point  and  the  angular  orientation  of  the 
vehicle  as  a  function  of  time.  This  com¬ 
mand  plots  the  most  recently  generated 
trajectory.  For  example,  if  "HULL"  was 
typed  just  prior  to  "GRAPH"  the  data 
plotted  would  be  for  the  remaining 
vehicle.  If  "ASCEND"  was  typed  im¬ 
mediately  prior  to  "GRAPH"  the  data 
plotted  would  be  for  the  capsule. 

Calculates  and  plots  the  behavior  of 
the  capsule  at  the  surface.  Both  the 
pitch  and  heave  response  spectrums  are 
plotted  together  with  the  capsule 
orientation,  the  average  wave  height, 
heave  amplitude  and  pitch  angle. 

Prints,  in  tabular  form,  the  dimensional 
and  non-di.mensional  results  of  the  last 
calculated  trajectory.  This  command  is 
used  primarily  when  running  the  simu¬ 
lation  from  an  ordinary  Model  33  tele¬ 
type.  Before  printing  the  data  the 
computer  asks  how  many  points  the  user 
wishes  to  print  out.  The  maximum  number 
of  points  you  can  print  is  101. 
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** CHECKS 


**NEW  CASE^ 


Prints  the  current  status  of  the 
simulation  (what  commands  have  been 
executed  up  to  the  present  time)  and 
all  current  input  parameters.  The  user 
has  the  option  of  printing  out  the 
rather  lengthy  weight  and  volume 
distributions . 

Initializes  the  system  and  prepares  it 
for  a  new  simulation.  All  previous 
data  except  stored  weight  and  volume 
distributions  will  be  lost  to  the  user. 
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VIII.  SAMPLE  SIMULATIONS 


A.  Illustrative  example 

In  order  to  familiarize  the  user  with  MODSEC  system, 
a  short,  simple  example  will  be  set  up  and  executed. 

The  first  step  in  using  the  simulation  system  is  to 
choose  a  length/diameter  ratio,  and  overall  length  of  the 
submarine: 

L/D  =10,0 
Length  =  300.0  ft. 

The  position  of  the  separation  plane  is  selected  next. 

This  is  one  of  the  variables  that  the  user  will  probably 
wish  to  vary  during  the  simulation  to  see  its  effect  on  the 
capsule  motion.  Initally  let 

0-4- 

the  weight  distribution  is  the  next  and  most  difficult 
parameter  to  select.  The  more  accurate  and  detailed 
the  weight  distribution  is,  the  more  meaningful  the  simu¬ 
lation  will  be.  Ideally,  the  user  should  perform  a  de¬ 
tailed  weight  analysis  on  the  proposed  concept  and  generate 
a  list  of  components,  their  weights,  and  locations  in  the 
submarine.  This  list  is  then  used  as  input  data  to  the 
simulation  system.  For  the  first  example,  we  will  assume 
that  the  weight  distribution  is  as  given  in  the  following 
table.  This  distribution  jls  approximately  parabolic  over 


the  length  and  yields 
long  tons. 

a  total  submarine 

weight  of  4,046 

Weight 

Long ' 1  Cg 

Vertical 

56 

.05 

0 

10  8 

.10 

.15 

155 

.15 

-.30 

19  2 

.20 

.20 

275 

.25 

-.05 

254 

.30 

0 

271 

.35 

.5 

288 

.40 

-  .45 

298 

.  45  „ 

-.20 

301 

.50 

-.45 

29-7 

.55 

-.9 

289 

.60 

-.8 

272 

.65 

.2 

255 

.70 

.35 

225 

.75 

-.45 

19  2 

.  80 

-.20 

154 

.  85 

.10 

10  8 

.90 

0 

56 

.95 

0 
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The  vertical  centers  are  selected  at  random  and  vary  from 
-1.0  to  1.0  depending  on  what  portion  of  the  local  radius 
each  weight  is  from  the  centerline. 

The  thickness  of  the  capsule  shell  is  assumed  to  be 
2.0  inches.  The  thickness  of  the  main  hull  is  assumed  to 
be  1.7  inches.  Both  shells  are  assumed  to  be  constructed 
of  steel  with  a  density  of  487  lbs/ft3. 

The  size  and  location  of  the  spaces  which  are  likely 
to  be  flooded  in  the  event  of  a  casualty  are  computed 
next.  In  this  simple  example,  we  assume  that  there  are 
only  two  floodable  spaces: 

Volume  Long'l  Center  Vertical  Center 

8000  ft3  .70  .1 

2000  ft3  .20  -.3 

We  now  have  enough  data  to  begin  the  simulation  of 
the  modular  submarine  shown  in  Figure  VIII-1. 


Figure  yi£I-l.  Modular  submarine.  Example  1. 
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We  begin  by  dialing  the  telephone  number  of  the 
computer.  The  user  dialog  then  proceeds  as  follows: 
(NOTE:  User  responses  are  underlined.  The  RETURN  key 
must  be  depressed  after  each  line  of  user  response.  In 
the  documentation,  the  carraige  return  key  is  denoted  by 


ON-LINE  SYSTEMS  16 
1324  9-20-70 

♦BENCHMARK  SERVICE 
USER  NUMBER?. .99,88) 


INFORM  UPDATED  AUG  4 


SYSTEM?.  . 


•RUN  DSK:H0DSEC*SAV 


The  computer  prints  a  heading. 


The  user  supplies  his  assigned 
number.  The  comma  is  necessary. 

A  mask  is  printed  by  the  system. 
The  user  enters  his  assigned 
password  over  the  mask. 

If  the  user  responds  INFORM  to 
SYSTEM?..,  information  will  be 
printed  at  his  terminal. 

The  time-sharing  system  is  ready 
to  accept  commands  when  a  period 
(.)  is  printed.  If  an  error 
occurs  during  the  simulation,  a 
period  will  be  printed.  Type  STi 
to  continue  the  simulation.  The 
command  RUN  D3K:MODSEC^  loads  the 
MODSEC  System  from  a  disc  file 
and  begins  execution. 


IF  YOU  REQUIRK  INSTRUCTIONS.  TYPE  "HELP". 
**HELP 


**  ABBREVIATED  USER  INSTRUCTIONS: 


THE  PROGRAM  YOU  ARE  RUNNING  WILL  HELP  YOU  DESIGN  AN  ELLIPSOIDAL  SUB¬ 
MARINE  WITH  A  SEPARABLE  CAPSULE*  EACH  COMHAND  BELOW  EITHER  PERFORMS  A 
COMPLETE  FUNCTION,  PROMPTS  ADDITIONAL  QUESTIONS  FROM  THE  PROGRAM, 

OR  PROVIDES  YOU  WITH  INFORMATION*  THE  PROGRAM  SIGNALS  THAT  IT  IS 
READ!  FOR  A  COMMAND  BY  PRINTING  **  AT  THE  LEFT  MARGIN  OR  FOR  A  RE¬ 
SPONSE  BY  PRINTING  —  * 


« 
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COMMAND 

HELP 

INPUT  GEOMETRY 

INPUT  WEIGHTS 

INPUT  VOLUMES 

INPUT  DEPTH 

INPUT  SPEED 

INPUT  TRIM 

INPUT  SEAS 

UPDATE 

WEIGHTS 

DAMAGE 

TRIM 

SEPARATE 

ASCEND 

HULL 


DESCRIPTION 

PRINTS  USER  INSTRUCTIONS 

LETS  USER  SPECIFY  GENERAL  GEOMETRIC  CHARAC¬ 
TERISTICS  OF  VEHICLE 

LETS  USER  SPECIFY  WEIGHT  DISTRIBUTION  OF  VEHICLE 

LETS  USER  SPECIFY  FLOODABLE  VOLUME  DISTRIBUTION 

LETS  USER  SPECIFY  INITIAL  DEPTH 

LETS  USER  SPECIFY  INITIAL  VERTICAL  AND  HORIZONTAL 
VELOCITIES 

LETS  USER  SPECIFY  INITIAL  TRIM  ANGLE 

LETS  USER  SPECIFY  SURFACE  SEA  STATE 

ALLOWS  CHANGES  TO  BE  MADE  TO  WEIGHT  OR  VOLUME 
DISTRIBUTIONS 

CALCULATES  WEIGHTS  AND  CENTERS 

SIMULATES  SELECTIVE  FLOODING  OF  VOLUMES 

CALCULATES  SOU Y ANT  FORCE  AND  TRIM  OF  INTACT 
VEHICLE 

SIMULATES  CAPSULE  SEPARATION  FROM  HULL 

CALCULATES  CAPSULE  TRAJECTORY 

CALCULATES  REMAINING  HULL  TRAJECTORY 
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GRAPH  PLOTS  BEHAVIOR  OF  SEPARATED  CAPSULE  OR  REMAINING 

VEHICLE  (WHICHEVER  HAD  THE  HOKE  RECENT  TRAJECTORY 
COMPUTED) 


MOTION  DETERMINES  THE  ATTITUDE  OF  CAPSULE  ON 

THE  SURFACE  AND  PLOTS  THE  SEA  SURFACE  SPECTRUM 

PRINT  PRINTS  TRAJECTORY  DATA  IN  TABULAR  FORM 

CHECK  PRINTS  CURRENT  STATUS  OF  SIMULATION 

NEW  CASE  INITIALIZES  PROGRAM  FOR  A  NEW  SIMULATION 


FOR  ADDITIONAL  DETAILS.  CONSULT  PROGRAM  DOCUMENTATION 


-♦INPUT  GEOMETRY 

SPECIFY  L/D,  LENGTH.  SEPARATION  PUNE  —  10. 300.. 4 

—  INPUT  WEIGHTS 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IN).  AVERAGE 
HULL  SHELL  THICKNESS  (IN)  —  2*0,1 *7 

INPUT  SPECIFIC  WEIGHTS  OF  CAPSULE  MATERIAL.  AND 
OF  HULL  MATERIAL  (LBS/FT— 3)  —  487*0.487.0 


DO  YOU  WANT  TO  USE  OLD  WEIGHTS? 
YES  OR  NO  —  NO 


---  old  weights  might  have  been  stored 
on  a  disc  f i le 
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HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?  —  19 

INPUT  WEIGHTS  IN  LONG  TONS.  LONG“L  CG.  VERT  CG 
IN  GROUPS  OF  THREE 

56 . • 05 . o  ~"56  ton  weight t  ,  05L  fxoit  the  bow/  on  the  submaxine 

centexline 

108. *10. *15 


155, 

•  1 5  v  *  ♦  30 

192. 

•20,' 

•20 

275, 

•25,- 

-•05 

254. 

•  30»( 

3 

271. 

•  35.' 

•5 

288, 

•40, ' 

•  45 

298, 

•45, 

-•20 

301. 

•50, -•45 

297, 

•  55 « "  •  9 

289. 

•  60, 

•  8 

272. 

•  65, 

•  2 

255. 

•  70.' 

•  35 

225. 

•  75. 

•  45 

192. 

•80. 

-•20 

154. 

•85. 

•10 

108. 

•  9.0 

56.. 

95.  •( 

3 

store  weights?  (YES  OR  NO)  —  NO  — weights  Jtay  be  stored  pexmanently  on 

a  disc  file 
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♦♦UPDATE 


— used  to  change  weight  or  volume  array 

DO  YOU  WISH  TO  UPDATE  WEICHTS  OR  VOLUMES? 

TYPE  V  OP.  W  —  W 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED.  DELETED.  OR  ADDED? 

TYPE  C.D.  OR  A  — 

INPUT  PRESENT  WEIGHT.  LCG.  AND  VCC  —  192. -20, .20  - - weight  to  be  changed 

INPUT  NEW  WGT  OR  VOLUME.  AND  LCG.  VCG  —  192.  .20.->20 - and  replaced  by  this  new 

one 

DISTRIBUTION  UPDATED.  ANY  MORE  CORRECTIONS? 

TYPE  Y  OR  N  —  Y 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED.  DELETED.  OR  ADDED? 

TYPE  C.D.  OR  A  --  C_ 

INPUT  PRESENT  WEIGHT.  LCG.  AND  VCG  —  271.-35. -5 

INPUT  NEW  WGT  OR  VOLUME,  AND  LCG.  VCG  —  271. .35. -.6 

DISTRIBUTION  UPDATED.  ANY  MORE  CORRECTIONS? 

TYPE  Y  OR  N  —  _Y 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED.  DELETED.  OR  ADDED? 

TYPE  C.D.  OR  A 

INPUT  PRESENT  WEIGHT.  LCG.  AND  VCG  —  2 88.-4.. 45 

INPUT  NEW  WGT  OR  VOLUME,  AND  LCG.  VCG  —  288..4.-»45 

DISTRIBUTION  UPDATED.  ANY  MORE  CORRECTIONS? 

TYPE  Y  OR  N  —  _Y_ 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED,  DELETED.  OR  ADDED? 

TYPE  C.D.  OR  A  —  C 

INPUT  PRESENT  WEIGHT.  LCG.  AND  VCG  —  289.. 6. -8 

INPUT  NEW  WGT  OR  VOLUME.  AND  LCG,  VCG  —  289,'.6.~— 6 

DISTRIBUTION  UPDATED,  ANY  MORE  CORRECTIONS? 

TYPE  Y  OR  N  —  Y 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED,  DELETED.  OR  ADDED? 

TYPE  C.D.  OR  A  —  _£ 

INPUT  PRESENT  WEIGHT.  LCG.  AND  VCG  —  225, .75.-45 

INPUT  NEW  WGT  OR  VOLUME.  AND  LCG.  VCG  —  225. .75. -.45 

DISTRIBUTION  UPDATED.  ANY  MORE  CORRECTIONS? 

TYPE  Y  OR  N  —  Y 


changing  vertical 
centers  of 
selected  weights 
so  the  CG  will  be 
below  the  center- 
line  . 
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IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANCED.  DELETED.  OR  ADDED? 

TYPE  C.D.  OR  A  — 

INPUT  PRESENT  WEIGHT.  LCG.  AND  VCC  —  192. *8. -.20 
INPUT  NEW  WGT  OR  VOLUME.  AND  LCG.  VCG  —  195, .8.-»9 

\ 

DISTRIBUTION  UPDATED.  ANY  MORE  CORRECTIONS?  ' 

TYPE  Y  OR  N  —  N  | 

I 

DO  YOU  WISH  TO  STORE  THIS  DISTRIBUTION  FOR  FUTURE  USE  • 


TYPE  YOR  H  -M  _ we 

ENTER  SHI?  "IDENTIFIER”.  (5  CHARACTERS)  — 
DISTRIBUTION  HAS  BEEN  STORED  FOR  FUTURE  USE 


now  store  the  weights 
SUBOl - under  this  name 


♦♦ WEIGHTS  - calculates  total  weights  and  centers. 


★♦INPUT  VOLUMSaSaES  A  means  that  a 

—  character  just 

DO  YOU  WANT  TO  USE  OLD  FLuODABLE  VOLUMES? 
YES  OR  NO  —  NO 


"rubout"  was  typed  to  delete  the 
typed.  The  computer  then  typed  the 
deleted  character  and  another  h 


HOW  MANY  FLOODABLE  VOLUMES  — _2 

INPUT  FLOODABLE  VOLUMES.  LONG"L*  AND  VERTICAL  CENTROIDS 
IN  CROUPS  OF  THREE 

8000. »70, »i - 8000  cu.  ft.  of  space  is  located  . 7l  from  the  bow  and  . 1R 

2000. *20, -*3  above  the  centerline. 

STORE  VOLUMES?  (YES  Or  NO)  —  yes  - floodable  spaces  are  stored  on  a  disc 

file 

ENTER  SHIP  "IDENTIFIER".  (5  CHARACTERS)  --  SUBOl - under  this  name 


DISTRIBUTION  HAS  BEEN  STORED  FOR  FUTURE  USE  ’ 


★★TR I x  -“calculates  the  attitude  of  intact  vehicle.  Used  to  balance 
— — L  submarine 


INTACT  VEHICLE  IS  HEAVY  OVERALL  BY  1764619.  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  23*45  DEGREES 

♦♦INPUT  GEOMETRY 

SPECIFY  L/D.  LENGTH.  SEPARATION  PLANE  —  8. 8.300.. 4 
♦♦TRIM 

INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  671201.  POUNDS 
INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  23*45  DEGREES 
♦♦INPUT  GEOMETRY 

SPECIFY  L/D.  LENGTH.  SEPARATION  PUNE  —  a.a5.300..4 


changing  the 
L/D  ratio  to  get 
a  neutrally 
buoyant  submarine 
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i 

i 

i 

\ 


I 


1 1 

n 


♦♦TRIM 


INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  739556*  POUNDS 
INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  23*45  DEGREES 
♦♦INPUT  GEOMETRY 

SPECIFY  L/D.  LENGTH »  SEPARATION  PLANE  —  a»85,3Q0>*4 

♦♦..'EIGr.TS 


♦♦TRIM 


INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  495358*  POUNDS 
INTACT  VEHICLE  IS  TRIMMED  BON  DOWN  22*4?  DEGREES 
♦♦UPDATE 


adding  a 
weight  to  the 
stern  to  bring 
the  submarine 
to  an  even  keel 
attitude 


DO  YO'J  WISH  TO  STORE  THIS  DISTRIBUTION  FOR  FUTURE  USE 
TYPE  Y  OR  N  —  Y 

ENTER  SHIP  "IDENTIFIER",  (5  CHARACTERS)  —  SUB01 

IS  DISTRIBUTION  NEW  OR  AN  UPDATE  ? 

TYPE  N  OR  U 

DISTRIBUTION  HAS  BEEN  STORED  FOR  FUTURE  USE 
♦♦'..'EIGHTS 


DO  YOU  WISH  TO  UPDATE  WEIGHTS  OR  VOLUMES? 

TYPE  V  OR  W  --  _W_ 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED,  DELETED.  OR  ADDED? 
TYPE  C»D«  OR  A  —  _C 

INPUT  PRESENT  WEIGHT,  LCG.  AND  VCG  —  56, *95,0 

INPUT  NEW  WGT  OS  VOLUME,  AND  LCG.  VCG  --  100*  *95.0 

DISTRIBUTION  UPDATED,  ANY  MORE  CORRECTIONS? 

TYPE  Y  OR  K  —  N 


♦♦TP.  I M 


INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  397293*  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  8*58  DEGREES 

♦♦UPDATE 

DO  YOU  WISH  TO  UPDATE  WEIGHTS  Oa  VOLUMES? 

TYPE  V  OR  W  — _W 

IS  DISTRIBUTION  TO  BE  UPDATED  TO  BE  CHANGED.  DELETED,  OR  ADDED? 
TYPE  C.D.  OH  A  — •  C  V 
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INPUT  PRESENT  WEIGHT »  LCG,  AND  VCG  —  100, -9? 


INPUT  NEW  V.'GT  OR  VOLUME.  AND  LCG.  VCG  ~  120.  *95,0 


DISTRIBUTION  UPDATED,  ANY  MORE  CORRECTIONS? 
TYPE  Y  OR  M  —  N 


DO  YOU  WISH  TO  STORE'  THIS  DISTRIBUTION  FOR  FUTURE  USE 
TYPE  Y  OR  N  —  Y_ 

ENTER  SHIP  "IDENTIFIER" ,  (5  CHARACTERS)  —  SU301  -modific-  weight  distribution 

has  been  stored 

IS  DISTRIBUTION  NEW  OR  AN  UPDATE  ? 

TYPE  N  OR  U  --  jj 

IS  DISTRIBUTION  NEW  OR  AN  UPDATE  ? 

TYPE  N  OR  U  --  U_ 

DISTRIBUTION  HAS  BEEN  STORED  FOR  FUTURE  USE 


♦♦WEIGHTS 


**TRIM 


INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  352493-  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  1-73  DEGREES 

**IN?UT  GEOMETRY 

SPECIFY  L/D,  LENGTH,  SEPARATION  PLANE  —  8-9,300,-4 

**NSIGHTS 

**TRIM 

INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  235069-  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  1-68  DEGREES 

*WIS PUT  GEOMETRY 


L/d  changed 
some  more 
to  make  boat 
neutrally 
buoyant 


SPECIFY  L/D,  LENGTH,  SEPARATION  PLANE  —  8-95,300, -4 


**V;  EIGHTS 


INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  119558-  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  1-58  DEGREES 


Si'ilCI-TY  L/D*  LeJL GTH »  SLPALATIOM  PLA.»E  S*D*oOO»*4 


Best  Available  Cot 


**TRIM 


INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  235069*  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  1*68  DEGREES 

**I  N'PU  GEOMETRY 

SPECIFY  L/D.  LENGTH.  SEPARATION  PLANS  —  9.0,300. *4 


** WEIGHTS 
**TR I M 


INTACT  VEHICLE  IS  LIGHT  OVERALL  2Y  6221*  POUNDS 

INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  1*49  DEGREES 
**IN?UT  GEOMETRY 

final  L/~,  Length, 
Separation  plane 


SPECIFY  L/D.  LENGTH.  SEPARATION  PLANE  --  9*01,300, *4 — 
^.'EIGHTS 


**TRIM 


INTACT  VEHICLE  IS  HEAVY  OVERALL  BY 
INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN 
** INPUT  GEOMETRY 


16233*  POUNDS 
1*4?  DEGREES 


SPECIFY  L/D.  LENGTH,  SEPARATION  PLANE  —  9*0025,300,  *4 
WEIGHTS 


**TRIM 


INTACT  VEHICLE  IS  LIGHT  OVERALL  3Y 
INTACT  VEHICLE  IS  TRIMMED  30 W  DOWN 
DA MAGS 


600*  POUNDS 
1*48  DEGREES 


DO  YOU  NEED  A  LIST  OF  FLOODABLE  VOLUMES?  { YES  OR  NO)  —  YES 


SPACE  NO 

VOL  ( CU  *  FT  * ) 

LONG  CC- 

VERT  CG 

1 

8000*0 

0*700 

0*100 

2 

2000*0 

0*200 

-0*300 

HOW  MANY 

SPACES  DO  YOU  WANT 

TO  FLOOD?  — 

1 

INPUT  THE  NUMBERS  OF  FLOODED  SPACES 


— casualty  sustained  in 
one  space 


- space  no.  1  near  the  stern  has 

been  flooded 


Best  Available  Copy 


I 


; 

■i 

'  i 


**7RIM 


DAMAGED  INTACT  VEHICLE  IS  HEAVY  OVERALL  BY  511400* 
DAMAGED  INTACT  VEHICLE  WILL  EVENTUALLY  TRIM  BOW  UP 
INTACT  VEHICLE  WILL  TRIM  HYDROSTATICALLY  FROM  -1*43 
DEGREES  BEFORE  DAMAGE  TO  31*51  DEGREES  AFTER  DAMAGE 


POUNDS — intact  vehicle  will 
31*51  DEGREES  sink 


** INPUT  TRIM 


'i 

,  t 


'  J 

V 


SPECIFY  INITIAL  TRIM  ANGLE  (DEGREES)  —  31*51 - same  as  hydrostatic  trim  angle 

i  (it  need  not  be) 

** INPUT  SPEED 

SPECIFY  INITIAL  VERTICAL  AND  HORIZONTAL  SPEED  (FT/SEC!  —  -5*0, 0*0 — sinking  slowly 
**SSPARATE  capsule  has  separated  from  submarine 
**ASCSND  — trajectory  is  calculated 


we  are  interested  in  the  'j 

first  85  seconds  after 
**GP.A?H  —  we  would  like  to  separation 

see  the  results  in  graphical  form 

CHECK  PLOT  SWITCH  — Tvpagraph  plot  switch  moved  to  plot  mode, 

carraige  return  typed. 

I 


( 


‘INPUT  MAXIMUM  SIMULATION  TIME 


SECONDS  —  S5 


Best  Available  Cop\ 
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OM«i£H<:aow  in  t:: o  s  wm  n< «  •>; e-< « o m  o 


surma; :i:ne  escape  capsule  simulation 

CUR  ,  3-NOV-70 

CHARACTERISTICS  0?  CAPSULE  : 

LENGTH  (FT)  :  120*00 

LENGTH  TO  DIAMETER  RATIO  :  S-00 

SEPARATION  PUNE  LOCATION  :  0-40 

INITIAL  TRIM  ANGLE  (DEG)  :  31-51 

INITIAL  X  VELOCITY  (FT/SEC)  :  0-00 

INITIAL  Y  VELOCITY  (FT/SEC)  :  -5-00 

19*00j  X/L  :  + 

j  Y/L  :  * 

!  THETA  : 


16-00 


n-ooi 


7-00“ 


I 

I 

r» 


3-00H 


.  Y 


* 

j.* 


* 


jr 


?**•>*  J1 

'fit 

,*  - 


.* 


•*  .** 


,'jfr 


.***' 


- —  -  -  - i  ! t  i  i  i  i !  I  li  1  I )  I  li  I  I ■  "  i i-fr1 |  1  [ 


•^rr-K>. 


10 


20  30 

NON-DIMENSIONAL  TIME 


40 


300 


V 

E 

200 H 
T 

C 

L 

v 


1001 

m 

A 

U 

D 

E 

I 

'"N' 


G 

R 

E 

v 


-100 


S'  -200 

50 


5  2 


Best  Available  Copy 


4 


( 


- print  tarulsr  results  of  the  last  calculated  trajectory, 

in  this  case  the  capsule  trajectory 

KOVi  MANY  POINTS  CO  YOU  WISE  TO  PHIS?  <2  TO  14)  ~  _!£_ 


TIKE 

HOP.  1 2 
DISTANCE 

non- die 

OKTAL 

VELOCITY 

.ENSICNAL  RsSUihS 
VERTICAL 

DISTANCE  VELOCITY 

ATTITUDE 

DECREES 

A  1»  G  U  LA  R 
VELOCITY 

0-00 

0*00 

0*00 

0*00 

—G  •  Go 

61  -51 

0*00 

3*08 

-0-00 

-0*00 

-o*i  a 

-0*04 

-32*73 

-33*07 

6*16 

-0*00 

-0*00 

-0*21 

0*02 

-114*38 

-10*08 

9*24 

-0*00 

-0*00 

-0*03 

0*07 

-co  *  1 3 

24  •  62 

12*32 

-0*00 

-0*00 

0*23 

0*12 

-.7*22 

12*07 

13*40 

-0*01 

-0*00 

0*66 

0*16 

-29*13 

-17*53 

1S*43 

-0*02 

-0*00 

1*21 

0*20 

-65*35 

-11*93 

21*57 

-0*04 

-0*00 

1*93 

0*25 

-32*31 

13*44 

f 

24*65 

-0*05 

-0*01 

2*79 

0*30 

-32*72 

12*71 

27*73 

-0*10 

-0*02 

3*73 

0*31 

-25*45 

-8*34 

1 

30*81 

-0*18 

-0*03 

4-73 

0*34 

-65*74 

-12*73 

33*69 

-0*27 

-0*03 

5*84 

0*37 

-60*84 

4*39 

( 

35*3? 

-0*35 

-0*03 

7*03 

0*40 

-48*19 

12*96 

1  \ 

" 

40*03 

-0*47 

-0*05 

8*28 

0*41 

-23*48 

0*88 

j 

43*5? 

-0*68 

-0*07 

9*73 

0*42 

-48*27 

-11*70 

l 

DIMENSIONAL! ZED  ?. 

ES’JLTS 

TIKE 

HORIZONTAL 

VERTICAL 

ATTITUDE 

AMU u DAE 

DISTANCE 

VELOCITY 

DISTANCE 

VELOCITY 

DECREES 

VELOCITY 

1 

SECONDS 

FEET 

•T/SEC 

?  2S'V 

y T/SEC 

DEC/SEC 

0  •  00 

0*00 

0*00 

0*00 

-5*00 

31  *51 

0*00 

5*95 

-0-09 

-0*02 

-21*73 

-2*35 

-32*73 

-18*16 

11*30 

-0*18 

-0*01 

-25*77 

1*15 

-114*68 

-5*22 

17*35 

-0*26 

-0*01 

-9*09 

4*63 

-85*13 

12*69 

23*80 

-0*44 

-0*07 

27*67 

7*53 

-17*22 

6*25 

29*75 

-1*13 

-0*16 

73*42 

3*79 

-29*13 

-9*16 

35-70 

-2*65 

-0*29 

145*52 

12*65 

-65*95 

-6*20 

41  *6  5 

-4*32 

-0*27 

231*20 

16*03 

-62*31 

6*96 

47*o0 

-6*55 

-0*60 

334*35 

13*33 

-32-72 

6*53 

63  •  56 

-12*03 

-1  *25 

447*27 

19*57 

-26  *  46 

-4*32 

59*30 

-21-37 

-1  *77 

553*03 

21*17 

-65*74 

-6*62 

55*45 

-31-89 

-1*73 

700*21 

23*25 

“66  •  3-i 

2*27 

71  **10 

-42*43 

-1*95 

843*81 

24*87 

-48*19 

6*71 

( 

77-33 

-05* 99 

-2  •  99 

993*  73 

25*39 

-23* 4s 

0*45 

V 

34*15 

-61*29 

-4*07 

1167*70 

25*92 

-48*27 

-6*0  6 
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Since  the  capsule  tends  to  pitch  bow  down  when  separated 
and  move  backwards,  it  will  quite  likely  interfere  with 
the  remaining  vehicle,  k’ot  a  satisfactory  situation. 


RESET 
**KULL 
I  "PUT 


PLOT  SWITCH  TO  TEXT  MODE  AND  GIVE  A  CARRIAGE  RETURN 
— calculate  the  trajectory  of  the  remaining  hull 
MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  So 


CHECK  PLOT  SWITCH 
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0r  «t  **r" '  -  ft  I" 


n 


Z  b  1  .  C  mm  «  i  <  i  •%.*«  *  0 


A  CAi'u'lAtft  A&Twiv**  ••  f 


*;o*  kaky  points  do  you  wish  to  print  (2  to  mi  —  h 


NON— Di MENSI GNAL  RESULTS 

HORIZONTAL  VERTICAL  ATTITUDE  ANGULAR 

DISTANCE  VELOCITY  DISTANCE  VELOCITY  DEGREES  VELOCITY 


o*oo 

0*00 

0*00 

0*00 

-0*0? 

31  *31 

0*00 

2*52 

C  •  00 

0*00 

-0*39 

-0*25 

97*97 

-69*71 

o*Go 

“0*00 

— 0  •  oO 

-I  *26 

-0*44 

113*  38 

22*0? 

7*55 

-0*00 

-0*00 

-2*69 

-0*62 

73*60 

42  *'22 

I'J*  LO 

“0  •  01 

-0*00 

-4*3? 

-0*79 

67*28 

-16*03 

12*53 

~0*C3 

-0*01 

-6  •  54 

-0*94 

93*46 

-32*13 

15*03 

-0*06 

-0*01 

-9*07 

-1*07 

101*03 

8*57 

17*61 

-0*03 

-0*02 

-11*91 

-1*18 

84  *06 

27*49 

20*12 

-0*12 

-0*02 

-14*39 

-1*27 

73*13 

2*44 

22  *  O*; 

-0*1? 

-0*02 

-18*28 

-1*34 

82*24 

-21 *50 

25*15 

-0*22 

-0*02 

-21*73 

-1*40 

94*54 

-11*86 

27*67 

“0*28 

-0*03 

-25*31 

-1*44 

93*26 

11*89 

30  *1S 

-0*35 

“0*03 

-23*03 

-1*43 

82  *49 

16*72 

22*70 

-0-11 

-0  •  03 

-32*73 

-1*50 

77*51 

-1*07 

35*93 

-0*51 

-0*03 

-37*63 

-1  *53 

78*22 

1*92 

di::s::sic 

iNALIZED  2 

ESOLTS 

:ims 

HORIZON 

TAL 

VERT 

i  Cal 

ATT  ITU 02 

/i  *•  SUiiAr* 

DISTANCE  V 

ELOCITY 

DISTANCE 

VELOCITY 

DE02EES 

VELOCITY 

■coeds 

PT/SEC 

5  HOT 

I-T/3EC 

DEO/ 320 

G  •  vO 

0  *C3 

0*00 

0*00 

-5*C0 

31*51 

0*00 

5*95 

0*0-1 

0*05 

-70*61 

-18*57 

97*97 

-29*47 

11*30 

-G*30 

-0*10 

-226*43 

-33*44 

113*33 

9*33 

17*85 

-0*80 

-0*33 

-455 *43 

-46*85 

73*60 

17*55 

23*30 

-2*26 

-0*26 

-78o*9i 

-60* 14 

67*23 

-7*62 

23*75 

"6*35 

-0*42 

-1177*45 

-71*33 

93*45 

-13*53 

25*70 

~9  •  o4 

-0*87 

-1333*18 

-81 *45 

101*03 

3*62 

<.1  *65 

-15*00 

-1*25 

-2143*43 

-69  *69 

84*06 

1 •  62 

*17*60 

-21*32 

-1*32 

-2693*94 

-95*63 

73*16 

1*03 

Cw  •  C-'O 

-30*53 

-1*35 

-3290*81 

-102*13 

82*24 

-9*09 

S/J'wj 

-40 • 21 

“1  •  63 

—39 i 1 • o2 

-100*45 

1/4  •  64 

-5*01 

— oC • o? 

-1*35 

-4535*65 

-109*60 

9o*26 

6  •  03 

Vi  *40 

-32*23 

-2*16 

-5217*00 

-112*33 

82*43 

7*07 

77*35 

-  .  *  o5 

-2*12 

-5531 • £7 

-114*36 

77*51 

-0*46 

oJ  •  00 

-01*27 

-2*26 

-6774*05 

-116*18 

73-22 

0*31 
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The  re.~air.ir. 
staking  inter 


vehicle  trir.s  up  suddenly 
erer.ce  more  probable. 


and 


CiOveS 


forward. 


»*i;;?jt  seas 

I !;?'JT  SURFACE  SEA  STATE  — _6_  *c -  Capsuue  surfaces  in  a  state 

**:-‘.OiIO”S  — Calculates  capsule  response  to  a  state  6  sea. 

THE  CAPSULE  IS  FLOAT  I  KG  SOW  DOTH* 

SURFACE  SEA  STATE  IS  6 

SIGKIFIOAKT  LAV E  HEIGHT  =  12*542  FT 

SIGKIFICA'.T;’  HEAVE  AF.PL*  =  15*634  FT 

AVERAGE  PITCH  AM?L*  =  3*371  DEGREES 


sea 


CHECH  PLOT  SWITCH 


o  t*2  ^  O  l/J  t'J  t-’  ^ 


PLOT  SWITCH  TO  TEXT  MO 


DC  AMD  CIVi  A  CAK8IAG2  Hi 


+*2'.’3 


LXIiC-JTIOH  TIME: 
TOTAL  ELAPSED  TIMS: 


2  KIM*  51*62  SEC* 
6?  KIM*  4*32  SEC* 


Complete  simulation  took  l  hour  7  minutes 


HTUHM 


60 


4 


„  . .  «.  *•*/ 


.  ,  .-.-:^-»r  -‘T  7’  “1"  m 

JV  J\jy  *  *0  4  •  *•*-  ******* 

,>  *u  ■  ■■_•  bliili1 

INPUT'  AVERAGE  CAPSULE  SHELL  THlCRHtJi  tI..Si  i..-*..l.GL 

hull  shill  thickness  <IK»  —  QmO>  —  no  snell  weights 

IK?'J?  SPECIFIC  WEIGHTS  OF  CAPSULE  KATEP.IAL.  AND 
or  HULL  MATERIAL  (LliS/m*3)  —  Q».0‘ 

])0  YOU  WANT  '10  USE  OLD  WEIGHTS? 

YDS  Oft  NO  —  yes  —  mistake,  did  not  want  to  use  old  weights 

»»in?ut  weights 

I'UP'JT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IK I*  AVERAGE 

hull  shell  thick:: ess  uni  —  o- *o* 

IKPUT  SPECIFIC  WEIGHTS  0?  CAPSULE  MATERIAL*  AND 
OF  HULL  MATER I AO *L  (LBS/FT**3)  —  »Q*  . 

do  you  v:a:;t  to  use  old  weights? 

YES  OR  HO  —  no 

HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?  — _2_ 

IHPUT  WEIGHTS  IK  LONG  TOSS.  LONG" L  CG.  VEST  CG 
IS  GROUPS  0?  THREE 

200 •  1  •  2 » *  0 

200* ** 4*0  . 


STORE  WEIGHTS?  (YES  OR  SO)  ~  _N0 


SPECIFY  LKD*  LENGTH*  SEPARATION  PLAKx.  3*112**5 
SIGHTS 
**S5?A2A?S 

V-  -  uUSo  i.-sgo  AS  initial  TRIM  aSGLS — cannot  separate  without  initial 
OOl  WPUT*  trim  ‘  '  trir-  an^le 

SPECIFY  INITIAL  TRIM  ANGLE  (DEGREES)  — _0_ 

»*SE?ASATE 

WARNING— NO  CASULTY  SPECIFIED.  DO  YOU  WANT  TO  CONTINUE?  (YES  OR  NO)  -  YES 

»»A SCENE 

^ "*  V  J  Vy'*  I  \*y  1  Ort  T 1 1 lit  ii  ii’v  0.»  u  S  ** 


200 


CHECK  PLOT  SWITCH 


100  •OOf 


SO*OOt 


submarine  escape  capsule  simulation 

ONR  ,  3-N0V-70 

characteristics  of  capsule  : 

LENGTH  (FT)  :  S>6«00 

LENGTH  TO  DIAKKTKit  RATIO  :  a*00 

SEPARATION  PLANE  LOCATION  :  0*o0 

INITIAL  TRIM  ANGLE  (DEG)  :  0-00 

INITIAL  X  VELOCITY  CFT/SSC)  :  0-00 

INITIAL  Y  VELOCITY  (FT/SEC)  s  0-00 

X/L  :  + 

Y/L  :  * 

THETA  :  • 

Longitudinal  CG:  .3 
Vertical  CG:  0 

Center  of  Buoyancy:  .3125 
Hydrodynamic  Center:  .285 
Buoyant  Force:  10.176  Tons 


3  60*00j 


.  k*  *  ' 
/•  • 


I  40*00" 


*  */. 

«  £  •  . 
V  /  ” 


20*00? 


■K>  I  m*Wt*i  H  KiH-W  >  H+&)*A4r&r^\-?.W.<r+,nKy>  iWKUmi; 
40  80  120 

NON-DIMENSIONAL  TIME 


0 


RESET  PLOT  SWITCH  TO  TEXT  MODE  AMD  GIVE  A  CA  HR  I  AGE  RETURN 


♦  GEOMETRY 

SPECIFY  LAD*  LEH'GTH*  SEPARATION  PLANE  —  o»115t  -o 

**asce;:d 


INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  250 


**GiU?H 


CHECK  PLOT  SWITCH 
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RESET  PLOT  SWITCH  TO  TEXT  MODE  AND  GIVE  A  CARRIAGE  RETURN  —  J 
** INPUT  WEIGHT 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IN).  AVERAGE 
HULL  SHELL  THICKNESS  (IN)  —  Q..Q. 

INPUT  SPECIFIC  WEIGHTS  0?  CAPSULE  MATERIAL.  AND 
OP  HULL  MATERIAL  ( L3S/FT**3)  --  0-.0» 

DO  YOU  WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  HO  --  U0_ 

HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?  —  2_ 

INPUT  WEIGHTS  IN  LONG  TONS.  LONG"L  CG.  VERT  CG 
IN  GROUPS  OF  THREE 
190. -2.0 
190. *4.0 

STORE  WEIGHTS?  (YES  OR  NO)  —  NO 

** WEIGHTS 

**ASCEND 

INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  250 
**GRAPH 

CHECK  PLOT  SWITCH 
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v-3  in 


*10 

43*00+ 


suemarise  escape  capsule  simulation 

GEE  .  3-E0V-70 

CHARACTERISTICS  0?  CAPSULE  I 

LENGTH  (?T)  :  5? *50 

LENGTH  TO  DIAMETER  RATIO  :  5*00 

SEPARATION  PLANS  LOCATION  :  0*50 

INITIAL  TRIP.  ANGLE  (DEG)  :  0*00 

INITIAL  X  VELOCITY  I  FT/SEC)  s  0*00 

INITIAL  Y  VELOCITY  (rT/SEC)  :  0*00 

X/L  :  + 

Y/L  :  * 

TKETA  :  • 

Longitudinal  CG.:  0.3 
Vertical  CG:  0 

Center  of  Buoyancy:  0.3125 
Hydrodynamic  Ctr.:  0.2845 

Buoyant  Force:  47.522  Tons 


33*00-) 


l 

! 


I 


&SS2T  PLOT  SWITCH  TO  TEXT  MODE  AMD  GIVE  A  CARRIAGE  RETURN  —  l 
★★INPUT  Vi  EIGHT 


INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IK).  AVERAGE 
HULL  SHELL  THICKNESS  (IN)  —  Q«tO« 

INPUT  SPECIFIC  WEIGHTS  OF  CAPSULE  MATERIAL*  AND 
OF  HULL  MATERIAL  <LBS/FT**3)  --  Q. .0* 

DO  YOU  WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  NO  —  N0_ 

HOW  MANY  WEIGHTS  ARE  TO  HE  INPUT?  —  2 

INPUT  WEIGHTS  IN  LONG  TONS.  LONG“L  CG,  VERT  CG 
IN  GROUPS  OF  THREE 
105. *2.0 
195. *4.0 


STORE  WEIGHTS?  (YES  OR  NO)  —  KO 

★★weights 

★★ASCEND 

INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  250 


★★GRAPH 


CHECK  PLOT  SWITCH 


10 

45*00+ 


SUBMARINE  ESCAPE  CAPSULE  SIMULATION 

o:::;  .  s-i.ov-vo 

CHARACTERISTICS  0?  CAPSULE  : 

LENGTH  (FT)  :  57*50 

LENGTH  TO  DIAMETER  RATIO  :  5*00 

SEPARATION  PLANE  LOCATION  :  O'*  60 

INITIAL  TRIM  ANGLE  (DEG)  :  0*00 

INITIAL  X  VELOCITY  (FT/SEC)  :  0*00 

INITIAL  Y  VELOCITY  (FT/SEC)  :  0*00 

X/L  :  + 

Y/L  :  * 

THETA  : 


Longitudinal  CG: 
Vertical  CG: 

Center  of  3uoyancy: 
Hydrodynamic  Ctr: 
Buoyant  Force: 


0.3 

0 

0.3125 
C.2345 
37.52  Tons 


26*004 


^.•**** 


15*00+ 


*H'***+. 


‘+H+T-H-H+ 


^HH44+ 


NON-DIMENSIONAL  TIME 


To  investigate  the  capsule  behavior  for  a  longer  time  we 
need  only  to  increase  the  simulation  time 


S£SST  PLOT  SWITCH  TO  TEXT  MODE  AND  GIVE  A  CARRIAGE  RETURN  —  ). 

**asck:;c 


INPUT  MAXIMUM  SIMULATION  TIMS  IN  SECONDS  --  400 


**G?-A?n 


CHECK  PLOT  SV-ITCK 


R'ESET  PLOT  SWITCH  TO  TEXT  MOD'S  AIID  CIVS  A  CARRIAGE  RETURN* 
♦♦INPUT  WEIGHTS 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IN It  AVERAGE 
HULL  SHELL  THICKNESS  (IN)  —  Q«.Q. 

INPUT  SPECIFIC  WEIGHTS  OF  CAPSULE  MATERIAL.  AND 
OF  HULL  MATERIAL  (LBS/FT**3>  —  Q.  .0* 

DO  YOU  WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  NO  —  X0_ 

HOW  MANY  WEIGHTS  ARE  TO  32  INPUT?  —  2_ 

INPUT  WEIGHTS  IN  LONG  TONS.  LONG" L  CG.  VERT  CG 
IN  GROUPS  OF  THREE 
105. » 5.0 
195. •4S.0 

STORE  WEIGHTS?  (YES  OR  NO)  —  NO 

♦♦WEIGHTS 

♦♦ASCEND 

INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  250 
♦♦GRAPH 


CHECK  PLOT  SWITCH 
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m  o  rs  m  a,  •»;  r*:  c- 


*10 

oQ»OG| 

I 

I 


SuBMAiillii  ESCAPE  CAPS'JLE  SItULATIOi; 

o::h  t  o-:;ov-?o 

CHARACTERISTICS  0?  CAPSOLE  : 

LEHSTK  </";•>  :  5?.  50 

LEHOTH  TO  DIAMETER  RATIO  :  5«00 

Si?.w;iA .  10. <  PLin.r  LOCaiiu.  «  o*oO 

&  i  >  I •>  ( Cu  w  )  •  0  •  uu 

IIwTIi.i.  X  V iuOCIT Y  ifT/S^C)  :  0*00 

INITIAL  V  VELOCITY  (rT/Sl-X)  :  0*00 

X/L  :  + 

Y/L  :  * 

THETA  : 


4  Longitudinal  CG:  .4 

Vertical  CG:  0 

Center  of  Buoyancy  .3125 

Hydrodynamic  Cir. :  .2845 

*  Buoyant  Force:  37.52  Tons 

40 -OO-,- 


200 


V 

I 

160'" 

I 

C 

l 

A 

rp 

#7. 

120 1 
T 
ij 
D 
£ 

I 


c 


£ 

S 

40 


0 


KO:;-DIK£}JSIOrlAL  TIME 
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RESET  PLOT  STITCH  TO  TEXT  MODE  A  HD  GIVE  a  CARRIAGE  RETURN 
**It:?UT  HEIGHTS 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IK>.  AVERAGE 
HULL  SHELL  THICKNESS  .( IN)  —  O-.O- 

INPUT  SPECIFIC  HEIGHTS  OE  CAPSULE  MATERIAL.  AND 
0?  HULL  MATERIAL  <  LBS/*’T**3>  —  Q..0» 

DO  YOU  WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  NO  —  NO. 

hoy;  many  weights  are  to  be  INPUT?  —  2 

INPUT  WEIGHTS  IN  LONG  TONS.  10KG"L  CG,  VERT  CG 
IN  GROUPS  0?  THREE 
1-35.  '2bt0 
195. »45.0 

STORE  WEIGHTS?  <YES  OR  NO)  --  NO 
'EIGHTS 
**ASCEXD 

INPUT  MAXIMUM.  SIMULATION  TIME  IN  SECONDS  --  230 

*+ GRAPH 

CHECK  PLOT  SWITCH 
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SUBKAifINj  EsCAPf 

C«8iw  “L;.  *"»«■» 

•ssssffs&r' 

WlWAi  TflJM  ItT  OCAno?;  * 

I«TlAt  -  i rvrnlf  ulS  (D^J  , 
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ou°yant  PorJerUi‘; 


_v.  ***’ 
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■  *~WnU*ZLi, 


RESIT  PLOT  SWITCH  TO  TEXT  MODE  AMD  GIVE  A  CARRIAGE  RETURN  —  \ 


♦♦INPUT  Vi  EIGHTS 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IN).  AVERAGE 
HULL  SHELL  THICKNESS  (IK)  —  Q..O« 

INPUT  SPECIEIC  WEIGHTS  OF  CAPSULE  MATERIAL.  AND 
OF  HULL  MATERIAL  (L3S/FT**3)  --  Q..Q. 

GO  YOU  WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  NO  —  JiO_ 

HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?  —  2. 

INPUT  WEIGHTS  IN  LONG  TONS.  LONG’L  CG.  VERT  CG 
IN  GROUPS  OF  THREE 
1S3. »22»0 
ISo. *42.0 

STORE  WEIGHTS?  (YES  OR  NO)  --  NO 

♦♦WEIGHTS 

♦♦ASCEND 

INPUT  MAXIMUM  SIMULATION  TIMS  IN  SECONDS  --  250 
♦♦GRAPH 

CHECK  PLOT  SWITCH 
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NOH-DIKSSSIONAL  TIMS 


WWWWOMO  5^  •— i  c  C  H  ^  H  >  Mf*OM 


RESET  PLOT  SWITCH  TO  TEXT  MODE  AND  GIVE  A  CARRIAGE  RETURN  —  ) 
♦♦INPUT  WEIGHTS 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IN),  AVERAGE 
HULL  SHELL  THICKNESS  (IN)  —  Q.,Q. 

INPUT  SPECIFIC  WEIGHTS  OF  CAPSULE  MATERIAL,  AND 
OF  HULL  MATERIAL  (LBS/FT**3>  —  Q.,Q» 

DO  YOU  WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  NO  —  NO 

HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?  —  2_ 

INPUT  WEIGHTS  IN  LONG  TONS,  LONG"!  CG.  VERT  CG  .. 

IK  GROUPS  OF  THREE 
195, *19.0 
195, »42,0 


STORE  WEIGHTS?  (YES  OR  NO)  —  HO 
EIGHTS 
♦♦ASCEND 


INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  250 
♦♦GRAPH 


CHECK  PLOT  SWITCH 
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HESS'!  PLOT  S VI ITCH  TO 


TEXT  MODE  AND  GIVE  A  CARRIAGE  RETURN 


♦♦END 


2  MIN.  55-13  SEC- 
83  MIN-  13-?3  SEC- 


EXECUTION  TIME: 
TOTAL  ELAPSED  TIME: 


C.  Simple  Check  On  Capsule  ar.d  Hull  Behavior 


If  the  submarine  is  exactly  symmetrical,  the  capsule 
should  behave  just  like  the  remaining  vehicle.  That  is, 
if  the  capsule  rises,  moves  forward  and  pitches  up,  the 
remaining  vehicle  should  rise  move  backward,  and  pitch  down. 
In  other  words,  there  should  be  a  vertical  plane  of  sym¬ 
metry  for  the  motion. 

In  order  to  check  this,  the  same  geometry  and  weight 
distribution  used  in  The  Simple  Parametric  Analysis  will 
be  used.  Since  we  are  now  concerned  with  the  remaining 
vehicle  motions  as  well  as  the  capsule  trajectory,  two 
more  200-ton  weights  will  be  input  on  the  centerline  at 
0.6L  and  0.8L.  The  L/D  is  5,  the  length  is  115  feet  and 
the  separation  plane  is  located  at  0.5L. 

The  graphical  output  confirms  the  expected  symmetrical 
behavior. 

Note  that  the  command  abbreviations  are  used  extensively 
in  this  example. 


•run  ds:\:kcdssc.sav 
EXECUTION 

IF  YOU  REQUIRE  INSTRUCTIONS.  TYPE  "HELP". 

*»!N?UT  JEOXETHY 

SPECIFY  L/D.  LENGTH.  SEPARATION  PLANE  —  6U00..5 
** INPUT  WEIGHTS 

INPUT  AVERAGE  CAPSULE  SHELL  THICKNESS  (IN).  AVERAGE 
HULL  SHELL  THICKNESS  <IK)  —  0.0 

INPUT  SPECIFIC  WEIGHTS  OF  CAPSULE  MATERIAL.  AND 
OF  HULL  MATERIAL  ( LBS/FT**3)  —  0.0 

DO  YOU  ’WANT  TO  USE  OLD  WEIGHTS? 

YES  OR  NO  —  NO 

HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?  —  4 

INPUT  WEIGHTS  IN  LONG  TONS.  LONG“L  CG.  VERT  CG 
IN  GROUPS  0?  THREE 
200. >2.0 
200. »4.0 
200.-6.0 
200. -3,0 


STORE  WEIGHTS?  (YES  OR  NO)  --  NO 

»*trim 

TRIM  ANGLE  CANNOT  BE  FOUND  UNTIL  WEIGHTS  AND  CENTERS  ARE  FOUND 
♦^WEIGHTS 


**TRIM 


INTACT  VEHICLE  IS  HEAVY  OVERALL  BY  451588.  POUNDS 
CG  AND  CB  COINCIDENT.  TRIM  ANGLE  INDETERMINATE 

**I  C 

SPECIFY  L/D.  LENGTH,  SEPARATION  PLANE  --  5. 108, .5 
**T 

INTACT  VEHICLE  IS  HEAVY  OVERALL  BY  10346?.  POUNDS 
CC  AND  CB  COINCIDENT.  TRIM  ANGLE  INDETERMINATE 

*»I  G 


SPECIFY  L/D.  LENGTH,  SEPARATION  PLANS  —  5.110, «5 
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♦♦TRIM 


I  STAC?  VEHICLE  IS  HEAVY  OVERALL  BY  7912*  POUNDS 

CG  AND  CB  COINCIDENT,  TRIM  ANGLE  INDETERMINATE 

♦♦I  G 

SPECIFY  L/D,  LENGTH.  SEPARATION  PUNE  —  5.115, .5 
♦♦TRIM 

INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  2465S9*  POUNDS 
CG  AND  CB  COINCIDENT,  TRIM  ANGLE  INDETERMINATE 

♦♦SEPARATE 


YOU  MUST  INPUT  AN  INITIAL  TRIM  ANGLE 
♦♦INPUT  TRIM 

SPECIFY  INITIAL  TRIM  ANGLE  (DEGREES)  —  0 


♦♦SEPARATE 

WARNING— NO  CASULTY  SPECIFIED*  DO  YOU  WANT  TO  CONTINUE?  (YES  OR  NO)  —  YES 
♦♦ASCEND 

INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  200 
♦♦GRAPH 


CHECK  PLOT  SWITCH 


>■-.  <o  t-  •<  o  w  w  o  z  uhj  n,  (K  •<  h  m  o  «  &.  o  as  c-« 


*10 

20*00f 


SUBMARINE  escape  capsule  simulation 
qxr  .  3-kov-?o 
CHARACTERISTICS  OF  CAPSULE  t 

LENGTH  <FT>  :  57*50 

LENGTH  TO  DIAMETER  RATIO  !  5*00 

SEPARATIOtl  PLANE  LOCATION  s  0*50 

INITIAL  TRIM  /.HOLE  (DEO)  :  0*00 

INITIAL  X  VELOCITY  (FT/SEC)  :  0*00 

INITIAL  Y  VELOCITY  (FT/SECI  t  0*00 

X/L  s  ♦ 

Y/L  :  * 

THETA  i  • 


16*00* 

?  r. 


12*00f 


\  / 

\/ 


0 


NON-DIMENSIONAL  TIME 
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RESET  PLOT  SWITCH  TO  TEXT  MODE  AND  GIVE  A  CARRIAGE  RETURN 


♦♦HULL 

INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  200 
♦♦GRAPH 


CHECK  PLOT  SWITCH 


O  f«} 


RESET  PLOT  SWITCH  TO  TEXT  MODE  AMD  GIVE  A  CARRIAGE 
RETURN 

♦♦INPUT  SPEED 

SPECIFY  INITIAL  VERTICAL  AND  HORIZONTAL  SPEED  (FT/SEC)  —  0.10 
♦♦ASCEND 

INPUT  MAXIMUM  SIMULATION  TIME  IN  SECONDS  —  200 
♦♦GRAPH 

CHECK  PLOT  SWITCH 
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RESET  PLOT  SWITCH  TO  TEXT  MODE  AMD  GIVE  A  CARRIAGE  RETURN  — 
♦♦CHECK 


PRESENT  SIMULATION  STATUS 


GEOMETRICAL  PARAMETERS  HAVE  BEEN  CALCULATED 


WEIGHTS  AND  CENTERS  HAVE  BEEN  CALCULATED 

INTACT  VEHICLE  ATTITUDE  HAS  BEEN  CALCULATED 

CAPSULE  ACCENT  HAS  BEEN  CALCULATED 

LENGTH /DIAMETER  RATIO  » 

5*00 

TOTAL  VEHICLE  LENGTH (FT) 

«  116*00 

SEPARATION  PUNE  LOCATION 

= 

0-50 

CAPSULE  SHELL  THICKNESS!  INO 

= 

0*00 

HULL  SHELL  THICKNESS ( IN- ) 

= 

0*00 

INITIAL  HORIZONTAL  VELOCITY 

= 

10*00 

INITIAL  VERTICAL  VELOCITY 

= 

0*00 

INITIAL  DEPTH (FT  1 

0*00 

INITIAL  TRIM  ANGLE (DEGREES) 

= 

0*00 

LOG  INTACT  VEHICLE 

= 

0*600 

VCC  INTACT  VEHICLE 

= 

0*000 

LCG  CAPSULE 

C 

0*300 

VCG  CAPSULE 

s 

0*000 

DO  YOU  WISH  TO  INSPECT 
TYPE  Y  OR  N  -  Y 

THE  CURRENT 

WEIGHT  DISTRIBU 

CURRENT  WEIGHT  DISTRIBUTION  IS: 

200*00 

0*20 

0*00 

200*00 

0*40 

0*00 

200*00 

0*60 

0*00 

200*00 

0*60 

0*00 

I 

I 

I 

I 

I 


DO  YOU  WISH  TO  INSPECT  THE  CURRENT  VOLUME  DISTRIBUTION? 
TYPE  Y  OR  t;  -  KO 

**I KPUT  SEAS 

INPUT  SURFACE  SEA  STATE  —  5 
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Z2L2Z n  \rr-7i 


i 

3K 


RES ST  PLOT  SWITCH -TO  TEXT  MODE  AND  GIVE  A  CARRIAGE  RETUSK  - 


**END 


EXECUTION  TIME; 
TOTAL  ELAPSED  TIME; 


2  MIN.  51*05  SEC 
55  MIN.  ?«30  SEC 


FREQUENCY*  SAD/SEC 


IX.  COMPUTER  PROGRAM  DOCUMENTATION 


IX- 1.  Main  Program 

The  main  program  will  be  referred  to  hereafter  as 
the  monitor.  It  provides  an  interface  for  all  of  the 
subroutines,  accepts  commands,  accepts  data  from  the 
terminal,  prints  error  messages,  and  maintains  a  status 
array  to  assist  the  user  in  performing  a  simulation  in 
a  logical  fashion. 

The  monitor  primarily  uses  the  reserved  COMMON  areas 
as  the  medium  for  transferring  input/output  data  between 
itself  and  a  subroutine.  This  allows  many  of  the  sub¬ 
routines  to  be  called  without  using  any  subroutine 
arguments. 

The  monitor  notifies  the  user  that  it  is  ready  for 
a  command  by  printing  a  double  asterisk  (**) .  Similarly 
the  user  is  expected  to  respond  with  alpha  or  numeric 
data  (within  the  context  of  the  question)  when  a  double 
dash  ( — )  is  encountered.  Each  monitor  command  is  scan¬ 
ned  for  accuracy,  and  if  a  command  is  net  recognised,  a 
message  so  indicating  is  printed,  and  a  double  asterisk 
is  printed  indicating  that  the  user  is  expected  to  input 
the  proper  command. 

Since  several  subroutines  require,  as  input,  data 
that  are  computed  by  other  subroutines,  a  sequential  ex¬ 
ecution  hierarchy  must  be  followed  by  the  user.  An  array 
has  been  set  aside,  ITEST  (30) ,  to  contain  the  status  of 
those  subroutines  which  must  be  executed  sequentially. 

Before  a  subroutine  further  in  the  logical  sequence  can 
be  run,  a  test  is  made  of  the  appropriate  pre-reqaisite 
subroutines.  If  all  conditions  are  met,  execution  continues; 
if  not,  an  error  message  explaining  the  deficiency  is 
printed  at  the  terminal,  and  the  user  is  returned  to  the 
double  asterisk. 

Error  Messages 

1)  "WEIGHTS  AND  CENTERS  CANNOT  BE  FOUND  UNTIL 
GEOMETRY  IS  SPECIFIED"  -  The  user  has  attempted  to  ex¬ 
ecute  the  subroutine  which  computes  the  weights  and 
centers  before  specifying  the  geometry  of  the  vehicle. 
Corrective  action:  specify  vehicle  geometry  (i.e., 

INPUT  GEOMETRY) ) . 

2)  "CENTERS  CANNOT  BE  FOUND  UNTIL  THE  WEIGHTS  HAVE 
BEEN  INPUT"  -  The  user  tried  to  calculate  the  longitudinal 
ar.d  vertical  centers  of  gravity  (i.e.,  WEIGHTS))  but  he 
had  not  specified  a  weight  distribution.  Corrective 
action:  input  a  weight  distribution  through  keyboard  or 
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specify  a  previously  saved  distribution  through  the  INPUT 
WEIGHTS  command  and  then  give  the  command  WEIGHTS). 

3)  "YOU  HAVE  NOT  INPUT  A  LIST  OF  FLOODABLE  VOLUMES"  - 

The  user  requested  a  damage  simulation  (DAMAGE))  but  has 
not  yet  specified  the  floodable  volume  list.  Corrective 
action:  type  INPUT  VOLUME)  and  specify  the  floodable 

volumes. 

4)  "VEHICLE  HAS  ALREADY  SEPARATED"  -  User  requested 

the  hydrostatic  characteristics  of  the  intact  vehicle 
(TRIM) )  subsequent  to  commanding  SEPARATE) .  Corrective 
action :  run  NEW  CASE) . 

5)  “TRIM  ANGLE  CANNOT  BE  FOUND  UNTIL  WEIGHTS  AND 
CENTERS  ARE  FOUND"  -  User  attempted  to  find  trim  of 
intact  vehicle  (TRIM))  before  specifying  the  weight  array. 
Corrective  action:  use  the  INPUT  WEIGHTS)  command  and 
try  again. 

6)  "YOU  MUST  INPUT  AN  INITIAL  TRIM  ANGLE"  -  User 
gave  the  SEPARATE)  command  but  failed  to  specify  a  trim 
angle.  Corrective  action:  indicate  trim  angle  by  com¬ 
manding  INPUT  TRIM)  and  then  give  the  SEPARATE)  command. 

7)  "WARNING  —  NO  CASUALTY  SPECIFIED.  DO  YOU  WANT 

TO  CONTINUE?  (YES  OR  NO) "  -  This  is  not  an  error  mes¬ 

sage.  It  is  printed  to  remind  the  user  that  he  has  not 
yet  damaged  the  intact  vehicle.  He  has  the  option  of 
continuing  (answer  -  YES) )  or  of  specifying  a  casualty 
(answer  -  NO) ,  then  give  DAMAGE)  command) . 

8)  "CAPSULE  HAS  NOT  3EEN  SEPARATED  FROM  MAIN  HULL"  - 
The  command  ASCENT)  or  HULL)  has  been  given  before  the 
SEPARATE)  command.  Corrective  action:  give  the  SEPARATE) 
command  and  try  again. 

9)  "NO  DATA  TO  PLOT"  -  User  tried  to  plot  the  motion 
of  the  separated  capsule  or  hull  (GRAPH))  but  has  not  com¬ 
puted  the  capsule  or  hull  motions.  Corrective  action: 

use  the  ASCENT)  or  HULL)  command  and  try  GRAPH)  again. 

10)  "SORRY — CAPSULE  HAS  SUNK"  -  The  user  has  tried 
to  inspect  the  capsule's  behavior  at  the  surface  (MOTION)) 
but  the  capsule  has  been  found  not  buoyant.  Corrective 
action;  user  should  input  new  geometric  characteristics 
and/or  a  new  weight  distribution  (INPUT  WEIGHTS)),  (INPUT 
GEOMETRY) ) ,  and/or  specify  new  floodable  volumes  (INPUT 
VOLUMES)),  and  try  again. 


11)  "NOT  ENOUGH  INPUT  DATA  TO  CALCULATE  SURFACE  MOTIONS. 
INPUT  SEA  STATE"  -  User  attempted  to  investigate  capsule 
behavior  on  the  surface  (MOTIONS))  without  specifying  the 
sea  state.  Corrective  action:  input  the  desired  surface 
sea  state  (INPUT  SEAS)). 

12)  "ILLEGAL  INPUT  COMMAND,  TRY  AGAIN."  -  The  user's 
command  was  not  recognized  by  the  monitor  as  a  legal  com¬ 
mand.  Corrective  action:  try  again. 
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IX- 1. 1  Flow  Chart  for  Main  Program  (M0NITR) 
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IX- 1. 2  Listing  of  Main  Program  (M0NITR) 


( 
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+ 

* 

* 

* 


1 

1 

2 


3 

10 


20  OS 
15 
20 
25 


50 

1100 

1000 

1500 

9GD9 

660 

600 

6001 

6002 


6005 


COMMON/ I  NPUT/FLEN.SeP » i’LDHAI*  TRIM » VELX.  V  ELY » V  ELL 
COMMON/ I  HOUT/FLC  j»IX  »  FLCBC? » FLCBEV ,  CGCAP1 »  CGitVl .  VOLITV .  ElCMii  * 
SI  TV  i  YOLCAP.  KLCiiCP.  SCAP.HGCAP,  V OLit V »  SRV  *  HGitV  , BAREA.PROJCP. 
PKOJKV  .TOTWY » ii  v  ,-;t  .  CAPWT .  FLCOIK.  vcg  itv  .  flcghv  .VCGHV » FLCOCP. 
VCGCAP .TIME ( 101 ) » ANGL1 101 ) » VDI SI  101 )  .liDISt  101 )  , 

AVELUOi  )  .WELU01)  .HVELliOl) 

INTEGER  COM 
DIMENSION  IFLDS?  (10) 


DIMENSION  ITST (30) t COM 1 20) » ME SSI  15) , MKSS2I ?> • FLOOD <3. 


100)  tN; 

DATA  1 

’  S  *  ,  ’  A 


IG:iT  (3*100) 

ESS(I),I=1»16>,(MESS2(J)»J=1,7)/'E’» 

« '  H  *  » '  G  * » *  ?  *  *  ’  T  * » *  C  * .  '  M  * » *  * » 'N ' ,  'U ' 


•I1. 

,'G* 


V 


Vi\  *D*. 
•N*,*S* 


•D' 


CALL  EKESETIO) 
NUMr SP=0 
GO  TO  2008 
DO  10  1=1.30 
ITS'I'l  I  )=Q 


N  oM  /iG'i  =0 
NuMFLD=0 


CA?'i=0* 

HULL'i'=0* 

dinit-o* 
iLCGi::=o* 
VCGIT7=0* 
/ LCGC?=0* 


7CGCA?=0* 
i  Ls«l'i=0* 
5E?=0- 


FLD?.AT=0* 


TKiy.=0* 

TYPE  15 

FORMAT! '  IP  YOU  REQUIRE  INSTRUCTIONS.  TYPE  ** KELP”  •  * ) 
TYPE  25 

FORMAT  (lh  . '+*',S) 

ACCEPT  30. (COM! 1 1.1=1, 20) 
x  0 .'.MAT  1 20A1 ) 

DO  660  N-1,20 
DO  1000  K=1 ,15 

I F ( COM ( N ) • EG*  MESS ( M ) )  GO  TO'loOO 
CONTINUE 


GO  TO  711 

GOTO  1 9999 . 700 . 600,61 0 . 620.625 . 630 . 635 . 635 ,615.650, 

655,660.5.2000)  .11 

STOP 

CONTINUE 
GO  TO  711 

I? (ITST(l)-l)  6001.6005.6001 
TYPE  6002 

FORMAT  I’  Vi  EIGHT'S  AND  CENTERS  CANNOT  BE  FOUND  UNTIL* 
.*  GEOMETRY  IS  SPECIFIED*) 

CO  TO  20 


IF  1 1 T ST ( 2 ) -1 )  6006 . 6008 ,6006 
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6006  TYPE  600? 

6007  FORMAT! '  CENTERS  CANNOT  BE  FOUND  UNTIL  THJS  WEIGHTS  • 

1  'HAVE  BEEN  I N  PUT '  ) 

6010  TYPE  6011 

6011  FORMAT ('  DO  YOU  WANT  TO  USD  OLD  WEIGHTS?  (YUS  OK  NO)'.'  --  '  ,*) 
ACC LPT  60 t  J 

IF  ( J  •  Nil*  *  Y '  1  GO  TO  20 

CALL  hXii ( WEIGHT .  I DENT *  KUMWG? . 0) 

I73T(2)=i 
CO  TO  20 

6008  N  U  X ! )  U  M =  N  U  :i  C  T  +1;  U  M  F  S  ? 

CALL  •.'GTS <  NUMDUM. WEIGHT • CAPT .  HULLT . CAPMAT. RVMAT  . 

*  C I  NEXT »R INERT) 

ITST (121=1 

GO  TO  20 

610  I?tITST(6)-l)  6110.611.6110 

6110  TYPE  6111 

6111  FORMAT t '  YOU  HAVE  NOT  INPUT  A  LIST  OF  FLOODABLE  VOLUMES' ) 

TYPE  6112 

6112  FORMAT ( '  DO  YOU  WANT  TO  USE  OLD  FLOODABLE  VOLUMES?  (YES' 

1  '  OK  NO) '.'  —  '»$) 

ACCEPT  60.  J 

I F  ( J  •  N  3  • '  Y ' )  GO  TO  20 

CALL  KDV',7  (FLOOD.  IDENl'.NUMFLD.l ) 

611  TYPE  6100 

6100  FORMAT ('  DO  YOU  NEED  A  LIST  OF  FLOODABLE  VOLUMES?  (YES' 

1  '  OK  KOI'.'  —  '.Si 

ACCEPT  60.  J 

I; < J-NE* ' Y ' )  GOTO  6105 

TYPE  6107 

6107  FORMAT  I*  SPACE  NO'.oX.'VOL  ( C U* FT* ) ' . oX, ' LONG  CG'.oX. 

*  '  VERT  CG',/,’  - ' .  6X  > ' - '  ,6X. 

*  ' - ' . 5X . ' - 'i/1 

DO  6103  J»1 , NUMFLD 

6103  TYPE  6i0i  .  J.  ( FLOODUC.  J  )  *.\=1 ,3) 

6101  FORMAT  1 1 H  .3X. I2.10X.F8-1 .7X.F7.3.6X.F7.3) 

6106  TYPE  6102 

6102  FORMAT! 'OHON  MANY  SPACES  DO  YOU  WANT  TO  FLOOD?','  —  '.*) 

ACCEPT  7203* NUMFSP 

TYPE  6106 

6103  FORMAT! '  INPUT  THE  NUMBERS  OF  FLOODED  SPACES'/) 

ACCEPT  6104 « ( IFLDSP ( J ) , J=1 .NUMFSP) 

5104  FORMAT (101) 

CALL  CAS  ( WEIGHT  *  NUM'NGT .  FLOOD. NUMFLD »  IFLDSP.  NUMFSP) 

6105  r:'ST(7)=l 

IF ( NUMFSP* EQ-O) ITST (7 1=0 
GO  TO  20 
615  CONTINUE 
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I F < ITS T ( B > )  6160,6153,6150 

6150  TYPE  6151 

6151  FORMAT  (lii  .’VEHICLE  HAS  ALREADY  SEPARATED  * ) 

6153  IF 1  ITST  1 1 1 ~1  t  ul52.0166.6162 

6152  TYPE  6002 

6156  It (ITST(12»“1)  6157.6150.61o7 
61o7  TYPE  6103 

6153  FORMAT ('  TRIM  ANGLE  CANNOT  BE  FOUND  UNTIL  WEIGHTS  AND* 

1  '  CENTERS  ARE  FOUND* ) 

GO  TO  20 

6159  CALL  HST1 ( ITST ( 7 ) ) 

GO  TO  20 

620  IF (ITST(l)-l)  6200.6202.6200 
6200  TYPE  6002 
GO  TO  20 

6202  I F I ITST (12 ) ”1 )  6157.6204.6157 

6204  I  FI ITST (51-1 )  6205.6207,6205 

6205  TYPE  6206 

6206  FORMAT  I '  YOU  MUST  INPUT  AN  INITIAL* 

1  *  TRIM  ANGLE*) 

GOTO  20 

6207  IF( ITST (7 1-1 >  6208.6210,6203 

6203  TYPE  6209 

6209  FORMAT! *  WARNING— NO  CASULTY  SPECIFIED*  DO  YOU  WANT  TO* 

1  '  CONTINUE?  I  YES  OR  NO)*,*  —  *,S) 

ACCEPT  30, J 

IF! J*NE* *Y* )  CO  TO  20 

6210  CONTINUE 
ITST (3 1=1 
GO  TO  20 

525  I F ( I T  ST  1 3 ) - 1 ) 6250 , 625 2 .6250 

6250  TYPE  6251 

6251  FORMAT! *  CAPSoLE  HAS  NOT  BEEN  SEPARATED  FROM  MAIN  HULL*) 
GO  TO  20 

6252  CA LL  MOT < CAPLEN, CAPWT , VOLCAP . C INERT. RGCAP » SCAP. FLCGCP. 

*  YCGCAP . PLCBCP.PROJCP .NPTS  > 

FL=CA?LEM 

ITST ( 9 ) =1 
ITST(30)=1 
GO  TO  20 

630  I F ( COM  1 N+3) • 2Q- '?* )  GO  TO  9500 

IF ( ITST 1 8 >-l )  6250.6300.6250 

6300  CALL  MOV ( RVLKN . RV NT . VOLKV « RI NERT » RGRV , SHV > FLCGRV , 

*  V  CGR v , FLCBRV . PRO JR V , NPTS ) 

FL=RV  LiN 

ITST(10)=1 
ITST 130) =2 
GO  TO  20 
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635  IF(  ITS':' (9 i  -i  16350.6332 .6360 

6350  IFlITSTUOl-l)  6353.6352. 6 3o3 

6352  GO  TO  (6355.6401 1 »M-7 

6353  GO  TO  16351 .6401 .6451 1 .M-? 

6351  TYPE  6354 

6354  FORMAT ! *  NO  DATA  TO  PLOT*) 

GO  TO  20 

6355  CALL  PLOT! ITS? <301 .NPTS) 

TYPE  2003 

2003  FORMAT ( *  RESET  PLOT  SWITCH  TO  TEXT  MODE  AND*. 

*'  GIVE  A  CARRIAGE  RETURN  --  *.S) 

ACCEPT  30.JJJJ 
ITST (101=0 
GO  TO  20 

6401  IF(ITSTlQ)-l)  6451 .6405.6451 
6405  CALL  PRNT(FL.NPTS) 

GO  TO  20 

6451  TYPE  6452 

6452  FORMAT C  NO  PATH  TO  PLOT”! 

GO  TO  20 

o50  CALL  P DA  i A  ( ITST.F LOOD.  a ~ IG HT .  1* UM «*GT. NUM?  L u . CAPT » HULLT . 

^  DINIT. FLCGIN.VCGITV .FLCGCP. VCGCAP1 
GOTO  20 

655  I F ( V  0  LCA?*o4 • -C APWT 1  6550.6530.6552 

65o0  TYPE  65ol 

6551  FORMAT ( 1  SORRY— CAPSULE  HAS  SUNK*) 

GO  TO  20 

6552  IF (IT ST (11 1-1)  6553.6555.6553 

6553  TYPE  7250 

ACCEPT  7151  .SUP.SEA 
1 SEA=5URSEA 
its:  in  1=1 

6555  IF( IT5T(G)”1 1  6250.6556.6250 

6556  IFi ITST  < 12 1 ~1 )6co7  »65o3.65o7 
6357  TYPE  6559 

6559  FORMAT ( ’  SURFACE  MOTION  CANNOT  BE  PLOTTED  UNTIL  WEIGHTS'/ 

*  '  AND  CENTERS  ARE  FOUND-'/) 

GO  TO  20 

655S  IFI ITSTf 1 1-1 16560.6561 .6560 

6560  TYPE  6562 

6562  FORMAT! '  SURFACE  MOTION  CANNOT  BE  PLOTTED  UNTIL  GEOMETRY'/ 

*  '  IS  SPECIFIED-'/) 

GO  TO  20 

6561  CALL  Vi  VE  ( FLBP.AT .  VCGCAP .  FLEN » BA  REA.  CAPLEN  •  CAP  WT*  FLCGCP.  1  SEA  1 
TYPE  2009 

ACCEPT  30.JJJJ 
GO  TO  20 
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700  U0  710  I=N'+1  »20 
7100  I  F  t  COM  ( I ) • EO • '  *)  00  TO  7110 

710  continue 

711  type  oo 

to  FORMAT*'  ILLEGAL  IK  PUT  COMMAS!) »  THY  AGAIN*  * > 

GO  TO  20 

7110  DO  7111  p'.=l  »? 

I F ( COM ( I +1 ) • EC •  MES  S  2  ( X  J  »  GO  TO  7112 

7111  COST I  So a 
GO  TO  711 

7112  GOTO  <  7 1;> ,  720 » 725 >  730 . 735 . 740 .  ?4o  > . M 
710  TYPE  71o0 

7150  FORMAT!  *  SPECIFY  L/l).  LENGTH.  SEPARATION*  PLASH  —  *.S) 

ACC  HPT'  7151  .FIDRAT.FLEN.SE? 

7151  FORMAT  ( 3F ) 

CALL  C-HOM 
ITSTIl )=1 

CAP  LLi»=i  LEM*SEP 
KV  LEN=F  LEU-CAPLEM 
GO  TO  20 
720  TYPE  7200 

7200  FORMAT  ( *  ISPUT  AVERAGE  CAPSULE  SKULL  THICKNESS  (IN),* 

*  ’  AVER.  AGE' »/« '  HULL  SHELL  THICKNESS  (IN)  —  *.S> 

ACCEPT  7151.CAPT.HULLT  , 

TYPE  7201 

7201  FORMAT ( *  ISPUT  SPECIFIC  WEIGHTS  OF  CAPSULE  MATERIAL.  ASD' ./» 

*  '  OF  HULL  MATERIAL  ( LBS/FT**3)  —  *  .i) 

ACCEPT  7151 .CAPMAT.RVMAT 
TYPE  2002 

2002  FORMAT < '  UO  YOU  WANT  TO  USE  OLD  WEIGHTS?' / 

*  •  YES  OH  SO  --  *  S  J 
ACCEPT  30.  J 

IF( J*NE* 'Y* )G0  TO  2003 

CALL  HDVV;  (  WEIGHT  •  I  DENT » NUMV/GT.Q ) 

GO  TO  2004 

2003  .  TYPE  7202 

7202  FORMAT!  '  HOW  MANY  WEIGHTS  ARE  TO  BE  INPUT?*.  *  '.*) 

ACCEPT  7203»tiUM'.’iGT 

7203  FORMAT ( I ) 

TYPE  7204 

7204  FORMAT <  *  ISPUT  WEIGHTS  IN  LONG  TOSS.  LOKO”L  CG.  VERT  CC* 
i  r  i:i  groups  of  three*/) 

ACCEPT  7161.  (  (WEIGHT (K.J1 )  »K=1  .5) .  Jl=l . iiUMWGT  > 

! : ( J • EO • *  Y ' ) GO  TO  2004 
2005  DO  7210  Jl=i .SUMWGT 
7210  WEIGHT  ( 1 » J1 )  =2240  •♦'WEIGHT  II  »J1 ) 

TYPE  7206 

7205  FORMAT <*  STORE  WEIGHTS?  (YES  OR  SO)'.*  —  *.J) 

ACCEPT  30. J 

I F ( J • EO • *  Y  * )  CALL  ADDOS ( WEIGHT. 0.21 . KUMWGT ) 

2004  IT5T  <2 ) =1 
GO  TO  20 
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IV ( COM  ( 1+2 )  •  Nt>  *  a*  I  CO  TO  7255 
l'YPt  7250 

input  sump a cl  ska  status*  » '  —  • , * j 

ACC KPT  7151,  SUHSKA 
i  S.\A=$  JAS  LA 
ITS'i'dl  1=1 
CO  TO  2C 

if  i co:-: 1 1 -+-2 » . e- * ? •  >  co  to  7ii 

TYP'D  72oo 

fO.iMATC  SPECIFY  INITIAL  VLHTICAL  AKD  iiOiilZONTAL  * 
•SPi'Ki)  (f’i’/SKC)  —  ‘ » $ ) 

A CO KPT  Vial ,Y£LY» VELX 
ITS? ( 3 ) =1 
CO  TO  20 
TYPL  7301 

FOHMA? ( '  SPECIFY  INITIAL  DaPTH  ( FT  I  * , •  —  ',$> 

ACCaPT  7151, DIN  IT 
ITS? <4>=1 
GO  TO  20 
TIPS'  7550 

f  0  A  i  f '  SPaCIFY  INITIAL  THIM  ANGLa  (  DLGRSaS ) ' , '  --  •,$) 
fiCCiPT  7lol,  Tp. IX 
T  ii  I  »Tii  I K/57  •  2S5S 
IT  ST  (5 )  =1 
CO  TO  20 
TYPL  2006 

FO.iMATC  DO  YOU  TAUT  TO  U3S  OLy  FL00DA3L2  VOLUMES?*/ 

'  YDS  OP.  iiO  —  *$) 

ACC KPT  30, J 

IF {J.Ka.'Y*  >00  TO  2007 

CALL  P.DV.M  FLOOD,  IDaiIT,;!uMFLD»l ) 

I  VST (o ) =1 
GO  TO  20 

TV PL  7400 

fO.’i.-.Ait  0  ’1  MA.IY  I  LOOLABLL  VOLij.XLS* , •  •—  ',5) 

ACCaPT  7203, N UMFLD 
TYPa'  7401 

FORMAT ( '  INPUT  FLOOJAiLF  V0LUV2S*  LONG"L,  AND  V2HTICAL* 
Or... . :-.0 IDS' / '  IN  GROUPS  OF  IHRLa'/J 
ACCaPT  7151,  ( < FLOOD (S,J ) ,K=1 ,3) » J=1 ,NJXFLD) 

TYPa  7402 

FORMAT  1  *  ST'Oua  VOLUMaS?  <YLS  Oil  NO)1,*  —  ’,$> 

ACCaPT  30,  J 

I  ?  ( J  •  a<:  •  ’  Y '  )  CA LL  ADDON ( FLOOD  ,1,23,  Nu'MFLD  > 

IT'S?  (61=1 
GO  TO  20 
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745  1=1+1 

GO  TO  7iCO 

2000  m >>;  2001 

2001  rOKMAVC  DO  YOU  DISK  TO  UPDATE  WEIGHTS  OH  VOLUMES? 1  / 

*  ’  TYPE  V  OK  —  •  ,?) 

ACCEPT  30  iJ 

ir  ( J  •  LO  • '  V '  )  GO  TO  1003 
IP  CJ.  KG-  ICO  TO  1903 
GOTO  711 

1999  CALL  UPDATE  (  PLOOD»  1  *  iliiMPLD  i 
I'YST ( 25 ) =1 
G  0  TO  20 

105B  CALL  UPDATE!  v.'EIGHT»0*NUM-.'GT) 

ITST(25)=i 
GO  TO  20 

9500  CALL  A ST MLS 
GOTO  20 
LHD 


IX- 2 . 


Subroutine  GE0 M 


The  volume,  surface  area,  volume  moment  of  inertia, 
radius  of  gyration,  longitudinal  center  of  bouyancy,  and 
the  longitudinal  center  of  gravity  of  the  shell  are  cal¬ 
culated  in  the  subroutine  GE0>3 .  The  subroutine  accepts 
the  values  of  the  separation  point,  length-to-diameter 
ratio  and  vehicle  length.  The  non-dimensional  variables 
defining  the  geometry  of  the  intact  vehicle  IR1=0,  Sl=2] 
are  set  and  the  program  passed  to  the  basic  geometry 
calculation  routine  where  non-dimensional  values  for  the 
required  parameters  are  calculated.  The  program  is  then 
returned  to  Statement  1,  where  these  values  for  the  intact 
vehicle  are  saved.  The  non-dimensional  variables  de¬ 
fining  the  capsule  [R1=0,  Sl=  the  separation  point]  are 
then  set  and  the  program  is  returned  to  Statement  2  where 
these  values  for  the  capsule  are  saved.  The  non-dimen¬ 
sion  variables  defining  the  remaining  vehicle  ]R1  = 
separation  point,  SI  =  2.0]  are  set  and  the  program  again 
passed  to  the  basic  geometry  calculation  routine.  When 
the  required  parameters  for  the  remaining  vehicle  have 
been  calculated  the  program  then  passes  to  Statement  30 
where  the  required  dimensional  parameters  are  calculated. 
The  program  then  returns  to  the  monitor. 

This  subroutine  contains  no  error  messages. 


Input  Variables 


XLEN 

SEP 

XLD 


Length  of  intact  vehicle  (ft.) 

Separation  point  from  nose  (decimal  fraction) 
Length  to  diameter  ratio 


Output  Variables 


C 

C0 

D 

CGCAP1 

CGRV1 

V0LITV 

LCBITV 

SITV 

V0LCAP 

LCBCAP 

SCAP 

RGCAP 

VOLRV 

SRV 

RGRV 

BAREA 

PR0JCP 

PR0JRV 


Longitudinal  center  of  buoyancy /XLEN,  intact  vehicle 
Longitudinal  center  of  buoyancy /XLEN ,  capsule 
Longitudinal  center  of  buoyancy /XLEN ,  remaining  vehicle 
Center  of  gravity /XLEN,  capsule  shell 
Center  of  gravity /XLEN ,  remaining  vehicle  shell 
Volume  of  intact  vehicle  (ft.3) 

Longitudinal  center  of  buoyancy  of  intact  vehicle 
(ft.  from  nose) 

Surface  area  of  intact  vehicle  (ft.2) 

Volume  of  capsule  (ft.3) 

Longitudinal  center  of  buoyancy  of  capsule  (ft.  from  nose) 
Surface  area  of  capsule  (ft.2) 

Radius  of  gyration  of  capsule  (ft.) 

Volume  of  remaining  vehicle  (ft.3) 

Surface  area  of  remaining  vehicle  (ft.2) 

Radius  of  gyration  of  remaining  vehicle  (ft.) 

Base  area  of  capsule  and  remaining  vehicle  (ft.2) 

Center  of  pressure  of  projected  capsule  area/XLEN 
Center  of  pressure  of  projected  remaining  vehicle 
area/XLEN 


IX-2.1  Flow  Chart  for  Subroutine  GE0M 
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l>t  'X/OC  G e,  d> m 


Best  Available  Copy 


IX-2.2  Listing  of  Subroutine  GE0M 


SUBROUTINE  GEOM 
REAL  LCB ITV » LCBCAP>  LCBRV 
COMMON/ INPUT/XLEN  t SEP #XLD*  SKIP<4) 

COMMON/ INOUT/C. CO. DiCGCAPl*CGRVl . VOLITV .LCBITV. SITV. 
*VOLCAP*LCBCAP»SCAP*RGCAP«VQLRV»  SRV»RGRV«BAREA*PRQJCP» 
*PROJRV. EXTRA (716) 

T1  «»  2**S£P 
XLO=  1/XLD 
C 

C  CALCULATION  FOR  INTACT  VEHICLE 

C 

XLSQ=XLO*XLO 
R1  =  0-0 
SI  =  2.0 
GO  TO  10 

1  VO  =  VI 
C  =  XI 
B1  =  A1 
G1  =  XM1 

FI  =  SORT  (XJ1) 

CGITV1  =  CG 
C 

C  CALCULATION  FOR  CAPSULE 

C 

R1  =  0 
SI  *  T1 
GO  TO  10 

2  D1  =  VI 
CO  =  XI 
S  =  A1 
HI  =  XM1 

F  =  SORT  (XJ1) 

CGCAP1  =*  CG 
PR0JCP=CGPR0J 
C 

C  CALCULATION  FOR  REMA I NRN4HUGHICLE 

C 

R1  =  T1 
SI  =  2* 

GO  TO  10 

3  El  *  VI 
D  =  XI 

Cl  =  A1 
XII  =  XM1 
FO  =  SORT  (XJ1) 

CGRV1  =  CC 
PROJRV=CGPROJ 
•GO  TO  30 
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c 

;;  BASIC  GEOM.  CALC*  PACKAGE 

C 

10  PI  *  3*14159 

VI  «  IP1/24* )*XLSQ*(3.+( Sl-Rl 1+tRl-l* )**3»( Sl-1* )**3) 

XI  «  (2./3.)*<Sl+*3-Rl**3>-<l./4.)*<Sl'M‘4-Rl**4) 

XI  *»  (XI 1/t  <2*/3» )*(3«*( Sl-Rl )  +  (Hl-l* )**3-( Sl-1* )**3) ) 
XM1=  (Pl/320*  )*XLSQ*(5«*(  Sl**4-Rl**4)-2**(  S1**5-B1 
1  **5 1 ) 

XJ1  =  (XM1/V1 )  -X1*X1 

A1  =  ( ( Sl-1 • 1/2- >*SQRT(XLSQ+2.*<1— XLSQJ*S1-(1— XLO 

1  *XL0M‘S1*S1  )  +  (l./(2.*SQRT(l— XLSQ) )  )*ASINt  <  (Sl-l.)*(l  — 

2  XLSQ) )/(SQRT(l— XLSQ) ) )-< (Rl-l»}/2* )*SQRT(XLSQ+2* 

3  *<1— XLSQ)*R1-(1— XLSQ)*Rl*Rl)-<l./I2.*SQRT(i— XLO+ 

4  XLO) )  )*ASIN(  ( (Sl-1*  )*(1— XLSQ)  >/(SQRT<l— XLSQ) )  > 

AREA=A1 

A1  =  <Pl/2» )*A1/XLD 
C 

CGPR0J=PR0J(S1.R1) 

C  CALCULATE  C*G*  OF  SHELL 

C 

25  CG=( ( 1 XLSQ+2- * ( 1  — XLSQ ) *S1- ( 1  — XLSQ )*S1*S1 ) **1 *5) 

*  -  ( ( XLSQ+2  •  * <  1  •  -XLSQ )  *R1-  <  1  •  -XLSQ  > *R1*R1  > **1  •  5 ) )' 

*  /<3.*(XLSQ-1*>*(AREA>) 

CG=(CG/2»)+.5 

IF(R1)20.20.3 

20  IFCS1-2*)  2.1.3 

C 

30  CONTINUE 

C  OUTPUT 

C 

VOLITV  =  V0*XLEN**3 

LCBITV  =  C*XLEN 

SI TV  =  B1  *XLEN*XLEN 
RGITV  =F1  *  XLEN 
CGITV  *  CGITV1  *  XLEN 
VOLCAP  =  D1  *  XLEN  **3 
LCBCAP  =  CO  *  XLEN 
SCAP  =  S*XLEN*XLEN 
RGCAP  =  F*XLEN 
CGCAP  =  CGCAP1*XLEN 
VOLRV  =  E1+XLEN**3 
LCBHV  -  D*XLEN 
SRV  a  C1*XLEN*XLEN 
RGriV  =  FO*XLEN 
CGRV  =  CGRV1*XLEN 

BAREA= ( ? 1 *  <  X  LO ) **2* <  2 • *T 1 -Tl+Tl ) *XLEN*XLEN ) /4 « 

TYPE  500. VOLCAP. LCBCAP. CGCAP 

500  FORMAT ( ’  CAPSULE  VOLUME  «*.F12.2./ 

*  •  CAPSULE  LONG.  CENTER  OF  BOUYANCY  -’.F10.3./ 

*  *  HYDRODYNAMIC  CENTER  OF  CAPSULE  .FIO'.'S.//) 

RETURN 

END 
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FUNCTION  PROJ(X.Y) 

Q=2**X“X*X 
P=X-1 • 

P1=Y~1 • 

Q1=2*+Y-Y*Y 

Z=.5*?*Q**.5+ASINf?) 

Z1“»5*P1*01**»5+ASIN(P1 ) 

PROJ=  (  t  <  “0**1  *5/3+2 )  - 1  -Ql**l  .5/3+Z1 ) )  /( Z-Zl J>  **5 
RETURN 
END 
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IX-3 . 


Subroutine  WGTS 


As  inputs  the  subroutine  accepts  the  number  of 
weights,  the  weight  array,  the  separation  point,  the 
vehicle  length,  the  length- to-diameter  ratio,  the 
thickness  of  the  capsule  and  remaining  vehicle  shell 
material  as  well  as  the  longitudinal  center  of  gravity 
of  the  capsule  and  remaining  vehicle  shell. 

The  program  then  calculates  the  total  weight  of 
the  intact  vehicle,  the  capsule,  and  the  remaining 
vehicle  components.  Subsequently  the  non-dimensional 
longitudinal  and  vertical  centers  of  gravity  of  the 
intact  vehicle,  the  capsule,  and  the  remaining  vehicle 
components  are  computer.  The  shell  and  plate  effects 
are  then  added  and  the  moments  of  inertia  of  the  capsule 
and  remaining  vehicle  are  calculated.  Finally,  the 
dimensional  results  are  determined  and  the  program  returned 
to  the  monitor. 
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Input  Variables 


10 

G 

TCAP 

TRV 

MATCAP 

MATRV 

SEP 

XLEN 

XLD 

CGCAPl 

CGRV1 

SCAP 

SRV 


Output 

T0TWT 

RVWT 

CAPWT 

LCGITV 

VCGITV 

LCGRV 

VCGRV 

CAPM0I 

RVM0I 


Number  of  weights  in  distribution 
Weights  and  locations  array  (lbs.,  dec.  frac. , 
dec.  frac.) 

Thickness  of  capsule  shell  (in.) 

Thickness  of  remaining  vehicle  shell  (in.) 
Specific  weight  of  capsule  material  (lbs/ft.3) 
Specific  weicht  of  remaining  vehicle  material 
(lbs. /ft. *) 

Separation  point  from  nose  (decimal  fraction) 
Length  of  intact  vehicle  (ft.) 

Length  to  diameter  ratio 

Center  of  gravity /XLEN ,  capsule  shell 

Center  of  gravity /XLEN,  remaining  vehicle  shell 

Surface  area  of  capsule  (ft.2) 

Surface  area  of  remaining  vehicle  (ft.2) 


Variables 


Weight  of  intact  vehicle  (lbs.) 

Weight  of  remaining  vehicle  (lbs.) 

Weight  of  capsule  (lbs.) 

Longitudinal  center  of  gravity  of  intact  vehicle 
Vertical  center  of  gravity  of  intact  vehicle 
Longitudinal  center  of  gravity  of  remaining  vehicle 
Vertical  center  of  gravity  of  remaining  vehicle 
Moment  of  inertia  of  capsule  (slug  *ft.2) 

Moment  of  inertia  of  remaining  vehicle  (slug  •  ft.2) 


IX-3.1  Flow  Chart  for  Subroutine  WGTS 
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IX-3.2  Listing  of  Subroutine  WGTS 
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C)  CJ 


:•  .-.Ou  i  i.  ,\  ..  GVO  (  i  0 ,0  »  T  CA  t  »  »  a  v  .  I’.A ,  CAft  .'■! A  t  n»  »  Cnl'XOl  •  n V  KOI  ) 

I'  l  ii  ci .  «. a L  <  K » L ,  ) 

COW.  UN/ 1  TruT  /X  !■*• . » Ufai* •  X Lu»  VO*  14*  Yo*  Yo  - 

COMMON/  ilaUJi  f  I  ait  »C  GCAT!  ♦  t  ialt  1  *  i>UN<  <U)  *  .iCAPt  X IT ( 2 1 

♦  PAGO  14 1  > . CT .  *  a »  .« i  t CAr;.  t  •  LCGlT  Y.VCGiii,  LCGnV • i  CGuV  •  L CGCAP, 

*  i  Caf  CA  P  *  La  i  A  (  /O  A ) 

0 

J  I  WT  .1  I  I  >  <  (0  1  ,0  (  aGOO  I  t  <  L  <  0  1  t  :1  I  a5U0  ) 

C  CA  LCUiaA  i  10. a  Of  TOTAL  r.  [  a#  r.  1 

C 

CAP  J'0 1=0*0 
;(v.'vGI=G*0 

:;0  10  .1  ;:=i,  0 

x  ( j  i  =o  • 

io  k  n  j  i  =o  • 

J0=3*IO 

ilG  20  I=i,J0*15 

ill  1+2  1=0  ( 1+2 l  +  U/XLO)  *(0(1+1  1-0 ( I+l)=*2)*+.b 
20  KID*  0(I)+K(1I 

CALCULATION  Of'  CAPGULL  aifclGHT 

00  00  i=l.J0,G 
If  ( SLP-G < i+1 1 1  00. >50. 20 
2o  K ( 0 )  =0(1 )+K  <01 

00  CO XT  IT CL 

C  CALC  o  LA I  Oil  Of  liLKAINING  vLilCLL  ..'LIGHT 

C 

K(2)  =  K  ( 1  1  -  X  ( 3 1 

C  CALCULATION  Of  vLAT*  A  No  LONG  •  CO  Of  INTACT  VLHICLL  LXCLPT 

C  SilaLLL 

If ( K ( 1 ) 1100,100.110 
100  .<(41=0* 

K ( b 1 =0 • 

00  TO  120 

no  oo  40  1=1, jo. o 

c 

K  <4  )  =  0(l)=  0(1+11/  .<(11  +K  (4 ) 

40  K  ( b  I  =  0(1  )*  .4(1+21/  Kill  +K(o) 

C  CALCULATION  Of  VLAT*  ANO  LONO*  CO  Of  CAPSULL  LXCLPT 

C  SiiKLL 

C 

120  If'IKIb)  1100.  130.  140 

130  K(yi=o* 

K ( 0 1 =0* 

00  TO  loo 
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“  i  «  ■ 


4*tJ  .  )  -.1 

if  -)(!+ll)  >.;»■;;> 

40  riiiii  -  i  it  „ (  1  *i  l/jt(3)  + 

;:(./»  *>  ju)  *  ti<  i+2»/\Mt+  fitoi 

Cr. rMO  I  *0 A  VXOi ( 1  I  *  (0  t  I  +1  )  ++2+a  (  1+2 )  **2  1 

ao  coi>.i!iur. 

c 

C  Ur.  I>i>  J  Lit',  i  *  I  0.»  0*  *  r.:\  V*  La.'.O*  C  0  Of  ;{£MAltiIN0  VKhlCLf! 

C  f.f'OM-i  Olti.l.i. 

loO  I .-'(  K  (2  )  )  100»160»17u 
ICO 

»  m so¬ 
il  O',  o  1  1  vj 

c 

IV 0  i;0  00  i-l»J0*O 

I  :  (0(1+1)  *3>.V)  CO.OOtOa 
bo  KIO)  =  0(  1  )*0(  1+1  )/X(2)  +  K  ( 6 ) 

K  ( V )  =  0(1)*  n(l+2)/X(2l  +  K  (V) 
k  V  *00  [  V  ;'.01+0  ( I  )  *  (  0  <  1+1  )**2+Hl  I  +2 )  **<£ ) 

60  CO!.'.  I!. or. 

C 

C  CALCuLAvIO..  o:  Sill.LL  at'fiCVS 

C 

ISO  VI  =  o. 14100 

Vi  =  2*  v.Skr1 

A iA  =  PI  *  ( ( 1  ~  ( Vl  -1  )*+2  *  I  *XI.ri:(*I  Liii. )  /  ( (  XLu*XLi>)  *4* ) 

■  PC/.P  =  |  r.n  *'f  0  A  P*X  A'i’C  At*/lJi» 

=  A  ;i  '.A  *  V V  *  On  V  k  V  / 1  2  * 

■\  SOAP  =  3CAP*TGA?+MAvCAiV12« 

;. Sr.V  =  0:-.  r  *V.;V  *;-!AV  r.V/12  • 

..  VOV  =  K(l)  +  .vpCAr  +  .iPuY  +.’.00  A  V  +  ..S:w 
K  (4  )  =  (  K  (4  )  *«i  ( 1  )  +  < OPCAP+vPn V  )  *SrY+  (  ,iSCAP*CCCAPl  )  +  (  Won v* 
i  COr.vl  i  l/OiOV 

I  r  (  K  ( .^ )  +  YCA  P+..SCA P ) 200 » 200  ?210 

210  K  (  o)  =  (K(.0)*.<(3)+..P0AV*00i'+..SCA?X'C0CAPl )  /(  K  ( 3 1 +.JPCA?+.)3CAP ) 
200  I K  ( !\  ( 2 )  +  ..  PitV+.’iSit  V  )  220 » 220  »23u 

200  Klo)  =  (rvlu)*K(2 1  +  .VPA  v  *se.2+.i5ii » +C0KV1 1  /  (  K 1 2  )  +  iHp:{V+.iSrtV  ) 

220  Kill  =  ..  i  0 . 

X  (2)  =  X  (2 )  +  +  vJSnV 

K  ( 3 )  =  K(3)  +  .\PCAP  +  0 OCA  V 

r 

V/ 

C  ALL  0?  ArfOVri  CO'S  A:ti)  ri.it  uiil'V  LtNOVr!  BLLO.'i  .ii;  CALC 

C  iJ I S  IOiiAi.  iii.SuLVS 

C 

Cnr::Ol=(CA2:!Ol  +  .,.SCAr*COCA21*CCCAPl+WPCA2!t'SKP*SL?“K(3) 

+•*  ( ( tj)  **2+K  ly )  **2  )  )  +/ LKO+XLc-M 


i 

0 

'{ 

i 

i 

i 

i 

i 


i 

i 
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( if »  KOI+.ioit  V  II 1  *Ciia  if  1+ ■tt’CAi'*  ( 2 ) 

**(  X  <u  I **2+K  I  *xun*xu:< 


do  vo  I =4 i y 
<  L( i )  =  k<i> 

VO  K 1 1 )  =  K ( I )  *  XLtN 
i-O-i'  .'.T=K  ( 1 ) 
i(V:,'i=KU> 

CAP.w  =«<;•>  I 
[,0iIi'V=KL(4l 
VWITV«KL(a) 
LCOj'.VsKHCj) 

VCGn  V=KL( V ) 
i<CG  CaP“K  LIB) 

V  COCA  :'“=K  L  ( 9 ) 
iitTUaN 

e;id 
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IX- 4.  Subroutine  CAS 

The  CAS  subroutine  is  called  by  the  coonand 
DAMAGE)  and  simulates  flooding  the  specific  flood able 
volume (a)  with  sea  water.  These  new  weights  are  then 
added  to  the  present  distribution  along  with  their 
respective  logitudinal  and  vertical  oenters. 

There  are  no  error  messages. 


Input  Variables 


WEIGHT 

NUMWGT 

FL00D 

NUMFLD 

IFLDSP 

NUMFSP 


Weight  array  (lbs.,  dec.  frac.,  dec.  frac.) 

Number  of  weights  (triples)  in  array,  WEIGHT 
Floodable  volume  array  (ft.1,  dec.  frac.,  dec.  frac.) 
Number  of  volumes  (triples)  in  array  FL^0D 
Array  containing  identification  of  spaces  to  be  flooded 
Number  of  spaces  to  be  flooded 


Output  Variables 

WEIGHT  Extended  weight  array  (lbs.,  dec.  frac.,  dec.  frac.) 
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IX- 4. 1  Flow  Chart  for  Subroutine  CAS 


-Sju:-0uV1  K  CAS  I KiGr-T  •  i.  Ji-i  .«iV»  KI.OOU,  liuMUO.  Il'l0Si>.NUM*'Si' ) 

.-.tiioriviaiioo)  *i'LOoma.ioo) * ivu>sp(10) 

*jO  100  J  =  1  * 

JslijM/.JV+I 

..SJariill  t  J  )  =cj4 •  *r'LOGi; ( l ,  Ir'LOSrM  I )  ) 
iiSICil'i' (2>  J)  =*’1,000  (2.  IFLOSP  ( I )  i 
100  sBIt»ri'i’<3»J ) •FLOOD (3*  ( I ) ) 

;<  KVuitr: 
c.UO 
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IX- 5. 


Subroutine  HST1 


Inputs  for  the  hydrostatics  program  include  the 
weight,  volume,  longitudinal  and  vertical  centers  of 
gravity,  longitudinal  center  of  buoyancy,  and  initial 
trim  of  the  intact  vehicle.  If  the  vehicle  is  damaged 
(IDMSW=0) ,  calculations  are  made  in  the  first  half  of 
the  program  primarily  between  statements  2000  and  2080. 
If  the  vehicle  is  damaged  (IDMSW=1) ,  the  calculations 
are  performed  between  2110  and  2200. 

The  subroutine  is  called  in  response  to  the 
monitor  command  TRIM}  to  evaluate  the  hydrostatic 
characteristics  of  the  undamaged  or  damaged  intact 
vehicle  designed  by  the  user.  It  outputs  no  data  for 
use  by  any  other  subroutine. 

There  are  no  error  messages  as  such;  however,  the 
messages  printed  out  are  meant  as  an  intuitive  design 
guide  for  the  user. 


Input  Variables 

IDMSW  Damage  switch 

FLCBIN  Longitudinal  center  of  buoyancy /L,  intact  vehicle 

V0LITV  Volume  of  intact  vehicle  (ft.3) 

WTINT  Weight  of  intact  vehicle  (lbs.) 

FLGGIN  Longitudinal  center  of  gravity/L,  intact  vehicle 

VCGIN  Vertical  center  of  gravity/L 


Output  Variables 

BFNBD  Net  buoyant  force  of  intact  vehicle  before  damage  (lbs.) 

TT  Bow  down  trim  of  intact  vehicle  (deg.) 

T  Bow  up  trim  of  intact  vehicle  (deg.) 

BFNAD  Net  buoyant  force  of  intact  vehicle  after  damage  (lbs.) 

TAD  Bow  down  trim  of  damaged  intact  vehicle  (deg.) 

TADIN  Bow  up  trim  of  damaged  intact  vehicle  (deg.) 


IX-5.1  Flow  Chart  for  Subroutine  HST1 
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SUBROUTINE  iiSVl  ( IUMSW) 

COMMON  / 1 NPUT/  FLOA .  PAS  S  ! 2 ) .  H 

C0MM0N/IN0UT/FLCBIM»SKIP!4).V0LITV.STUFF!12).WTINT. 

*  FOOL !  2 ) « F LCG I N .  VCG I N »  EXTRA (711) 
C*++*******it>iiitt,(>*,t',t'*,t‘i«*,*‘*******+*,tt 

c 

C  CALL  TO  OVERFLOW (“I )  SUPPRESSES  DIVIDE  CHECK  ERROR  HESSAGKS* 

C 

C******************************** 

CALL  OVERFL(O) 

B0YINT=64«*V0LITV 
GO  TO  (2110) * IDMSrf 
BFNBD=BOYINT-WTINT 
B=ABS(BFNBD> 

I F ( BFNBD ) 2000 .2010.2020 
2000  TYPE  2990. B 

2990  FORMAT (36H  INTACT  VEHICLE  IS  HEAVY  OVERALL  BY  .F10«0. 

1  7H  POUNDS) 

CO  TO  2030 
2010  TYPE  2980 

2980  FORMAT (36H  INTACT  VEHICLE  IS  NEUTRALLY  BUOYANT) 

GO  TO  2030 

2020  TYPE  2970. BFNBD 

2970  FORMAT (36H  INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  .F10»0. 

1  ?H  POUNDS) 

2030  IF(ABS(FLCGIN-FLCBINI-.0001 >3000.3000.2031 
3000  IFUCGIN-.  0001 13010.3010.2031 
3010  TYPE  3020 

3020  FORMAT*  *  CG  AND  CB  COINCIDENT.  TRIM  ANGLE  INDETERMINATE’/) 
RETURN 

2031  T=57*2968*ATAN2! (FLCGIN-FLCBI N ) . (-VCGIN) ) 

TT=ABS<7> 

2040  IF! IFIX(T*100* ) 12050.2060.2070 
2050  TYPE  2960. TT 

2960  FORMAT (36H  INTACT  VEHICLE  IS  TRIMMED  BOW  DOWN  »F7*2» 

1  8H  DEGREES)  r 

RETURN 

2050  TYPE  2950 

2950  FORMAT <27H  INTACT  VEHICLE  HAS  NO  TRIM) 

RETURN 

2070  TYPE  2940. T 

2940  FORMAT ( 34.4  INTACT  VEHICLE  IS  TRIMMED  BOW  UP  .F7.2. 

1  ’DEGREES’) 

2080  RETURN 

THE  ABOVE  ESTABLISH  THE  BEFORE  DAMAGE  CONDITION  OF  THE 
C  INTACT  VEHICLE 
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2110  BFNAD=BOYINT-WTINT 
B=ABS(BFNAD> 

I? ( IFIXIBFNAD+100*  t  >2120.2130. 2140 
2120  TYPE  2920. B 

2920  FORMAT <44H  DAMAGED  INTACT  VEHICLE  IS  HEAVY  OVERALL  BY  « 

1  F10.0.  7H  POUNDS) 

GO  TO  2150 
2130  TYPE  2910 

2910  FORMAT (44H  DAMAGED  INTACT  VEHICLE  IS  NEUTRALLY  BUOYANT) 

GO  TO  2150 

2140  TYPE  2900.BFNAU 

2900  FORMAT I 44H  DAMAGED  INTACT  VEHICLE  IS  LIGHT  OVERALL  BY  » 

1  F10*0»  7H  POUNDS) 

2150  IF(ABS(FLCGIN-FLCBIN>--0001 >4000.4000,2160 
4000  I F ( VCC I N- • 0001 14010. 4010.2160 

4010  TYPE  4020 

4020  FORMAT) *  CG  AND  CB  COINCIDENT.  TRIM  ANGLE  INDETERMINATE’/) 

RETURN 

2160  TAD I N=57 • 2958*ATAN2 < (FLCGIN-FLCBIN) , (-VCGIN) ) 

TAD=ABSITADIN ) 

IF( IFIX (TADIN+100* ) 12170.2180,2190 
2170  TYPE  2890. TAD 

2890  FORMAT (54H  DAMAGED  INTACT  VEHICLE  WILL  EVENTUALLY  TRIM 
1  BOW  DOWN  ,  F7*2,8H  DEGREES) 

CO  TO  2200 
2180  TYPE  2380 

2880  FORMAT (50H  DAMAGED  INTACT  VEHICLE  WILL  EVENTUALLY  REACH 
1  EVEN  KEEL) 

GO  TO  2200 

2190  TYPE  2870.TADIN 

2870  FORMAT <52H  DAMAGED  INTACT  VEHICLE  WILL  EVENTUALLY  TRIM 

1  BOW  UP  ,  F7*2.8H  DEGREES) 

2200  TYPE  2860.  T.TADIN 

2860  FORMAT I47H  INTACT  VEHICLE  WILL  TRIM  HYDROSTATICALLY  FROM 

1  .F7.2./.26H  DEGREES  BEFORE  DAMAGE  TO  ,F?*2.14H  DECREES  AFTER 

2  *  DAMAGE’) 


2210  RETURN 
END 
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IX-6.  Trajectory  Solution  (M0T) 

This  subroutine  calculates  the  trajectory  of  either 
the  escape  capsule  or  the  remaining  vehicle  depending  on 
whether  it  was  called  by  typing  the  command  ASCEND)  or  ' 

the  command  HULL) , 

The  subroutine  allows  the  user  to  specify  the  upper 
limit  of  the  time  interval  in  seconds.  This  value  is 
non-dimer.sionalized  internally  and  calculations  proceed 
for  evaluating  the  distances,  angles  and  velocities  of 
the  vehicle  from  the  point  of  separation  subject  to  the 
initial  conditions  of  speed,  trim  and  physical  characteristics 
of  the  vehicle  input  or  calculated  in  the  monitor  or  in 
preceding  subroutines. 

Error  Messages 

1)  "THE  VEHICLE  IS  DYNAMICALLY  UNSTABLE ...  INTEGRATION 
HAS  BEEN  TERMINATED.  RESULTS  MAY  BE  PRINTED  OR  GRAPHED."  - 
An  instability  has  been  detected  and  time  histories  have 
been  terminated.  The  histories  may  be  printed  or  plotted. 

Corrective  Action:  Change  input  parameters  to  try  to  obtain 
stable  solution  if  desired. 
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XLEN 

WT 

FV0L 

FINERT 

RAD 

SAREA 

FLCG 

FVCG 

FLCB 

FLHC 

BAREA 

Y3 

Y4 

Y5 


Length  of  vehicle  (ft.) 

Weight  of  vehicle  (lbs.) 

Volume  of  vehicle  (ft3.) 

Moment  of  inertia  of  vehicle  (slug*  ft^.) 

Radius  of  gyration  of  vehicle  (ft.) 

Surface  area  of  vehicle  (ft2.) 

Longitudinal  center  of  gravity  of  vehicle/XLEN 
Vertical  center  of  gravity  of  vehicle/XLEN 
Longitudinal  center  of  buoyancy/XLEN 
Center  of  pressure  of  projected  vehicle  area/XLEN 
Basic  area  of  vehicle  (ft^.) 

Initial  trim  angle  (deg.) 

Initial  vertical  velocity  (ft. /sec.) 

Initial  horizontal  velocity  (ft. /sec.) 


Output  Variables 


X 

z 

D 

E 

V 

G 

W 

NPTS 


Time  array,  non-dimensional 
Attitude  array  (deg.) 

Vertical  distance  array,  non-dimensional 
Horizontal  distance  array,  non-dimensional 
Angular  velocity  array  (deg. /non-dimen. time) 
Vertical  velocity  array,  non-dimensional 
Horizontal  velocity  array,  non-dimensional 
Number  of  points  in  each  array 
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IX-6.1  Flow  Chart  for  Trajectory  Solution  (M0T) 
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■  »■  *  t*  *  «lw:  at  r<  t&  .■ 


SUBROUTINE  M0T<XLENtVT>FY0L»FINEKTtRGtSARKA»FLCQtFVCG»fLCB*7UiC* 
*NPTS) 

REAL  HI *H2*  M3* II 

COMMON / I NPUT/FLEN t  XO • FLO RAT » T3* Y4  *  Y5*  Y6 
COMMON /RUNG/ I NIT IA 

COMMON / 1 NQUT/ EXTRA  1 15) •BAREAtSKIPlllI »XI101 l tZ(101 ) * 

*D<101 >.EE< 101 ».V<1011»C<101» *¥1101 > 

DIMENSION  F(6)tY<6) 

INTEGER  T.C.Q.RUNGE 
KOUNT-1 
INITIA-0 
PI-3.14159 
30  N=6 

TYPE  3000 

3000  FORMAT ( '  INPUT  MAXIMUM  SIMULATION  TIMS  IN  SECONDS  —  *»> 

ACCEPT 3001 tX9 

3001  FORMAT (F) 

X9-X9/1XL£N/32.17)**.5 

H=X9/100. 

HH-H 

10  IF(*5001.GE.HH)C0  TO  60 

HH=HH--5 
K0UNT=K0UNT+1 
GO  TO  10 
60  0=1 

KKNT-KOUNT-l 

H=H/FL0AHK0UNT» 

T=1 

M3=0*2 

M1=2.*WT/(32«17*1. 9&*XLEN*BAREA ) 

H2=FY0L/<XLKN*BARBA» 

D1=FLCG/X0 

B2=FLCB*FLKN/XLEN 

D3=FLHC*FLEN/XLEN  _  _ _ 

RHO-1.995 

Y1=0* 

Y2=0. 

Y6=0* 

D0L=SQRT(BAREA*4«/3. 14159 l/XLEN 
E-SQRTU  — <D0L>**2» 

AL=2.*<1.-E*E)*<.5*AL0GH1.+E)/<1— El)-E»/E**3 

B  E*1 • / ( E**2 ) - < 1 • -E**2 1 *A  LOG l(l*+£)/(l*~EI »/(2.*E**3» 

FKA=AL/(2«-AL+ 

-FKC=BE/(2.-BE> 

FKTHET=E**4*(BE-AL>/H2.-E**2>*12.*a-a 

E**2-(2.-E**2)*«BE-ALm  . 

1 1-2 . *P INERT/ ( RH0*XLSN**3*SARSA ) *2 • *M2* 1 1 RC/XLSM  > **2 1 *FKTHET 
XX-O.  _ _ 


© 

© 

© 

o 

o 

o 

o 

o 

o 

~3 

o 

o 

O 

o 
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r 

360 

Y(1)-Y1 

?  ( 

370 

YI21-Y2 

V 

380 

Y(3)-Y3 

1 

390 

Y  f  4 ) -Y 4/ SORT < 32 • 1?*XLEN ( 

? 

400 

Y ( 5 ) -Y5/S0RT < 32 • 1?*XLKN > 

r. 

410 

Y(6)-0* 

430 

C-0 

■  o 

510 

DO  8  JNT-1.101 

J JHT- t JNT*KOUHT > - C  KXNT  > 

IF(T«LT. JJNT)  00  TO  555 

o 

IF<T*EQ* JJKTIGO  TO  520 

8 

CONTINUE 

520 

XIJNTI-XX 

o 

EEIJNT)-Y(l) 

525 

D(JNT>«Y(2) 

c 

530 

Z(JNT) 

«Y<3>*67.2958 

o 

635 

W(JNT)=Y<4> 

540 

G<JNT)=Y<5) 

545 

V ( JNT  >  =Y  <  6 )  *67  >2958 

c 

550 

NPTS«JNT-1 

565 

T=T+1 

560 

CALL  RUNGE(6.Y.F.XX»H.K) 

r. , 

580 

IF<K-1 >750.600.750 

( 

600 

F(1)=Y<4> 

610 

F ( 2 ) =Y  <  5 ) 

620 

F  <  3  >  =»Y  ( 6  > 

Cl=l« 

IF(Y(4).GT*50.)C0T0  2000 

- 

IF(Y(5).GT.50*>G0  TO  2000 

621 

VX=SQRT ( 32 • 17*XLEN* ( Y  <4  > **2*Y (51**2) > 

IFIVX--001 >633.633.622 

622 

Rl=VX*XLEN/l*0E-5 

624 

FDRAG=(3*46*AL0G(Rl>-5.6>**<— 5i 
*+C • 0016*  C  XLEN/SQRT  <  BAR£A*4 • /3* 14159  > ) *R1 ** ) -2 • /5* ) 

625 

DDRAG=*0875714-0* 119702*D0L+0» 1 33631 *D0L**2 

626 

BCOEF=FDRAG*SAREA/BAREA 

627 

TES=XLEN/FLEN 

628 

IF  (TES»LE*0«5)  GOTO  631 

629 

BDRAG=0 • 029* ( BCOEF** ( - .5  > 1 * <  ( DOL*XLEM > /(FLKM/XOS J 

- 

630 

GOTO  632 

631 

BDRAG=0 • 029* ( BCOEF** (“ «5 ) ) 

v 

632 

C1=FDRAG+DDRAG*BDRAG* ( BA  REA/ SAREA  > 

633 

CONTINUE 

EN1=FL0AT( INT(ABS<Y(3> )/(2**3*14159) I » 

(_ 

IF< Y (3> *LT»0« )  GO  TO  650 

( 

EN2=Y(3)-EH1*2.*3. 14159 

GO  TO  660 

650 

EH2»Y(3)+EN1*2. *3.14159 

660 

IF(Y(5>.EQ-0.0.A»D.Y(4>.GE‘0.0)  A-EN2 
IF<Y<5>.GT.0*0.AMD.Y<4>»£Q*p»0>  A-EM2-PIA2. 
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IF!Y(5)*LT*0*0*AND*Y<4>*EQ*0*C»  A-EN2+P1/2* 
IF(Y(5)*EQ*0*0*AND*Y!4)*LT*0*0)  A-EN2+PI 
IF(Y<6>*NB*0*0*AND*Y(4>*NE.0.0)  00  TO  666 
GO  TO  670 

666  ANS=ATAN21Y<5)»Y<4>» 

IF!ANS*LT*0*0>  ANS«AHS+2**PI 
A-EN2-ANS 

670  FKX“FKA*ABS ( COS ( EN2 ) )+FKC*ABS<SIN!EN2) ) 

FK Y=FK A*ABS ( S I N  <  EN2 ) ) +FKC*ABS <  COS (EN2I) 

P1«M1+M2*2**FKX 

C  DRAG  «C 1 *COS  <  A ) **2+*  2*ABS ! S I N I A ) **3  > 

CLIFTa*2*(SIN( A)**2)*ABS!C0S(AI >*SION( 1* tSIN(2* *A>> 

V2=Y <4 S  **2+Y<5>**2 

F ( 4  >  =- V2*  <  C  L I  FT*S  IH(EH2-A)  +C  DRAG*COS  ( EN2-A I  >  /PI 
P2“M1 +M2*2 • *FKY 

F ( 5 )  - ( 2 • *M2“M1+V2* <  CLI FT*COS< EN2-A ) »CDRAG*SIN! EN2-A ) » ) /P2 
P3-V2* < CDRAG*SIN <  A » +CLIFT4C0S! A » > 


A  DD=2 • +H2+S I N ! EN2  > ♦V2* ! C  LI FT*S I N  ( A  >  - CDRAG*COS ( A ) > 
DAMP«M3*Y<6)*ABS!YC6> J 

F!6)«<P3*<D1-D3>+2**M2*!D1-D2>*COS!EN2I+ADD*FVCO/XO>/I1-DAMP 

740  GO  TO  660 
750  I F ( XX-X9 ) 780* 780 *770 

770  GO  TO  815 
780  C=C+1 

790  IF(C-0)560.800»S60 

800  C=0 

810  GO  TO  510 

2000  TYPE  2001 

2001  FORMAT!'  THE  VEHICLE  IS  DYNAMICALLY  UNSTABLE* • • INTEGRATION  HAS' 
*/.'  BEEN  TERMINATED*  RESULTS  MAY  BE  PRINTED  OR  GRAPHED* ' .//I 

815  RETURN 

816  END 


820  SUBROUTINE  RUNGEIHtY.F.X.H.lCJ 
822  DIMENSION  P<50) tS«50) »Y(6I »F(6) 

COHMON/RUNQ/M 
828  M=M+1 

H5=H**5 

830  CO  TO  (S80»900»9SO«1040*U10>.N 

880  K»1 

890  RETURN 

900  DO  930  J=1.H 

910  S(J)*Y<J> 

920  PU)»F(J) 

930  Y<J)=S(J)+H5*F<J) 

950  X-X+H5 
960  K*1 

970  RETURN 
980  DO  1000  J«1 »N 
990  P< J>*?( Jt+2**F( J) 

1000  Y<J>»S<JKH5*F<J» 

1020  K-l 
1030  RETURN 
1040  DO  1060  J*1.N 
1050  P<J)*P(J)+2**FU» 

1060  YIJI«S(J)+H*FIJI 

1080  X«X+H5 

1090  K=1 

1100  RETURN 

1110  DO  1120  J-l.N 

1120  Y(J)-S(J>-HP(J)+F<JJ)*H/6*0 

1140  H=0 

1150  K«0 

1160  RETURN 

1170  END 


IX-7. 


Subroutine  PLftT 


This  package  is  called  in  response  to  the  monitor 
command  GRAPH) .  It  furnishes  the  user  with  annotated 
plot  of  the  vehicle's  non-dimensional  vertical  and  hori¬ 
zontal  distances  from  the  point  of  separation  (time  =  0) 
and  of  its  attitude  in  degrees  versus  the  independent 
variable,  non-dimens ionalized  time.  The  vehicle  may  be 
either  the  capsule  or  the  remaining  vehicle  depending 
upon  whether  ASCEND)  or  HULL) ,  respectively  was  the  more 
recent  ccmmancT. 

Determination  of  whether  the  vehicle  motions  being 
plotted  are  for  the  capsule  or  remaining  hull  is  made  by 
checking  whether  the  value  of  the  subroutine  argument  is 
1,  or  2,  respectively.  This  us  used  to  set  the  proper 
label  and  to  compute  vehicle  length  which  are  printed  as 
part  of  the  plot  information. 

The  data  arrays  which  are  plotted  are  passed  through 
the  COMMON  labelled  IN0UT.  Scaling  of  the  arrays  is  per¬ 
formed  by  Typagraph  supplied  software  as  is  the  drawing 
and  labeling  of  the  axes. 

The  subroutine  also  prints  a  table  of  the  vehicle's 
general  characteristics  for  identification  purposes.  The 
completion  of  the  subroutine  is  followed  by  a  return  to 
the  main  program  for  additional  monitor  commands. 

There  cure  no  error  messages. 
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Input  Variables 


VTST  Capsule  -  remaining  vehicle  switch 

XLEN  Intact  vehicle  length  (ft. ) 

SEP  Separation  point  from  nose  (decimal  fraction) 

EL0D  Length  to  diameter  ratio 

TRIM  Initial  trim  angle  (deg.) 

XVEL  Initial  horizontal  velocity  (ft. /sec.) 

YVEL  Initial  vertical  velocity  (ft. /sec.) 

T  Time  array,  non-dimensional 

D  Vehicle  attitude  array  (deg.) 

V  Vertical  distance  array,  non-dimensional 

H  Horizontal  distance  array,  non-dimensional 


Output  Variables 


DAT 

XLEN 

SEP 

EL0D 

TRIM 

XVEL 

YVEL 

T 

D 

V 

H 


Date  of  simulation 
Intact  vehicle  length  (ft.) 

Separation  point  from  nose  (decimal  fraction) 
Length  to  diameter  ratio 
Initial  trim  angle  (deg.) 

Initial  horizontal  velocity  (ft. /sec.) 

Initial  vertical  velocity  (ft. /sec.) 

Time  array,  non-dimensional 
Vehicle  attitude  array  (deg.) 

Vertical  distance  array,  non-dimensional 
Horizontal  distance  array,  non-dimensional 
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IX- 7.1  Flow  Chart  for  Subroutine  PL0T 
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SUBROUTINE  PLOT) ITST.NPTS! 

COMMON/INPUT/XLESi, SEP. ELOD. TRIM. XVEl.YVEL.DVKL 
COMMON’/INOUT/SKIP (2? )  .T<101  >  .DUOl  > » V  <  101 )  .H<  101>.EXt3Q3> 
INTEGER  TTITLE.DTITLE.HTITLE 

DIMENSION  TT IT LEI  121 *DTITLE( 12 ) .HTITLEI 12) > WORK <202) 
DIMENSION  DATI2) *VEH<4) 

DATA<T'i'ITLE(K)  »K=1 .12 1 /'NON -DIMENSIONAL  TIME*. 

*  •  •/ 

DATA (  DTITLEI K ) »K“1 • 12 )/' VEHICLE  ATTITUDE  IN  DEGREES  •» 
*  •  • •/ 

DATA  <  HT  IT  LE  <  K ) » K=1 *  12 ) / ' D I  STANCE  PROM  SEPARATION  POINT 
*  '  •/ 
i  F ( ITST-1 ) 1 > 1  «2 

1  VEh  <  1 )  =  *  CAPSU* 

VEH ( 2 )  =  *  LE  :  * 

VSH<3)=* 

V  EH ( 4 1  =  *  * 

VLEN=XLEN*SE? 

GOTO  3 

2  VEH<  1 )  =  *  REMAI ‘ 

VEH ( 2 )  =  *  NING  * 

VEH(3)  =  *  VEHIC' 

VEH(4>  =  *  LE  J  • 

VLEN=XLEN*(1*“SEPI 

3  CONTINUE 

C  “N"=  LUMBER  GF  POINTS  TO  BE  PLOTTED 
N=  NPTS 
DO  20  1=1. N 
WORK ( I ) =H  <  I ) 

20  WORK ( I +N ) =V ( 1 1 

N2=2*N 

Th=TR I M*57 • 2958 

C  ASSIGNMENT  OF  ARRAY  SYMBOLS  *.  a.  * 

VSYM=*  *• 

HSY.M='  +* 

DSYM=* • * 

C  SCALE  ALL  ARRAYS 

CALL  TGSCALEIT.N.6- .TMIN. TSCALE.l ) 

CALL  TGSCALEt WORK. N2, 7*75, VMIN.VSCALE.il 
CALL  TGSCALEID.N.6* .DMIN.DSCALE.l ) 

C  SET  INITIAL  ORIGIN  FOR  PLOTTING  AXES 
CALL  TGPL0TI9-0.0I 
C  DRAW  LEFTMOST  VERTICAL  AXIS 

CALL  TGAXI S < 0»0.Q.0t 7*75. VMIN.VSCALE. 1.10*2* HTITLEI 
CALL  TGAXIS(0.0.0.0.6‘O.TMIN.TSCALE.4.-10.2.TTITLE) 

CALL  TGAXIS(6.0.0-0.6.0.DMIN.DSCALE. 2. -10.2. DTITLE) 
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C  PLOT  HOK Z Z  DISTANCE 

CALL  TGCUKVE(T> V.N.TMIN. VMIN.TSCALE. VSCALE. VSYM. 1  ) 

C  PLOT  VEriT*  DISTANCE 

CALL  TCCURVEtT .H.N.TMIN.VMIN.TSCALE.  VSCALE.HSYM.il 
C  PLOT  ATTITUDE 

CALL  TCCURVEtT. D.N.TMIN.DHIN.TSCALE.DSCALE.DSYM.l) 

CALL  TO END 
TYPE  30 

30  FORMAT! *  +  * 10?  ??????????????????????????????????? HHQ  * ) 

CALL  DATE! DAT) 

TYPE  40.DAT.YEH.VLEN.EL0D.SEP.THiXYEL.YVEL 
40  FORMAT <1 HO. Tl?. 'SUBMARINE  ESCAPE  CAPSULE  SIMULATION* 

*  T23.'0NK  .  * 2A5./T17. * CHAKACTEHISTICS  OP  ‘4A5./. 

*  T2S. ’LENGTH  (PT)  :  * .F8«2./.T15» ‘LENGTH  TO  DIAMETER  RATIO  I  *. 

*  F3.2./.T14. ‘SEPARATION  PLANE  LOCATION  :  *.fa*2./. 

*  Tlo. ‘INITIAL  TRIM  ANGLE  (DEC)  :  *.P 8*2./. 

*  T12.’ INITIAL  X  VELOCITY  (FT/SEC)  :  ’.F8.2/. 

*  T12.* INITIAL  Y  VELOCITY  (FT/SEC)  :  *.F8.2./, 

*  T3o.*X/L  ;  +' ./.T36. * Y/L  :  ** ./.T34, 'THETA  t  •*» 

TYPE  60 

50  FORMAT! '+* JO?  ??????????????????????????????????????* ) 

CALL  TGEND 

RETURN 

END 


IX-8. 


Subroutine  PRNT 


The  data  arrays  describing  the  vehicle's  (capsule 
or  remaining  vehicle)  motions  following  separation  may 
be  plotted  (cf.  Subroutine  PL0T)  or  printed  out  in 
tabular  form  (Subroutine  PRNT) .  PRNT  allows  the  user  to 
sample  two  to  fourteen  points  of  each  array  (the  first 
and  last  points  of  the  array  are  alvrays  included)  for 
inspection  of  their  numerical  values.  The  user  calls 
the  subroutine  from  the  monitor  by  typing  PRINT) .  The 
subroutine  responds  by  requesting  the  number  of  points 
the  user  wants  to  print  and  uses  that  number  to  calculate 
an  approximately  equally  spaced  sample  for  tabulation. 
Because  of  integer  truncation  and  other  physical  limita¬ 
tions,  the  algorithm  that  selects  the  points  to  be 
plotted  begins  to  lose  its  accuracy  if  more  than  four¬ 
teen  points  are  requested.  No  restrictions,  however, 
are  placed  on  any  number  requested  between  2  and  101. 

The  error  message  "ERROR,  NUMBER  OF  POINTS  MUST  FALL 
WITHIN  THE  RANGE  [2,  101].  TRY  AGAIN."  is  given  if 
the  user  tries  to  print  less  than  2  or  more  than  101 
points.  The  user  is  then  again  asked  "HOW  MANY  POINTS 
DO  YOU  WISH  TO  PRINT  2  TO  14  — " . 

Non-dimensional  results  are  first  printed  out  as  a 
table.  Then  the  values  are  dimensionalized  and  printed 
out  as  a  second  table. 
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Input  Variables 


FI  Length  of  vehicle  (ft.) 

T  Time  array,  non-dimensional 

D  Vehicle  attitude  array  (deg.) 

V  Vertical  distance  array,  non-dimensional 

H  Horizontal  distance  array,  non-dimensional 

DD0T  Angular  velocity  array  (deg. /non-dim.  time) 

VD0T  Vertical  velocity  array,  non-dimensional 

KD0T  Horizontal  velocity  array,  non-dimensional 

N  Number  of  points  to  be  printed 

NPTS  Number  of  points  in  each  array 


Output  Variables 

T  Time  array,  non-dimensional 

D  Vehicle  attitude  array  (deg.) 

V  Vertical  distance  array,  non-dimensional 

K  Horizontal  distance  array,  non-dimensional 

DD0T  Angular  velocity  array  (deg. /sec.) 

VD0T  Vertical  velocity  array,  non-dimensional) 

KD0T  Horizontal  velocity  array,  non-dimensional 

ST  Time  array  (sec.) 

SV  Vertical  distance  array  (ft.) 

SH  Horizontal  distance  array  (ft.) 

SVD  Vertical  velocity  (ft. /sec.) 

SHD  Horizontal  velocity  (ft. /sec.) 

SDD  Angular  velocity  (deg. /sec.) 


I 


IX- 8. 1  Flow  Chart  xor  Subroutine  PRNT 


( 
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IX-8.2  Listing  of  Subroutine  PRNT 


( 


( 
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O  Ci 


surkout  i  :;s  psst «  fl.spts) 

DIMENSION  STdGll  .SHflOl)  .SHDdOU.SVdOl  >, 

*SVi)d01 1  tSDDdOl ) 

COMMON/INPUT/FLOA. SEP. FLDRAT. TRIM. XVELtYVEL.DVEL 
common/ iN0UT/$:<i?t2?>  .Tdou.DdoidVdon.Hdoi). 

♦  DDOT UOl ) .  VDOT  (101 ) .HOOT  < 101 ) 

i  type  5 

ACCEPT  6 iW 

I F l N • LT*  2 • OR • K *01  •  NPTS ) GO  TO  100 

TYPE  10 
TYPE  20 

10  FORMAT  ( /  //  •  1H-  .T24  *  *  NON-DIMENSIONAL  RESULTS') 

20  FORMA? (IK  . T4.' TIME', T16, ' HORI ZONTAL* .  T37, ' VERTICAL'. 

*To2 » 'ATTITUDE'  *T62»  'ANGULAR*  »/»T12»  'DISTANCE*  »T22, 

*'  VELOCITY' »T32» 'DISTANCE' *T42» ' VELOCITY* .T52# 'DEOBEES' » 

*T62, 'VELOCITY' ) 

TYPE  30 

30  FORMAT  dH  ,T3.' - '.T12.' - '.T22.' - '  ,T32, 

. . . . . '  ,T52.  * - * •  T62 • ' - '  ) 

3  FORMAT! *  HOW  MANY  POINTS  DO  YOU  WISH  TO  PRINT  12  TO  14)  “ 
*s ) 

6  FORMAT! I) 

KS?E?=( NPTS-l  )/CK-l ) 

DO  30  J=1 » tipTS-KSTE?  *  KSTEP 

50  TYPE  '10«T(  J )  .ri(  J )  .HDOTI J  ) » V  ( J ) »  VD0T(  J)  »D(  J )  *DDOT<  J) 

T  Y P  E  40 » T ( NPT  S ) . K ( NPTS ) . HDOT ( NPTS ) .  V ( N  PTS) . VDOT  <  NPTS  > . 

*D( NPTS I »DDO? ( NPTS ) 

40  FOr.MAT I1H  ,?(F3*2»2X)) 

CONVERSION  TO  DIXENSIONALIZED  VALUES  IS  COMPUTED  NEXT* 


F DO-SORT ( FL/32* 17 ) 
FTO=SQRT ( F L*32* 17  > 
C  DF=SGi<? ( 32 • 1 7/FL) 
DO  GO  1=1, NPTS 
ST ( I ) =?  < I ) *FDG 
SH( I )=K< I )*FL 
SHD (I ) =HDOT ( I ) *FTG 
SV(I)=Vd)*FL 
SDD(  I  )*ODO?d  l*GDF 
60  SVD( I )=VDOT 1 1 )*F7G 
C 
C 

c 
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TYPE  70 

70  FORMAT!  1H-,'T24» 1 D I  MENS IONA LI  ZED  RESULTS*) 

TYPE  20 
TYPE  SO 

SO  FORMAT (1H  ,T3. * SECOUUS* ,T14. 'FEET* ,T23. 'FT/SEC* ,T34, ’FEET* *T43» 

*' FT/SEC*  »‘i'63,  *  DEG/SEC*  ) 

TYPE  30 

DO  SO  J=1.K?TS-KSTE?.KSTEP 

00  TYPE  40»ST( J ) » SK( J) »SHD( J ) »SV( J 1 ,SVD( J) »D( J) «  SDD(  J ) 

TYPE  40,ST(i!PTS)  .SK(NPTS) ,SHD!N?TS) ,SV(NPTS) .SVD(NPTS) , 

*DC:?T5) *SDD(K?TSI 
GO  TO  120 

100  TYPE  llOiMPTS 

110  FORMAT (*  ERROR,  NUMBER  OF  POINTS  MUST  FALL  WITHIN*,/ 

*  '  THE  RANGE  C2.*.I3.’]*  TRY  AGAIN**,/) 

CO  TO  1 

120  TYPE  130 

130  FORMAT!//) 

RETURN 

END 
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IX-9. 


Subroutine  PDATA 


In  order  to  provide  the  user  with  a  current  status 
of  what  he  has  done,  the  command  CHECK)  has  been  in¬ 
corporated  into  the  monitor.  This  command  calls  sub¬ 
routine  PDATA  which  types  out  the  status  of  the  simu¬ 
lation.  It  does  this  by  testing  the  status  variables 
of  the  pertinent  subroutines  which  have  been  passed 
as  subroutine  arguments.  If  the  test  is  affirmative, 
a  message  is  typed  out  notifying  the  user  that  the  sub¬ 
routine  in  question  has  been  exercised.  The  general 
characteristics  of  the  intact  vehicle  are  next  listed 
and  the  user  has  the  option  of  inspecting  the  current 
weight  and/or  volume  distribution. 

Upon  completion  of  the  status  report,  control  re¬ 
turns  to  the  monitor. 

There  are  no  error  messages. 
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Input  Variables 

ITST  Subroutine  status  check  array 

FL00D  Floodable  spaces  distribution  array  (ft.s, 

dec.  frac. ,  dec.  frac.) 

WEIGHT  Weight  distribution  array  (lbs.,  dec.  frac., 

dec.  frac.) 

NUMWGT  Number  of  weights  in  WEIGHT  array 

NUMFLD  Number  of  spaces  in  FL00D  array 

CAPT  Capsule  shell  thickness  (in.) 

HULLT  Hull  shell  thickness  (in.) 

FLEN  Intact  vehicle  length  (ft.) 

SEP  Separation  point  from  nose  (decimal  fraction) 

FLDRAT  Length  to  dismeter  ratio 

TRIM  Initial  trim  angle  (radians) 

VELX  Initial  horizontal  velocity  (ft. /sec.) 

VELY  Initial  vertical  velocity  (ft. /sec.) 

DINIT  Initial  depth  (ft.) 


Output  Variables 


FLDRAT  Length  to  diameter  ratio 

FLEN  Intact  vehicle  length  (ft. ) 

SEP  Separation  point  from  nose  (dec.  frac.) 

CAPT  Capsule  shell  thickness  (in.) 

HULLT  Hull  shell  thickness  (in.) 

VELX  Initial  horizontal  velocity  (ft. /sec.) 

VELY  Initial  vertical  velocity  (ft. /sec.) 

DINIT  Initial  depth  (ft.) 

TR  Initial  trim  angle  (deg.) 

WEIGHT  Weight  distribution  array  (lbs.,  dec.  frac., 

dec.  frac.) 

FL00D  Floodable  spaces  distribution  array  (lbs., 

dec.  frac.,  dec.  frac.) 
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IX-9.1  Flow  Chart  for  Subroutine  PDATA 


( 
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IX-9.2  Listing  of  Subroutine  PDATA 


SUSKOUT  I NE  PDATA ( ITST, FLOOD,  WE  IGiiT » N  UM te’vJT « KUMFLJ) » CAPT» HULL?* 
trLOv}I.'<t  \  CG I V  V  t  FLOG OP  t  V  GG  CA  P ) 

CONNON  , SEP*  FLORA?* TRIM*  VELX*  Vi!LY * VELD 

dimension  i?Si  <  30  > » flood <  3. 100 » .  viei ghti 3. 100 » 

T'iVE  10 


FO;. ::/.;(// 1 lox. • present  simulation  STATUS’*/) 

I  F  ( IT  ST  (291*  EU  •  i  )  T  YPE  202 
I  F  ( IT S T 1 1  )  •  EG  •  i  )  L":  PE  201 
IF ( ITSTIV;  .EO.DTYPE  207 
IFi  US  1(12)  •FC-ilTYPS  212 
I F  (  I  ■.  si ( 8  > • EG • 1 )  TY  PL  203 
I F !  I TST  ( 9 )  •  EC • i  )  TY 22  209 
I F ( M  ST ( 1 0 ) • EG • 1 ) TYPE  210 
.  r.= V  A  I X+3  7  •  2  908 

TYPE  110 .  FLDAA? .  VLE!< .  SE?  *CAP?«  HULL?*  VELX* VELY  »OIMIT*TH 
t  r  L^G  1 1. » v  CG  IT  v  t  F  LCGCP  t  VCGCAP 
T:PF  Ill 
A  C  C  HPT  oO  1 1 S I G 

ifi i sic igo  to  is 

lo  ii  i=i,:;um.,gt 

V/SIGriT  1 1 1 1 )  ='.:SIGHT  11,1) /2240* 

TYPE  112.(1  ii'alGH? 1 1 ,  J ) ,  1=1 .3) , J=1,NUMWG?) 
lO  12  1=1 1  NuHKGT 


W2I GST  ( 1 » I )  =V;  FIG  H  i  (1,1)  *2240* 

TYPE  113 

.".LUr.Pi  oOtiLit. 

I F  ( I S IG  t  -:c  • '  Y  ’ )  T  i  PF  114 ,( ( FLOOD ( I ,  J ) ,  1=1 , 3)  ,  J=1  iNUMFLl) ) 

RETURN 


FOAM AT ( / * 


LENGTH/DIAMETER  HAT  10 
TOTAL  VEHICLE  LENGTH (FT) 
SEPARATION  PLANE  LOCATION 
CAPSULE  SHELL  THICKNESS!  IN’. ) 
Hlul  SHELL  THICKNESS! IN. ) 
INITIAL  HORIZONTAL  VELOCITY 


'  INITIAL  VERTICAL  VELOCITY 
'  INITIAL  DEPTH (FT) 

’  INITIAL  "HIM  ANGLE (DEGREES) 
’  LCG  INTACT  VEHICLE 
'  VCG  INTACT  VEHICLE 
'  LCG  CAPSULE 
'  VCG  CAPS JLt 


=’  F10*2,/ 
=’  F10*2 , / 
=’  F10.2,/ 
='  F10.2./ 
=’  F10.2,/ 
=*  F10«2./ 
='  F10«2,/ 
=•  F10.2,/ 
='  F10#2,/ 
=  *  F10*3/ 

=  ’Ft  0-3/ 

=  *  F10*3/ 
=  *  F10»3/ 


) 


ill  EUR MAT  I  /  1  DO  Y  OU  Vi  I  iH  TO  INSPECT  TjiE  CUKHKNT  HEIGHT1 

*  •  DISTRIBUTION1/1  TYPE  y  OH  -  ■  S ) 

US  :0 iiMAri/'  DO  YOU  WISH  TO  It! SPEC?  THE  CURRENT  VOLUME' 

*  •  BIS .’A  I  BUY 1 0»«? ' / '  TYPE  Y  OH  K  -  •$) 

IIP  *O..MA?</'  CURRENT  ..'EIGHT  DISTRIBUTION  IS:'//.(oX.SF12*2)) 

114  POL MAT  I / '  CURRENT  VOLUME  DISTRIBUTION  IS  i '//« (5X»3F12*2)  ) 

Ou  -  0  ...'Av  (  a  1  ) 

202  PUR MAT  </'  CURRENT  WEIGH?  DI STRIBUTION  OH  VOLUME  DIS* 

*  'TRIBUvION1/  •  HAS  BEEN  UPDATED  *  J 

201  PO..MAT C / '  GEOMETRICAL  ? A HA ME? EH S  HAVE  BEEN  CALCULATED') 

20?  PORM.VYI/'  A  CASUALTY  HAS  BEEN  SPECIFIED* ) 

212  PO.-.MA.  (/'  WEIGHTS  AND  CENTERS  HAVE  BEEN  CALCULATED1) 

2Ca  PO.-.MAV  (/*  INTACT  VEHICLE  ATTITUDE  HAS  BEEN  CALCULATED1) 

209  PJ.'.MAY  (  / 1  CAPSULE  ACCENT  HAS  BEEN  CALCULATED1) 

210  POEM AT ( / 1  REMAINING  VEHICLES  PATH  HAS  BEEN  CALCULATED* ) 

END 
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IX- 10.  Subroutine  ADD0N 

This  subroutine  updates  the  weight  and  volume  dis¬ 
tribution  files  maintained  on  disk.  ADDON  is  not  called 
directly  by  a  command;  however,  the  option  to  add  or  up¬ 
date  distributions  is  provided  when  a  distribution  is 
input  from  the  keyboard  or  when  a  distribution  has  been 
updated.  Distributions  may  be  added  to  the  file  or  in 
the  case  of  an  updated  distribution,  the  old  distribution 
on  file  is  replaced  by  the  updated  version.  The  user 
must  input  a  "distribution  identifier"  which  is  5  characters 
long;  it  may  be  alphanumeric  and  must  be  unique  in  the 
file.  If  the  user  inputs  an  identifier  for  which  there 
already  exsists  distribution,  he  may  not  input  a  new  identifier 
for  his  distribution.  Once  this  operation  is  completed 
the  message  "DISTRIBUTION  HAS  BEEN  STORED  FOR  FUTURE  USE" 
is  output  and  control  returns  to  the  monitor. 


There  are  no  error  messages. 


Inout  Variables 


ARRAY 
ITYP 
I  UN  IT 


Weight  or  volume  array  to  be  output  on  file 
Specifies  that  array  contains  weights  or  volumes 
Specifies  weight  or  volume  file  number 
Number  of  triplets 


Output  Variables 


None 


IX-10.1  Flow  Chart  for  Subroutine  ADD0N 


( 
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T>TSr  A fctJ  ote  uP’OftTs) 


IX-10.2  Listing  of  Subroutine  ADD0N 


( 


( 
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TYPE  ADDON'. EOS 

SUBROUTINE  ADDON (ARRAY, ITYP, IUNIT.N) 

DIMENSION  ARRAY (3,100), SCR<3) 

C*******ROUTINE  TO  ADD  DISTRIBUTION  TO  DATA  FILE* 
V0LM=5hT0LUM 
Vt’ATE=4HWATE 
NSV/=0 
ISIG=*  • 

CALL  0FILEI22, 'SCRAT' ) 

IFUTY?*EQ*1  )G0  TO  10 

C*+**+.i(<niYP=l  ADD  TO  'JOLUMESt  ITYPE=0  ADD  TO  WEIGHTS* 

20  CALL  IiTLE(21,WATE> 

CO  TO  15 

10  CALL  If ILE(23»V0LM) 

C****,**FIL2s  OPEN  20=WEIGHT  FILE,  21=VV0LUMS  PILE,  ^SCRATCH 
izri,  50 

ACCEPT  51  *  I DENT 
RSADt I UNIT, 52 JMTDENT.NO 
IF(EOFCdUNIT) )  GO  TO  36 
IF(MTDEST*EQ*IDENT)GO  TO  26 
Wit  I TE  ( 22  •  52 )  MTDE1IT  •  NO 
DO  30  1=1, NO 

R2AD< IUKIT.53) (SCR(J) ,J=l,3) 

WSITE(22.53)(SCR(J),J=i,3) 

GO  TO  25 
END  FILE  22 
GO  TO  101 
TYPE  54 

FORMAT ( '  is  DISTRIBUTION  NEW  OR  AN  UPDATE  ?V 
'  TYPE  N  OR  U  —  •*) 

ACCEPT  55.  ISIC 
FORMAT (All 
I?( ISIG*EG* 'N* )G0  TO  28 
C*******'//  rite  UPDDA  TEPD I STR I  BUT  I  ON 
WRITE<22,52)  IDENT.N 
DO  29  1=1, N 

29  WRITE (22 <53 > (ARRAY ( J, I ) , J=l ,3) 

C********SKIP  OLDD  DISTRIBUTION. 

DO  40  1*1, NO 

READ( IUKIT.53) ( SCR( J I , J=1 ,3) 

GO  TO  25 

NSW=1 

Vi?.  I T  E  ( 22 , 52 )  M  T  DENT  ,  NO 
DO  32  1=1, NO 

READ(IUNIT»53) (SCR ( J) , J=i ,3) 

WRITE (22, 53) (SCR( J),J=l,3) 

READ ( I UNIT » 52 ) KTDENT » NO 
■IF  (EOFCUUNIT)  )G0  TO  100 
GO  TO  31 


30 


27 


26 

54 


55 


40 


28 

31 


33 


32 
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100  IF(NSW*EQ*1)G0  TO  34 

101  END  FILE  22 

IF( ITY?*EQ*1 ICO  TO  39 
CALL  0FILE(21.WATS> 

CO  TO  103 

39  CALL  0FILE(23» VQLH) 

103  CALL  IFILE(22. ' SCRAT' I 

C*******START  COPY  FROM  SCRATCH  FILE  TO  ROIC  INAL  IUNIT* 

105  READ(22,52)MTDENT.NQ 

IFt E0FC(22) ICO  TO  102 
WRITE  ( IUNIT»52)MTDENT»N0 
DO  104  1=1. NO 

READ (22. 53 1 ( SCR ( J ) »J=1.3) 

104  WRITE! 1 UNIT .53 1 (SCR(J) »J=1.3I 

GO  TO  105 

102  END  FILE  IUNIT 
TYPE  5? 

5?  FORMAT!  *  DISTRIBUTION  HAS  BEEN  STORED  FOR  FUTURE  USEV) 
RETURN 

C*******WRITE  NEW  DISTRIBUTION  AT  END  OF  LOG*  4* 

34  TYPE  56 

56  FORMAT! *  PLEASE  INPUT  IDENTIFIER  OF  NEW  DISTRIBUTION*/ 

*  '  5  CHARACTERS  ALLOWED  ~  **) 

ACCEPT  51. I DENT 

36  IF(ISIG*EQ.*U')GO  TO  101 

WRITE  (22.52 JIDENT.N 
DO  35  1=1. N 

35  WP.ITEI22.53) (ARRAY ( J. 1 1 • J=l»3) 

N  SW=0 
GO  TO  101 

50  FORMAT!*  ENTER  SHIP  "IDENTIFIER".  (5  CHARACTERS)  —  **) 

51  FORMAT (A5 I 

52  F0RHAKA5. 13) 

53  FORMAT (3F I 
END 


IX- 11  Subroutine  UPDATE 

This  routine  is  called  in  response  to  the  command 
UPDATE).  The  subroutine  provides  the  user  the  option  of 
updating  either  the  weight  or  volume  distribution  once  it 
is  in  storage.  The  user  may  at  his  option  change,  delete, 
or  add  a  weight  or  volume  and  its  centers.  Depending  upon 
the  operation  the  user  must  input  the  old  triplet  that 
is  to  be  deleted  or  changed  and  he  must  also  input  the  new 
triplet  if  it  is  to  be  changed  or  added.  The  user  must 
take  caution  to  remember  that  the  weights  must  be  input 
in  pounds  for  this  routine.  Once  the  distribution  is  up¬ 
dated,  the  user  is  given  the  option  of  storing  the  dis¬ 
tribution  for  future  use.  Program  control  then  returns 
to  the  monitor. 

Error  Messages 

1)  "WEIGHT  OR  VOLUME  AND  POSITION  CANNOT  BE  FOUND”  - 
The  user  has  attempted  to  change  or  delete  a  weight  or 
volume  which  does  not  exist  in  the  distribution  or  which 
he  may  have  input  incorrectly.  (Note:  Weights  must  be 
expressed  in  pounds.) 

2)  "DETECT  ERROR  IN  INPUT  TRY  AGAIN"  -  The  command 
input  by  the  user  for  changing,  adding,  or  deleting  is 
incorrect. 


Input  Variables 

ARRAY  Weight  or  volume  array  to  be  updated 

ITYP  Specifies  whether  weights  or  volumes  are  to 

oe  updated. 

N  Number  of  weight  or  volume  triplets  in  the 

distribution 


Output  Variables 
ARRAY  Update  array 
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IX-11.2  Listing  of  Subroutine  UPDATE 

fr 

( 


( 
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SUBROUTINE  UPDATE! ARRAY* ITYP.N) 

DIMENSION  ARRAY (3.100) 

1  TYPE  50 

50  FORMAT! '  IS  WEIGHT  TO  BE  UPDATED  TO  BE  CHANGED*  DELETED. * 

*  •  OR  ADDED?*/'  TYPE  C.D.  OR  A  —  '.!) 

42  ACCEPT  51.ISIG 

51  FORMAT (All 

IF! ISIG*£Q* 'A* IGO  TO  30 
IF! ITYP*EQ*1 >G0  TO  10 

C*******ITYP=1  UPDATE  VOLUMES.  ITYPE-0  UPDATE  WEIGHTS* 

TYPE  54 

52  FORMAT!*  INPUT  PRESENT  VOLUME.  LCG.  AND  VCG  —  '.»> 

ACCEPT  53.  VO.FLC.VC 

53  FORMAT <3F> 

GO  TO  20 

10  TYPE  52 

54  FORMAT! *  INPUT  PRESENT  WEIGHT.  LCG.  AND  VCG  —  *.*) 

ACCEPT  53. VO.FLC.VC 

GO  TO  20 

C*******  ADD  NEW  WEIGHT  OR  VOLUME  TO  DISTRIBUTION* 

30  IF! ITYP*EQ*1 )G0  TO  11 
TYPE  55 

65  FORMAT! *  INPUT  WEIGHT.  LCG.  VCG  TO  BE  ADDED  —  '.*) 

ACCEPT  53.  VO.FLC.VC 
,  V0»V0*2240* 

GO  TO  15 

11  TYPE  56 

56.  FORMAT  <•  INPUT  VOLUME  .  LCG.  AND  VCG  TO  BE  ADDED  —  '.*) 
ACCEPT  53. VO.FLC.VC 
15  N=N+1 

ARRAY  <1.N)=.V0 
ARRAY(2.N)-FLC 
ARRAY <3«N)=VC 
GO  TO  100 

C  *«*4ut>*BEGIN  SEARCH  FOR  CHANGE  OR  DELETE  OPTION* 

20  DO  35  1=1. N 

IF! I FIX (ARRAY <1.1 1*100*  >  *NE* IFIX! V0*10Q*  ) )G0  T035 
IF! I FIX (ARRAY <2*1 )*100* l *NE*IFIX!FLC*100* ) IGO  TO  36 
IF!IFIX<ARRAY<3.I)*100*>*NE.IFIX!VC*i0O»))G0  TO  35 
MATCH=I 
GO  TO  36 

35  CONTINUE 
TYPE  57 

57  FORMAT!'  WEIGHT  OK  VOLUME  AND  POSITION  CANNOT  BE  FOUND') 
GO  TO  100 

36  IF!ISIG*EO*'C' )G0  TO  40 
IF<ISIG*EQ*'D')GO  TO  38 
TYPE  58 

58  FORMAT ('  DETECT  ERROR  IN  INPUT.' TRY  AGAIN') 

GO  TO  1 
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38  ¥0-0*0 

FLC-0*0 
YC-0*0 
GO  TO  3? 

40  TYPE  59 

59  FORMAT C  INPUT  NEW  WGT  OR  VOLUME*  AND  LCG.  VCG  —  *.*> 

ACCEPT  53* VO*FLC* VC 

I F t ITYP*£Q*0)  V0-V0*2240* 

37  ARRATI1. MATCH) -VO 

ARRAY  12. MATCHJ-FLC 
ARRAY )3.MATCH)»VC 
IF) ISIG*NE* 'D' >GQ  TO  100 
NEND-N-MATCH 
N=N-1 

DO  25  I-O.NEND 

AR  RA Y (1 . MATCH-*- 1 ) -ARRA Y ) 1 . MATCH+I +1 ) 

ARRAY (2.MATCH+I )-ARRAY (2.MATCH+I+1 ) 

ARRAY < 3.MATCH+I )-ARRA Y ( 3.MATCH+I+1 ) 

25  CONTINUE 
100  TYPE  60 

60  FORMAT) *  DISTRIBUTION  UPDATED.  AMY  MORE  CORRECTIONS?*/ 

*  '  TYPE  Y  OR  N  —  *.*> 

'  ACCEPT  51.ISIG 

IF) ISIG*EQ*  *Y* )GOTO  1 
TYPE  62 

62  FORMAT)*  DO  YOU  WISH  TO  STORE  THIS  DISTRIBUTION  FOR  FUTURE  USE* 

*  /*  TYPE  Y  OR  N  —  •*) 

ACCEPT  51.  ISIG 
IFJISIG'EO^'N* )G0  TO  43 
IF) ISIG»NE*  *Y* »G0  TO  43 

IF) ITYP»E0*1 1CALL  ADDON (ARRAY. I TYP* 23. N) 

IF(ITYP’EQ*0>CAU  ADDON (ARRAY. ITYP. 21 .Ml 
43  RETURN 
END 


a 
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IX-12 .  Subroutine  ASTMES 

System  messages  in  response  to  the  monitor  command, 
HELP} ,  are  read  from  a  data  file  and  printed  out  at  the 
user's  terminal  through  this  subroutine.  The  subroutine 
opens  the  previously  created  file,  reads  72  characters, 
checks  the  last  character,  the  next  to  last  character, 
the  second  from  last  character,  etc. ,  until  it  finds  a 
non-blank  character,  then  goes  to  a  printing  loop  and 
prints  one  line  up  to  and  including  the  first  non-blank 
character  mentioned  above. 

The  next  72  characters  are  then  read  and  a  similar 
procedure  is  followed  until  an  end  of  file  mark  is  en¬ 
countered  at  which  time  control  returns  to  the  monitor. 

There  are  no  error  messages. 


Input  Variables 


— None — • 


Output  Variables 

'  MESS  Message  array  containing  one  line  of  text 
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IX- 12.1  Flow  Chart  for  Subroutine  ASTMES 


( 
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IX- 12. 2  Listing  of  Subroutine  ASTMES 
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SUBROUTINE  ASTHES 
DIMENSION  MESS (72 1 
CALL  IFILEd.  'HELP*  I 

1  IF  (EOFl(ll)  GOTO  99 

READ( 1.10). (MESS! I) .1=1 .72) 
10  FORMAT (72A1) 

IFiMESSd ) *EQ» 'S' )  GOTO  00 

J=72 

2  IF(MESS(J>*NE*'  ')  GOTO  20 

J=J-1 

IF(J)  1.1.2 

20  TYPE  30»(HESSU>.I-1.J> 

30  FORMAT (1H  .72A1> 

GOTO  1  * 

60  TYPE  31 

31  FORMAT  < / 1 

,  GO  TO  1 
99  RETURN 

END 


IX- 13. _ Subroutine  RDVW 


RDVW  retrieves  weight  or  volume  distributions  from 
their  respective  files  if  the  commands  INPUT  WEIGHTS  or 
INPUT  VOLUMES  is  used  and  the  option  to~use  old  weights 
or  volumes  is  chosen.  The  user  must  input  the  identifier 
of  the  distribution  he  wishes  to  input  from  the  file. 

When  the  distribution  has  been  retrieved  the  message: 
"DISTRIBUTION  IS  NOW  AVAILABLE  FOR  COMPUTATION  AND  HAS 
BEEN  SAVED  FOR  FUTURE  USE"  is  output,  and  control  returns 
to  the  monitor.  — 

Error  Messages 

1)  "IDENTIFIER  XXXXX  NOT  ON  FILE,  RETURN  TO  MONITOR" 
The  identifier  XXXXX  has  been  input  incorrectly  or  that 
distribution  does  not  exsist  on  file. 

2)  "ERROR  IN  DATA  FILE,  RETURN  TO  MONITOR"  -  The 
number  of  weights  or  volumes  in  the  distribution  is  in¬ 
correct.  The  distribution  should  be  reentered  from  the 
keyboard  and  the  file  updated. 


Input  Variables 

ARRAY  Array  which  retrieved  distribution  is  to  be 

stored  in. 

IDENT  Identifier  of  distribution  to  be  retrieved 

Specifics  whether  weights  or  volumes  are  to 
be  input 


Output  Variables 

ARRAY  Array  containing  distribution 

N  Number  of  triplets  in  distribution 
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IX-13.1  Flow  Chart  for  Subroutine  RDVW  j 

;  I 
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IX- 13. 2'  Listing  of  Subroutine  RDVW 


SUBROUTINE  RDVW( ARRAY *IDENT*N*ITYP) 

DIMENSION  ARHAY (3*1001 
WATE«4HWATE 
V0LM=5HV0LUM 
61  FORMAT (All 

TYPE  53 

53  FORMAT < *  PLEASE  INPUT  DISTRIBUTION  IDENTIFIER'/ 

*  •  5  CHARACTERS  ALLOWED  —'ll 
ACCEPT  52* I DENT 

52  FORMAT «A5I 

IF<ITYP<EQ<OICO  TO  10 
IUNIT=23 

CALL  IFILE(IUNITtVOLM) 

CO  TO  12 
10  IUNIT=21 

14  CALL  IFILE(IUNIT.WATE) 

12  READi IUNIT*54)  MTDENTtN 

IF(EOFC( IUNITt (CO  TO  100 
IF(MTDENT<EQ<IDENT>CO  TO  101 
DO  15  I*1.N 

15  READ( IUNIT*51 (DUMMY 
CO  TO  12 

101  <  DO  102  1=1. N 

'  READ(IUNIT.55)<ARRAY<J.I>.J-1.3) 

102  IF(EOFC(  IUNITMCO  TO  99 
TYPE  56 

55  FORMAT (3F) 

54  FORMAT (A5.I3) 

56  FORMAT ( *  DISTRIBUTION  IS  NOW  AVAILABLE  FOR  COMPUTATION  AND*/ 

*  •  HAS  BEEN  SAVED  FOR  FUTURE  USE'/) 

RETURN 

100  TYPE  57  . 1  DENT 

57  FORMAT <  *  IDENTIFIER  'A5.'  NOT  ON  FILE*  RETURN  TO  MONITOR') 

REWIND  IUKIT 
RETURN 

99  TYPE  58 

56  FORMAT ( '  ERROR  IN  MUCH?  DATA  FILE*  RETURN  TO  MON* 

*  ' ITOR* ) 

RETURN 

END 
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IX-14.  Surface  Motion  Predection  (WVE) 


This  subroutine  calculates  the  heave  and  pitch 
response  of  the  capsule  as  a  function  of  the  sea  state. 
Pierson-Muskowitz  fully  developed,  long-crested  sea  spectra 
are  assumed.  The  capsule  is  assumed  to  float  in  one  of 
three  attitudes;  vertical  bow  up,  vertical  bow  down,  or 
horizontal.  The  attitude  which  yields  the  largest  meta- 
centric  height  is  selected  as  the  most  stable.  Once  the 
most  stable  attitude  is  selected  the  response  is  calcu¬ 
lated.  The  subroutine  prints  out  the  attitude  of  the 
capsule  on  the  surface.  The  subroutine  also  plots  graphs 
of  the  pitch  angle  and  heave  amplitude  as  a  function  of 
wave  length. 


’  20$ 


IX- 14. 1  Flow  Chart  for  Subroutine  WVE 
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SUBROUTINE  N V E ( FLD. VCG » FLEN . BAREA . XLEN. TWOT. FLCC» I SEA » 
DIMENSION  H(10ltNAMH(12>*NAMS<12).NAHP(12).l:AMSP<12). 
*NAKPP(12I 

DIMENSION  S 1100 )»SPR( 100) tN(100) 

DATA ( H ( I ) t I«ltl0)/«12**8S*2*25»4»23»6»78tl0»»lfl*» 
*32.5.&2*5t8S./ 

DATA (NAMW(K) *K"1 »12)/’FREQUENCY*  RAD/SEC  ’tfl*'  V 

DATA(NAMSIK>.K«1.12»/‘HEAVE  AHPL  SQRD/FRSQ* »8*‘  */ 

DATA(NAMP(KJ  »K*1»12)/' (HEAVE  AMPLITUDE ) **2/FBX6  VS*  ?REQ  •# 
•6*'  V 

DATA(NAMSP(K) *K»1*12)/'PITCH  ANCLE  SQRD/FBEQ  **7*'  V 

DATA(NAMPPtK) »K“1#12I/'PITCH  ANCLE  SQ RD/FREQ  '•7*'  \t 

DIMENSION  SUM(IOO) fROV(lOO) 

DIMENSION  MESS(3.3> 

DATA((MESS(I*J)*J-1»3>.1-1*3>/'  BON  ,.,DOWM••.,  *»'  BON  % 

*  *UP.  HORI*  **Z0NTA*» *LLY*  */ 

REAL  ISAVE* ISAVE2»ISAV£3*JSAVE* jSAVE2* JSAVE3 
INTEGER  X  *  Y  •  Z  • RON » SUM *TCELLS* SUMX 
DO  301  1-1.100 
301  SUM( I )**0 

CALL  ERRSET(O) 

PI«3*141S9 

C  COMPUTE  CONSTANT  PORTIONS  OF  EQUATION* 

"  I SAV  E» ( FLEN/2 • > -XLEN 

ISAVE2«>(  FLEN/2*  )**2 
ISAVE3-XLEN/99* 

IF(XLEN/FLEN.LT.*5IY0-S0RT(BAREA/PI» 

IFIXLEN/FLEN.GE**5>Y0-FLEN/(2**FLD) 

JSAVE-YO 

JSAVE2«( ( FLEN/2 )/FLDI**2 
JSAVE3*2*/99« 

SUMX-0 

TCELLS-0 

NC  ELLS- 1 F IX  <  TNGT/ ( 256 • * ( YO* YO ) *XLEN ) *10 • **6 • ) 

C  FIND  REGION  OCCUPIED  BY  CAPSULE (VERTICAL  ATTITUDE) 

DO  5  X-ltlOO 

IX-X 

SX«X-1 

ALOC-( ( (SX*ISAVE3)+ISAYE)**2>/ISAVE2 
DO  10  Z»1.50 
SZ-Z-1 

CLOC-I ( ( <  SZ*JSAVE3>«‘JSAVE>-Y0>**2>/ISAYE2 
SAVLOC»ALOC+CLOC 

DO  20  Y-1.50  \ 

SY-Y-1 

BLOC-( ( ( (SY*JSAVE3»x>JSAVE)-Y0)**2)/USAVE2 
IF(AL0C+BL0C+CL0C*GT*1*)C0  TO  20 
SUM(X>-SUM(X)+4 
ROW(Z)«ROW(Z)+2 
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ROtf(100-Z)-ROV<100-ZI+2 
20  CONTINUE 

10  CONTINUE 

TCELLS-TC  EllS+SUH ( X ) 

IF(TCELLS*GE*NCELLS)CO  TO  30 
DO  4  J-lflOO 
4  ROHUJ-O 

5  CONTINUE 

30  A¥U-SUH(IX>*4*<  T0+1-E-4 

DO  31  1-1. IX 

SUHX-SUMX+ ( SUM i I » » ( I -1 n 

31  SUM(I>«0 

FM  VU-  ( <  256*  *M YO*YO  t ) *XLEN**2« ) /10*  **8* ) *FLOAT< SUMX ) 

BVU-FHVU/TWGT 

DO  35  Z-1.100 

BMVU=>BMVU+<ROV<Z>*<  «Z-60»**2» ) 

35  CONTINUE 

BMVU-< (4**(Y0*Y0) l/(100**XLEN) )*«BMVU/FLQAT<TCELLS>  I 
GMVU=BMVU+(BVUHXLEN-(FLCG*XLEM)  J ) 
DVU-XLEN*(FL0AT<X)/100* ) 

TCELLS-0 

SUMX-0 

C  COMPUTE  BOV  DOWN  ATTITUDE* 

DO  50  X-100.1.-1 

IX-X 

SX-X-1 

ALOC- ( << SX*I SAYE3 1 ♦ I SA VEI  «*2 ) /I SAYE2 

DO  60  Z-1.50 

SZ-Z-1 

CLOC-< ( <  <SZ*JSAVB3>*JSAVE>-Y0>«*2J/JSAVE2 

SAVLOC-ALOC+CLOC 

DO  70  Y-1.50 

SY-Y-1 

BLOC- ( < ( <  SY* J  SA VE3  >  *JSA VE ) - YO  »  **2 1 /J  SA YS2 
IF(SAVL0C+BL0C*GT*1*»G0  TO  70 
SUM(X»*SUM<X)+4 
ROW! Z >=ROW( Z)+2 
ROV< 100-Z )-ROM( 100-Z »+2 
70  CONTINUE 
60  CONTINUE 

TC  ELLS-TC  ELLS+SUM ( X ) 

IF(TCELLS*G£*NCELLS)GO  TO  90 
DO  65  J-1.100 

65  ROW( J)-0 

60  CONTINUE 

90  AVD-SUM<IX)*4**Y0*Y0*l*X-4 

DO  95  I-100.IX.-1 
SUHX-SUMX+ISUM < I ) *< 100*1 > I 
96  SUM(I)-0 

FM YD- ( ( 266 • *  C YQ*YO  >  *XLEM**2  * ) /1 0*  **0 • > *FLOAT I SUMX ) 
BVD-FMYD/TVGT 
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DO  100  2-1.100 

100  BHVD-BMVD+4R0W4Z)*44Z-50>**2I> 

BMVD-<  44**<YO*YOI)/4100**XLSN)  )*<BMVD/FL0AT4TCSLI.S> > 

CMV D-BH VD+4 BVD" 4  XLEN+FLCG ) ) 

DVD-XLBN* ( FLOAT (100-X) /tOO* > 

TCELLS-0 

SUMX-0 

C  FIND  AREAS  IN  HORIZONTAL  ATTITUDE* 

DO  110  X-1.100 

IX=X 

SX-X-1 

CLOC-f 4 ( ISX*JSAVE3)*JSAf£>-Y0)**2>/JSAVE2 

DO  120  Z-1.100 

SZ-Z-1 

ALOC-4 4  4SZ*ISAVE3)+ISAV£)**2)/ISAYB2 

SAVLOC-AlOC+CLOC 

DO  130  Y-1.50 

SY-Y-1 

BLOC=( ( I <SY*JSAVE3I*JSAFE)-Y0)**2)/JSAVE2 
IF(SAVLOC+BLOC*GT*l* )G0  TO  130 
SUM4XI-SUM4XJ+2 
R0W4YI-R0¥4Y)+2 
130  CONTINUE 
120  CONTINUE 

TCELLS-TCELLS+SUM4XI 
IF(TCELLS*GE*NCELLS)COTO  140 
DO  150  J-1,50 
ROW ( J ) -0 

150  CONTINUE 

110  CONTINUE 

140  AVH-SUM4IX)*2**YO*XLEM*l*B-4 
DO  190  N=1 . IX 
SUHX-SUHX+ 4  SUM  4  N ) *  4 H-l » > 

190  SUH(N>-0* 

DO  205  Y-1.50 

205  BH VH-BM VH+4  ROM 4  Y  •  *  4 Y”50  )  *+2 ) 

BMVH-4 42  **Y0  I /100* ) *4  BMVH/FL0AT4  TCELLS  I I 
FM VH- 4  4 512 • *4 Y0**3*  >  *XLEN ) /1 0**8 » *FLOAT 4  SUMX  > 

BVH=FMVH/TWGT 

DVH-42**Y0)*4FL0AT<X)/100* ) 

700  CONTINUE 

IF(GMVD*CT*CHVH*AND*GNVD*GT*GHVU)TYPE  300. IMESS41.il. 1-1.3) 
IF4GMVU*GT*GMVD.AND*GMVU*GT*GMVH> TYPE  300. 4 MESS42.I I • 1-1 .3) 
IF<CMVH*GE*GMVD*AND*GMVH.GE*CMVU)TYPE  300.4MESS43.D.I-1.3) 
IF4VCG*CT*0* •AND*GMVH*GT*CMVD*AND*GHVH*GT*GMYU)TTPB  400 
300  FORMAT (//*  THE  CAPSULE  IS  FLOATING*  »3A6. /) 

'  400  FORMAT! *  AND  HAS  TURNED  UPSIDE  DOVE**//) 

ISYM-**' 
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CTYPB  302 » I SKA 

'302  FORMAT <  *  SURFACE  SEA  STATE  IS  M2./I 

GSQT-2*/ (32.17*32.17) 

A«« 0081 *32 • 1 7*32* 1 ? 

B«. 74* ( 32 • 17/ ( H <1 SEA+1 » *2 .44 >) **4- 
IF(GMVD»GT*GKVH* AND»GMVD.(JT.QMVU)GO  TO  SOO 
IF<GHVU.OT.OMVH.AMD.OMVU.OT.<JH?D»CO  TO  810 
If (GMVH.CE.GMVO.AND.OMVR.CE.CNVU)OO  TO  820 
600  DRAFT-DVD 

SLOPE-0. 

C-64.*AVD 

B1 - 1 • E-2* SORT  <  4  *  *P I *C ) * DRAFT 
XM-TWGT/32‘2 
Al— 5*XM 
CO  TO  530 
610  DRAFT-DVU 

SLOPE-0. 

C-64»*AVU 

B1 -1  *  E-2*SQRT«4«*PI*C ) *DRAFT 

XM-TWGT/32*2 

Al-XM/2* 

CO  TO  530 
820  DRAFT-0. 

SLOPE-l • 

XM-0- 

C-64**AVH 

Bl*l .E-2*SQRT(4.*PI*C»*DBAFT 
Al-1* 

630  SSLOPE-O* 

SHAMP-O. 

HBAR-O. 

DO  11  1-1,100 
FREQ=(FLOAT< I 1/2.5I/10* 

S 1 1 )-( A/FREQ**5) *EXP< ~B/FREQ**4 ) 

HBAR-HBAR+SI I l*«04 

SPR ( I ) =GSQT*S ( I) *FREQ**4 

S<  I  )-S(  1 1*(  f  (C-Al*FRBQ*FREQ)**2-*-Bl*Bl*FRSQ*PREQ>/ 

*  <  (C-(XM+A1 l*FREQ*FBE0J**2+Bl*Bl*FRE0*FRE0l )* 

*  EXP ( -2  *  *FR  EQ*FREQ*DRAPT/32  *2 1 
SPR  < I ) -SPR  < I ) *SLOPE 
SSL0P£=SSL0PE+.04*FRBQ**4.*S< I)*2«15 
SHAHP-SHAMP+S< 1 1*. 04 
U  Vm-FREQ 

SSL0PE-SQRT(SSL0PE»*57.3*2.5/«2.*32.17) 

SHAMP-4. *SORT ( SHARP ) 

HBAR-4«*SQRT(HBAR> 

. .PE  303.HBAR, SHARP, SSLOPE 

303  FORMAT ( *  SIGHIFICAMT  WAVE  HEIGHT  -  SFIO.S,*  FT*,/ 

*  •  SIGHIFICAMT  HEAVE  AMPL*  -  FT*./ 

*  *  AVERAGE  PITCH  AMPL*  -  *,F10.3»»  DEGREES’,/) 
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IF( SLOPE* EQ*0* ) SPR( 1 )»1 ♦ 

CALL  TGSCALE(W«100t6*»WMINtWSCALEtO) 

CALL  TCSCALE(Stl00«6«  tSNINtSSCALEtO) 

CALL  TCSCALE(SPRtl00*6* •SPRMIKtSPRSCLiOl 
CALL  TGPL0T(7.,0) 

CALL  TCAXIStO* iO*  >6* • SHIN tSSCALEtl* 10*2 • HAMS) 

CALL  TGAXISJO*  *0* >6* i WHIN* VSCALE»4»10tZtKAMR) 

CALL  TGCURVEt VtStlOOtUMINtSMIKf  tfSCALEtSSCALEt ISYMtl ) 

CALL  TGEND 

CALL  TGPL0T(7..0) 

CALL  TGAXIS10* »0* >6* •S?RKIR»SPRSCL>ltlO»2»KAMSP) 

CALL  TGAXISIO*  «0*  »6*  »VMIK» VSCALS«4tlO«2tRAMW) 

CALL  TGCURV£(W»SPRilOO»WMIMtSP&MIRtVSCALS»S?RSCL> 
♦ISYM.1J 

CALL  TCEKD 

RETURN 

END 


IX- 15  Data  Files 


The  messages  accessed  in  response  to  the  monitor 
command  HELP  reside  in  a  data  file  called  HELP.DAT 

Permanent  storage  of  weight  distributions  and  volume 
distributions  is  provided  by  the  data  files  WATE.DAT  and 
V0LUM.DAT,  respectively.  Since  these  files  must  contain 
some  "data"  in  order  to  exist,  an  identification  name, 
A0001,  and  a  single  null  triplet  have  been  stored  in  each 
file . 


These  data  files  are  called  automatically  during 
loading  (provided  they  exist)  and  therefore  are  not  in¬ 
cluded  in  the  LOAD  or  EXECUTE  command  string. 

A  listing  of  the  contants  of  all  three  data  files  is 
shown  on  the  following  pages. 


**  ABBREVIATED  USER  INSTRUCTIONS; 

t 

THE  PROGRAM  YOU  ARE  RUNNING  HILL  HELP  YOU  DESIGN  AN  ELLIPSOIDAL  SUB¬ 
MARINE  WITH  A  SEPARABLE  CAPSULE*  EACH  COMMAND  BELOH  EITHER  PERFORMS  A 
COMPLETE  FUNCTION.  PROMPTS  ADDITIONAL  QUESTIONS  FROM  THE  PROCRAMt  OR  PRO 

VIDES  YOU  NITH  INFORMATION*  THE  PROGRAM  SIGNALS  THAT  IT  IS  READY  FOR  A 
COMMAND  BY  PRINTING  **  AT  THE  LEFT  MARGIN  OR  FOR  A  RESPONSE  BY 
PRINTING  —  * 


COMMAND 


DESCRIPTION 


PRINTS  USER  INSTRUCTIONS 


INPUT  GEOMETRY 


LETS  USER  SPECIFY  GENERAL  CEOMETRIC  CHARAC¬ 
TERISTICS  OF  VEHICLE 


INPUT  WEIGHTS 


LETS  USER  SPECIFY  HEIGHT  DISTRIBUTION  OF  VEHICLE 


INPUT  VOLUMES 


LETS  USER  SPECIFY  FLOODABLE  VOLUME  DISTRIBUTION 


INPUT  DEPTH 


LETS  USER  SPECIFY  INITIAL  DEPTH 


INPUT  SPEED 


LETS  USER  SPECIFY  INITIAL  VERTICAL  AND  HORIZONTAL 
VELOCITIES 


INPUT  TRIM 
* 

INPUT  SEAS 

t 

UPDATE 


LETS  USER  SPECIFY  INITIAL  TRIM  ANGLE 

LETS  USER  SPECIFY  SURFACE  SEA  STATS 

ALLOWS  CHANGES  TO  BE  MADE  TO  HEIGHT  OR  VOLUME 
DISTRIBUTIONS 


WEIGHTS 

* 

DAMAGE 


CALCULATES  HEIGHTS  AND  CENTERS 


SIMULATES  SELECTIVE  FLOODING  OF  VOLUMES 


CALCULATES  BOUYANT  FORCE  AND  TRIM  OF  INTACT 
VEHICLE 


SEPARATE 


SIMULATES  CAPSULE  SEPARATION  FROM  HULL 


ASCEND 


CALCULATES  CAPSULE  TRAJECTORY 


222 


HULL 

I 

.  GRAPH 


t 

MOTION 

$ 

PRINT 

I 

CHECK 

* 

NEW  CASS 

S 

t 


CALCULATES  REMAIN  INC  HULL  TRAJECTORY 

PLOTS  BEHAVIOR  Of  SEPARATED  CAPSULE  OR  REMAINING 
VEHICLE  (WHICHEVER  HAD  THE  MORE  RECENT  TRAJECTORY 
COMPUTED > 

DETERMINES  THE  ATTITUDE  OF  CAPSULE  ON 

THE  SURFACE  AND  PLOTS  THE  SEA  SURFACE  SPECTRUM 

PRINTS  TRAJECTORY  DATA  IN  TABULAR  FORM 

PRINTS  CURRENT  STATUS  OF  SIMULATION 

INITIALIZES  PROGRAM  FOR  A  NEW  SIMULATION 
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